The evolution of the azoospermia factor region AZFa in higher primates

Clones of a PAC contig encompassing the human AZFa region in Yq11.21 were comparatively FISH mapped to great ape Y chromosomes. While the orthologous AZFa locus in the chimpanzee, the bonobo and the gorilla maps to the long arm of their Y chromosomes in Yq12.1-->q12.2, Yq13.1-->q13.2 and Yq11.2, respectively, it is found on the short arm of the orang-utan subspecies of Borneo and Sumatra, in Yp12.3 and Yp13.2, respectively. Regarding the order of PAC clones and genes within the AZFa region, no differences could be detected between apes and man, indicating a strong evolutionary stability of this non-recombining region.

Evolutionary conserved chromosomal segments in the human karyotype are bounded by unstable chromosome bands

In this paper an ancestral karyotype for primates, defining for the first time the ancestral chromosome morphology and the banding patterns, is proposed, and the ancestral syntenic chromosomal segments are identified in the human karyotype. The chromosomal bands that are boundaries of ancestral segments are identified. We have analyzed from data published in the literature 35 different primate species from 19 genera, using the order Scandentia, as well as other published mammalian species as out-groups, and propose an ancestral chromosome number of 2n = 54 for primates, which includes the following chromosomal forms: 1(a+c(1)), 1(b+c(2)), 2a, 2b, 3/21, 4, 5, 6, 7a, 7b, 8, 9, 10a, 10b, 11, 12a/22a, 12b/22b, 13, 14/15, 16a, 16b, 17, 18, 19a, 19b, 20 and X and Y. From this analysis, we have been able to point out the human chromosome bands more "prone" to breakage during the evolutionary pathways and/or pathology processes. We have observed that 89.09% of the human chromosome bands, which are boundaries for ancestral chromosome segments, contain common fragile sites and/or intrachromosomal telomeric-like sequences. A more in depth analysis of twelve different human chromosomes has allowed us to determine that 62.16% of the chromosomal bands implicated in inversions and 100% involved in fusions/fissions correspond to fragile sites, intrachromosomal telomeric-like sequences and/or bands significantly affected by X irradiation. In addition, 73% of the bands affected in pathological processes are co-localized in bands where fragile sites, intrachromosomal telomeric-like sequences, bands significantly affected by X irradiation and/or evolutionary chromosomal bands have been described. Our data also support the hypothesis that chromosomal breakages detected in pathological processes are not randomly distributed along the chromosomes, but rather concentrate in those important evolutionary chromosome bands which correspond to fragile sites and/or intrachromosomal telomeric-like sequences.

New insights into the evolution of chromosome 1

A complex low-repetitive human DNA probe (BAC RP11-35B4) together with two microdissection-derived region-specific probes of the multicolor banding (MCB) probe-set for chromosome 1 were used to re-analyze the evolution of human chromosome 1 in comparison to four ape species. BAC RP11-35B4 derives from 1q21 and contains 143 kb of non-repetitive DNA; however, it produces three specific FISH signals in 1q21, 1p12 and 1p36.1 of Homo sapiens (HSA). Human chromosome 1 was studied in comparison to its homologues in Hylobates lar (HLA), Pongo pygmaeus (PPY), Gorilla gorilla (GGO) and Pan troglodytes (PTR). A duplication of sequences homologous to human 1p36.1 could be detected in PPY plus an additional signal on PPY 16q. The region homologous to HSA 1p36.1 is also duplicated in HLA, and split onto chromosomes 7q and 9p; the region homologous to HSA 1q21/1p12 is present as one region on 5q. Additionally, the breakpoint of a small pericentric inversion in the evolution of human chromosome 1 compared to other great ape species could be refined. In summary, the results obtained here are in concordance with previous reports; however, there is evidence for a deletion of regions homologous to human 1p34.2-->p34.1 during evolution in the Pongidae branch after separation of PPY.

Molecular evolution of the human chromosome 15 pericentromeric region

We present a detailed molecular evolutionary analysis of 1.2 Mb from the pericentromeric region of human 15q11. Sequence analysis indicates the region has been subject to extensive interchromosomal and intrachromosomal duplications during primate evolution. Comparative FISH analyses among non-human primates show remarkable quantitative and qualitative differences in the organization and duplication history of this region - including lineage-specific deletions and duplication expansions. Phylogenetic and comparative analyses reveal that the region is composed of at least 24 distinct segmental duplications or duplicons that have populated the pericentromeric regions of the human genome over the last 40 million years of human evolution. The value of combining both cytogenetic and experimental data in understanding the complex forces which have shaped these regions is discussed.

Population differences in DNA sequence variation and linkage disequilibrium at the PON1 gene

Polymorphisms of the promoter region (-108C/T) and the coding region (192Q/R) of the paraoxonase 1 gene (PON1) showed differences in association with cardiovascular disease risk in various populations. To characterize the genetic variation underlying these important polymorphisms, we examined DNA sequence variation both in a 1.3-kb promoter region 16.5 kb from codon 192, and in a 1.7-kb region centered on the 192Q/R polymorphic site of the coding region of PON1, in 30 Africans, 30 Europeans and 64 Japanese. We found 10 polymorphic sites and 11 haplotypes in the 1.3-kb promoter region and 10 biallelic polymorphic sites and 10 haplotypes in the 1.7-kb region. From the PON1 sequences of chimpanzees and an orangutan, the ancestral type of codon 192 was found to be R. The number of pairs of polymorphic sites between the promoter and 1.7-kb regions that were in significant linkage disequilibrium was much higher in a Japanese population than in African and European populations. In addition, the pairs of polymorphic sites in linkage disequilibrium differed among the three populations. These results suggest that some of the population differences in association with risk for coronary heart disease can be explained by population differences in haplotype frequency of PON1 haplotypes.

Evolution of beta satellite DNA sequences: evidence for duplication-mediated repeat amplification and spreading

In this article, we report studies on the evolutionary history of beta satellite repeats (BSR) in primates. In the orangutan genome, the bulk of BSR sequences was found organized as very short stretches of approximately 100 to 170 bp, embedded in a 60-kb to 80-kb duplicated DNA segment. The estimated copy number of the duplicon that carries BSR sequences ranges from 70 to 100 per orangutan haploid genome. In both macaque and gibbon, the duplicon mapped to a single chromosomal region at the boundary of the rDNA on the marker chromosome (chromosome 13 and 12, respectively). However, only in the gibbon, the duplicon comprised 100 bp of beta satellite. Thus, the ancestral copy of the duplicon appeared in Old World monkeys ( approximately 25 to approximately 35 MYA), whereas the prototype of beta satellite repeats took place in a gibbon ancestor, after apes/Old World monkeys divergence ( approximately 25 MYA). Subsequently, a burst in spreading of the duplicon that carries the beta satellite was observed in the orangutan, after lesser apes divergence from the great apes-humans lineage ( approximately 18 MYA). The analysis of the orangutan genome also indicated the existence of two variants of the duplication that differ for the length (100 or 170 bp) of beta satellite repeats. The latter organization was probably generated by nonhomologous recombination between two 100-bp repeated regions, and it likely led to the duplication of the single Sau3A site present in the 100-bp variant, which generated the prototype of Sau3A 68-bp beta satellite tandem organization. The two variants of the duplication, although with a different ratios, characterize the hominoid genomes from the orangutan to humans, preferentially involving acrocentric chromosomes. At variance to alpha satellite, which appeared before the divergence of New World and Old World monkeys, the beta satellite evolutionary history began in apes ancestor, where we have first documented a low-copy, nonduplicated BSR sequence. The first step of BSR amplification and spreading occurred, most likely, because the BSR was part of a large duplicon, which underwent a burst dispersal in great apes' ancestor after the lesser apes' branching. Then, after orangutan divergence, BSR acquired the clustered structural organization typical of satellite DNA.

Myosin gene mutation correlates with anatomical changes in the human lineage

Powerful masticatory muscles are found in most primates, including chimpanzees and gorillas, and were part of a prominent adaptation of Australopithecus and Paranthropus, extinct genera of the family Hominidae. In contrast, masticatory muscles are considerably smaller in both modern and fossil members of Homo. The evolving hominid masticatory apparatus--traceable to a Late Miocene, chimpanzee-like morphology--shifted towards a pattern of gracilization nearly simultaneously with accelerated encephalization in early Homo. Here, we show that the gene encoding the predominant myosin heavy chain (MYH) expressed in these muscles was inactivated by a frameshifting mutation after the lineages leading to humans and chimpanzees diverged. Loss of this protein isoform is associated with marked size reductions in individual muscle fibres and entire masticatory muscles. Using the coding sequence for the myosin rod domains as a molecular clock, we estimate that this mutation appeared approximately 2.4 million years ago, predating the appearance of modern human body size and emigration of Homo from Africa. This represents the first proteomic distinction between humans and chimpanzees that can be correlated with a traceable anatomic imprint in the fossil record.

Human-specific amino acid changes found in 103 protein-coding genes

We humans have many characteristics that are different from those of the great apes. These human-specific characters must have arisen through mutations accumulated in the genome of our direct ancestor after the divergence of the last common ancestor with chimpanzee. Gene trees of human and great apes are necessary for extracting these human-specific genetic changes. We conducted a systematic analysis of 103 protein-coding genes for human, chimpanzee, gorilla, and orangutan. Nucleotide sequences for 18 genes were newly determined for this study, and those for the remaining genes were retrieved from the DDBJ/EMBL/GenBank database. The total number of amino acid changes in the human lineage was 147 for 26,199 codons (0.56%). The total number of amino acid changes in the human genome was, thus, estimated to be about 80,000. We applied the acceleration index test and Fisher's synonymous/nonsynonymous exact test for each gene tree to detect any human-specific enhancement of amino acid changes compared with ape branches. Six and two genes were shown to have significantly higher nonsynonymous changes at the human lineage from the acceleration index and exact tests, respectively. We also compared the distribution of the differences of the nonsynonymous substitutions on the human lineage and those on the great ape lineage. Two genes were more conserved in the ape lineage, whereas one gene was more conserved in the human lineage. These results suggest that a small proportion of protein-coding genes started to evolve differently in the human lineage after it diverged from the ape lineage.

The SPANX gene family of cancer/testis-specific antigens: rapid evolution and amplification in African great apes and hominids

Human sperm protein associated with the nucleus on the X chromosome (SPANX) genes comprise a gene family with five known members (SPANX-A1, -A2, -B, -C, and -D), encoding cancer/testis-specific antigens that are potential targets for cancer immunotherapy. These highly similar paralogous genes cluster on the X chromosome at Xq27. We isolated and sequenced primate genomic clones homologous to human SPANX. Analysis of these clones and search of the human genome sequence revealed an uncharacterized group of genes, SPANX-N, which are present in all primates as well as in mouse and rat. In humans, four SPANX-N genes comprise a series of tandem duplicates at Xq27; a fifth member of this subfamily is located at Xp11. Similarly to SPANX-A/D, human SPANX-N genes are expressed in normal testis and some melanoma cell lines; testis-specific expression of SPANX is also conserved in mouse. Analysis of the taxonomic distribution of the long and short forms of the intron indicates that SPANX-N is the ancestral form, from which the SPANX-A/D subfamily evolved in the common ancestor of the hominoid lineage. Strikingly, the coding sequences of the SPANX genes evolved much faster than the intron and the 5' untranslated region. There is a strong correlation between the rates of evolution of synonymous and nonsynonymous codon positions, both of which are accelerated 2-fold or more compared to the noncoding sequences. Thus, evolution of the SPANX family appears to have involved positive selection that affected not only the protein sequence but also the synonymous sites in the coding sequence.

Amino-acid changes acquired during evolution by olfactory receptor 912-93 modify the specificity of odorant recognition

The sense of smell in mammals can perceive and discriminate a wide variety of volatile odorants. Odorants bind to specific olfactory receptors (ORs) to initiate an action potential that transduces olfactory information to the olfactory cortex. We previously identified the structural motifs of odorant molecules (aliphatic 2- or 3-ketones) required to activate mouse OR912-93 by detection of the odorant response using calcium measurement in transfected cells. In order to study changes in the specificity of this receptor that might have occurred during evolution, we cloned the orthologous genes from six primate species and pig and assayed the encoded receptors for responses to odorants. Primate OR912-93 orthologs share 88-97% sequence identity. All the receptors responded to 2- and 3-heptanone except the squirrel-monkey OR, which responded only to 3-heptanone, and the human and orangutan ORs, which were not functional. Directed mutagenesis allowed us to convert the squirrel-monkey response to that of the other functional 912-93 ORs by substituting three amino acids in the second extracellular loop. Orangutan and human 912-93 ORs regained function after restoration of the arginine residue in the DRY motif required for G-protein activation. However, the human receptor was constitutively activated in the absence of ligand stimulation. Using natural mutants of the OR912-93 receptor, we provide evidence that squirrel-monkeys evolved towards a restriction of the specificity of this receptor and therefore that slight alterations in the sequence of a receptor can induce subtle changes in recognition specificity.

High mutation rates in human and ape pseudoautosomal genes

It has been suggested that recombination may be mutagenic, which, if true, would inflate intraspecies diversity and interspecies silent divergence in regions of high recombination. Here, we test this hypothesis comparing human/orangutan genome-wide non-coding divergence (K) to that in the pseudoautosomal genes which were reported to recombine much more frequently than the rest of the genome. We demonstrate that, compared to the average human/orangutan non-coding divergence (K=3%), the substitution rate is significantly elevated in the introns of SHOX (K=5.7%), PPP2R3L (K=8.7%) and ASMT (K=6.5%) genes located in the human and orangutan Xp/Yp pseudoautosomal region (p-PAR), where recombination is over 20-fold higher than the genomic average. On the other hand, human/orangutan non-coding divergence at the Xp/Yp pseudoautosomal boundary (K=3.5%) and in the SYBL1 gene (K=2.7%), located in the human Xq/Yq pseudoautosomal region (q-PAR), where recombination is known to be less frequent than in p-PAR, was not significantly higher than the genome average. The data are consistent with the hypothesis that recombination may be mutagenic.

Recombination has little effect on the rate of sequence divergence in pseudoautosomal boundary 1 among humans and great apes

Recent studies indicated that recombination is strongly mutagenic. In particular, data from the mouse pseudoautosomal boundary (PAB) suggested that locally intensive recombination increased the nucleotide substitution rate by more than 100-fold and greatly increased the GC content. Here we study the rates of nucleotide substitution in eight introns of the human and great ape XG gene, which spans the boundary between the pseudoautosomal region 1 (PAR1) and the X-specific region. Contrary to what is expected under the above hypothesis, our sequence data from humans and great apes reveal that the PAR1 introns of XG have actually evolved slightly slower than X-specific introns. Only when a New World monkey was compared with hominoids were the rates slightly increased in the PAR1 introns. In terms of base composition, although the intergenic regions of the human PAR1 show a significant increase of G and C nucleotides, the base composition of the surveyed PAR1 introns is similar to that of the X-specific introns. Direct and indirect evidence indicates that the recombination rate is, indeed, much higher in PAR1 introns than in X-specific introns, and that the present PAB has persisted since the common ancestor of hominoids. Therefore, the mutagenic effect of recombination is far weaker than previously proposed, at least in hominoid PABs.

Mixed-Model Reanalysis of Primate Data Suggests Tissue and Species Biases in Oligonucleotide-Based Gene Expression Profiles

An emerging issue in evolutionary genetics is whether it is possible to use gene expression profiling to identify genes that are associated with morphological, physiological, or behavioral divergence between species and whether these genes have undergone positive selection. Some of these questions were addressed in a recent study (Enard  et al. 2002) of the difference in gene expression among human, chimp, and orangutan, which suggested an accelerated rate of divergence in gene expression in the human brain relative to liver. Reanalysis of the Affymetrix data set using analysis of variance methods to quantify the contributions of individuals and species to variation in expression of 12,600 genes indicates that as much as one-quarter of the genome shows divergent expression between primate species at the 5% level. The magnitude of fold change ranges from 1.2-fold up to 8-fold. Similar conclusions apply to reanalysis of Enard  et al.'s (2002) parallel murine data set. However, biases inherent to short oligonucleotide microarray technology may account for some of the tissue and species effects. At high significance levels, more differences were observed in the liver than in the brain in each of the pairwise species comparisons, so it is not clear that expression divergence is accelerated in the human brain. Further, there is an apparent bias toward upregulation of gene expression in the brain in both primates and mice, whereas genes are equally likely to be up- or downregulated in the liver when these species diverge. A small subset of genes that are candidates for adaptive divergence may be identified on the basis of a high ratio of interspecific to intraspecific divergence.

Chromosomal phylogeny and evolution of gibbons (Hylobatidae)

Although human and gibbons are classified in the same primate superfamily (Hominoidae), their karyotypes differ by extensive chromosome reshuffling. To date, there is still limited understanding of the events that shaped extant gibbon karyotypes. Further, the phylogeny and evolution of the twelve or more extant gibbon species (lesser apes, Hylobatidae) is poorly understood, and conflicting phylogenies have been published. We present a comprehensive analysis of gibbon chromosome rearrangements and a phylogenetic reconstruction of the four recognized subgenera based on molecular cytogenetics data. We have used two different approaches to interpret our data: (1) a cladistic reconstruction based on the identification of ancestral versus derived chromosome forms observed in extant gibbon species; (2) an approach in which adjacent homologous segments that have been changed by translocations and intra-chromosomal rearrangements are treated as discrete characters in a parsimony analysis (PAUP). The orangutan serves as an "outgroup", since it has a karyotype that is supposed to be most similar to the ancestral form of all humans and apes. Both approaches place the subgenus Bunopithecus as the most basal group of the Hylobatidae, followed by Hylobates, with Symphalangus and Nomascus as the last to diverge. Since most chromosome rearrangements observed in gibbons are either ancestral to all four subgenera or specific for individual species and only a few common derived rearrangements at subsequent branching points have been recorded, all extant gibbons may have diverged within relatively short evolutionary time. In general, chromosomal rearrangements produce changes that should be considered as unique landmarks at the divergence nodes. Thus, molecular cytogenetics could be an important tool to elucidate phylogenies in other species in which speciation may have occurred over very short evolutionary time with not enough genetic (DNA sequence) and other biological divergence to be picked up.

Gene diversity patterns at 10 X-chromosomal loci in humans and chimpanzees

We have investigated the pattern and extent of nucleotide diversity in 10 X-chromosomal genes where mutations are known to cause mental retardation in humans. For each gene, we sequenced the entire coding region from cDNA in humans, chimpanzees, and orangutans, as well as about 3 kb of genomic DNA in 20 humans sampled worldwide and in 10 chimpanzees representing two "subspecies." Overall nucleotide diversity in these genes is about twofold lower in humans than in chimpanzees, and nucleotide diversity within and between species is low, suggesting that a high level of functional constraint acts on these genes. Strikingly, we find that a summary of the allele frequency spectrum is significantly correlated in humans and chimpanzees, perhaps reflecting very similar levels of constraint at these genes in the two species. A possible exception is FMR2, which shows a higher number of nonsynonymous than synonymous substitutions on the human lineage, suggesting the action of positive selection.

Divergence of the genes on human chromosome 21 between human and other hominoids and variation of substitution rates among transcription units

The study of genomic divergence between humans and primates may provide insight into the origins of human beings and the genetic basis of unique human traits and diseases. Chromosome 21 is the smallest chromosome in the human genome, and some of its regions have been implicated in mental retardation and other diseases. In this study, we sequenced the coding and regulatory regions of 127 known genes on human chromosome 21 in DNA samples from human and chimpanzees and a part of the corresponding genes from orangutan, gorilla, and macaque. Overall, 3,003 nucleotide differences between human and chimpanzee were identified over approximately 400 kb. The differences in coding, promoter, and exon-intron junction regions were 0.51 +/- 0.02%, 0.88 +/- 0.03%, and 0.85 +/- 0.02%, respectively, much lower than the previously reported 1.23% in genomic regions, which suggests the presence of purifying selection. Significant variation in substitution rate among genes was observed by comparing the divergence between human and chimpanzee. Furthermore, by implementing a bioinformatics-based approach, we showed that the identification of genetic variants specific to the human lineage might lead to an understanding of the mechanisms that are attributable to the phenotypes that unique to humans, by changing the structure and/or dosage of the proteins expressed. A phylogenetic analysis unambiguously confirms the conclusion that chimpanzees were our closest relatives to the exclusion of other primates and the relative divergence of the Homo-Pan and that of (Homo-Pan)-Gorilla are 4.93 million years and 7.26 million years, respectively.

Neutral substitutions occur at a faster rate in exons than in noncoding DNA in primate genomes

Point mutation rates in exons (synonymous sites) and noncoding (introns and intergenic) regions are generally assumed to be the same. However, comparative sequence analyses of synonymous substitutions in exons (81 genes) and that of long intergenic fragments (141.3 kbp) of human and chimpanzee genomes reveal a 30%-60% higher mutation rate in exons than in noncoding DNA. We propose a differential CpG content hypothesis to explain this fundamental, and seemingly unintuitive, pattern. We find that the increased exonic rate is the result of the relative overabundance of synonymous sites involved in CpG dinucleotides, as the evolutionary divergence in non-CpG sites is similar in noncoding DNA and synonymous sites of exons. Expectations and predictions of our hypothesis are confirmed in comparisons involving more distantly related species, including human-orangutan, human-baboon, and human-macaque. Our results suggest an underlying mechanism for higher mutation rate in GC-rich genomic regions, predict nonlinear accumulation of mutations in pseudogenes over time, and provide a possible explanation for the observed higher diversity of single nucleotide polymorphisms (SNPs) in the synonymous sites of exons compared to the noncoding regions.

Inactivation of the EGF-TM7 receptor EMR4 after the Pan-Homo divergence

We here report on the identification of a novel human EGF-TM7 receptor, designated EMR4. Like most EGF-TM7 receptor genes, EMR4 is localized on the short arm of chromosome 19, in close proximity to EMR1. Remarkably, due to a one-nucleotide deletion in exon 8, translation of human EMR4 would result in a truncated 232-amino acid protein lacking the entire seven-span transmembrane region. This deletion is not present in nonhuman primates, including chimpanzees, suggesting that EMR4 became nonfunctional only after human speciation, about five million years ago. Thus, EMR4 surprisingly accounts for a genetic difference between humans and primates related to immunity.

Human specific loss of olfactory receptor genes

Olfactory receptor (OR) genes constitute the basis for the sense of smell and are encoded by the largest mammalian gene superfamily of >1,000 genes. In humans, >60% of these are pseudogenes. In contrast, the mouse OR repertoire, although of roughly equal size, contains only approximately 20% pseudogenes. We asked whether the high fraction of nonfunctional OR genes is specific to humans or is a common feature of all primates. To this end, we have compared the sequences of 50 human OR coding regions, regardless of their functional annotations, to those of their putative orthologs in chimpanzees, gorillas, orangutans, and rhesus macaques. We found that humans have accumulated mutations that disrupt OR coding regions roughly 4-fold faster than any other species sampled. As a consequence, the fraction of OR pseudogenes in humans is almost twice as high as in the non-human primates, suggesting a human-specific process of OR gene disruption, likely due to a reduced chemosensory dependence relative to apes.

Genomic DNA insertions and deletions occur frequently between humans and nonhuman primates

Comparative DNA sequence studies between humans and nonhuman primates will be important for understanding the genetic basis of the phenotypic differences between these species. Here we compare approximately 27 Mb of human chromosome 21 with chimpanzee DNA sequences identifying 57 genomic rearrangements (deletions and insertions ranging in size from 0.2 to 8.0 kb) between the two species. These rearrangements are distributed along the entire length of chromosome 21, with approximately 35% found in genomic intervals encoding genes (genic intervals), and have occurred in the genomes of both humans and chimpanzees. Comparison of approximately 9 Mb of human chromosome 21 with orangutan, rhesus macaque, and woolly monkey DNA sequences identified a combined total of 114 genomic rearrangements between humans and nonhuman primates. Analysis of these rearrangements revealed that they are randomly distributed with respect to genic and nongenic intervals and identified one deletion that has likely resulted in the inactivation of a gene (beta1,3-galactosyltransferase) in the woolly monkey. Our data show that genomic rearrangements have occurred frequently during primate genome evolution and significantly contribute to the DNA differences between these species. These DNA rearrangements are commonly found in genic intervals, and thus provide natural starting points for focused investigations of qualitative and quantitative gene expression differences between humans and other primates.

Genetic sex identification in orangutans

To date, no established protocol for genetic sex identification in orangutans (Pongo pygmaeus) exists. In nearly all apes (gibbons, gorillas, chimpanzees, and humans), genetic sex identification is possible using the amelogenin gene because copies located on X and Y chromosomes have different sizes. Here we report that orangutan sex identification can be resolved through multiplex polymerase chain reaction (PCR) of the Y-linked SRY locus and the amelogenin locus. PCR amplifications of orangutan amelogenin produces one fragment size in both sexes, while SRY amplifies only in males. This protocol will allow primatologists to identify the sex of orangutans through genetic analysis.

Analysis of primate genomic variation reveals a repeat-driven expansion of the human genome

We performed a detailed analysis of both single-nucleotide and large insertion/deletion events based on large-scale comparison of 10.6 Mb of genomic sequence from lemur, baboon, and chimpanzee to human. Using a human genomic reference, optimal global alignments were constructed from large (>50-kb) genomic sequence clones. These alignments were examined for the pattern, frequency, and nature of mutational events. Whereas rates of single-nucleotide substitution remain relatively constant (1-2 x 10(-9) substitutions/site/year), rates of retrotransposition vary radically among different primate lineages. These differences have lead to a 15%-20% expansion of human genome size over the last 50 million years of primate evolution, 90% of it due to new retroposon insertions. Orthologous comparisons with the chimpanzee suggest that the human genome continues to significantly expand due to shifts in retrotransposition activity. Assuming that the primate genome sequence we have sampled is representative, we estimate that human euchromatin has expanded 30 Mb and 550 Mb compared to the primate genomes of chimpanzee and lemur, respectively.

Large-scale variation among human and great ape genomes determined by array comparative genomic hybridization

Large-scale genomic rearrangements are a major force of evolutionary change and the ascertainment of such events between the human and great ape genomes is fundamental to a complete understanding of the genetic history and evolution of our species. Here, we present the results of an evolutionary analysis utilizing array comparative genomic hybridization (array CGH), measuring copy-number gains and losses among these species. Using an array of 2460 human bacterial artificial chromosomes (BACs) (12% of the genome), we identified a total of 63 sites of putative DNA copy-number variation between humans and the great apes (chimpanzee, bonobo, gorilla, and orangutan). Detailed molecular characterization of a subset of these sites confirmed rearrangements ranging from 40 to at least 175 kb in size. Surprisingly, the majority of variant sites differentiating great ape and human genomes were found within interstitial euchromatin. These data suggest that such large-scale events are not restricted solely to subtelomeric or pericentromeric regions, but also occur within genic regions. In addition, 5/9 of the verified variant sites localized to areas of intrachromosomal segmental duplication within the human genome. On the basis of the frequency of duplication in humans, this represents a 14-fold positional bias. In contrast to previous cytogenetic and comparative mapping studies, these results indicate extensive local repatterning of hominoid chromosomes in euchromatic regions through a duplication-driven mechanism of genome evolution.

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.