High diversity in functional properties of melanocortin 1 receptor (MC1R) in divergent primate species is more strongly associated with phylogeny than coat color.

We have characterized the biochemical function of the melanocortin 1 receptor (MC1R), a critical regulator of melanin synthesis, from 9 phylogenetically diverse primate species with varying coat colors. There is substantial diversity in melanocyte-stimulating hormone (MSH) binding affinity and basal levels of activity in the cloned MC1Rs. MSH binding was lost independently in lemur and New World monkey lineages, whereas high basal levels of MC1R activity occur in lemurs and some New World monkeys and Old World monkeys. Highest levels of basal activity were found in the MC1R of ruffed lemurs, which have the E94K mutation that leads to constitutive activation in other species. In 3 species (2 lemurs and the howler monkey), we report the novel finding that binding and inhibition of MC1R by agouti signaling protein (ASIP) can occur when MSH binding has been lost, thus enabling continuing regulation of the melanin type via ASIP expression. Together, these findings can explain the previous paradox of a predominantly pheomelanic coat in the red ruffed lemur (Varecia rubra). The presence of a functional, MSH-responsive MC1R in orangutan demonstrates that the mechanism of red hair generation in this ape is different from the prevalent mechanism in European human populations. Overall, we have found unexpected diversity in MC1R function among primates and show that the evolution of the regulatory control of MC1R activity occurs by independent variation of 3 distinct mechanisms: basal MC1R activity, MSH binding and activation, and ASIP binding and inhibition. This diversity of function is broadly associated with primate phylogeny and does not have a simple relation to coat color phenotype within primate clades.

Chimpanzee, orangutan, mouse, and human cell cycle promoters exempt CCAAT boxes and CHR elements from interspecies differences.

Mechanisms regulating the cell division cycle are well conserved among all eukaryotes. Consistently many proteins regulating the cell cycle are functionally interchangeable between many organisms. Cell division control is regulated on different levels of which the transcriptional level appears to be particularly important for controlling synthesis of many cell cycle proteins. We had earlier described transcription factor-binding sites essential for regulating genes important for the transition from the G(2) phase to mitosis. A tandem repressor site named cell cycle-dependent element (CDE) and cell cycle genes homology region (CHR) are responsible for the correct expression during the cell cycle. Another feature of these G(2)/M-specific promoters is the activation through 2 or 3 CCAAT boxes binding the transcription factor nuclear factor-Y (NF-Y). These major activating sites have to be spaced 32 or 33 bp apart to be fully functional. We were interested in looking at the evolutionary changes in regulatory elements and overall promoter structure of 3 well-characterized cell cycle genes. Here, we compare the DNA sequences and functional features of the cdc25C, cyclin B1, and cyclin B2 promoters from humans, mouse, chimpanzee, and orangutan. We find numerous differences in the nucleotide sequence between mouse and primate promoters. However, CHR and CCAAT boxes stand out in that they are perfectly conserved in all promoters tested. The CDE site contains nucleotide exchanges between mouse and primate promoters. Comparing sequences and functions of chimpanzee, orangutan, and human promoters, we observe a complete conservation in nucleotide sequence of the regulatory elements. Functional assays of the cyclin B1, cyclin B2, and cdc25C promoters yield moderate variations in activity and thereby a good conservation of function. Although we find nucleotide differences in cell cycle promoters between orangutan and humans of about 5%, there are never changes in any of the CCAAT boxes or CDE/CHR sites in the cyclin B1, cyclin B2, and cdc25C promoters. Furthermore, we describe the influence of the tumor suppressor p53 and the transcriptional activator NF-Y on regulation of the newly cloned primate promoters.

Haplotype structure of FSHB, the beta-subunit gene for fertility-associated follicle-stimulating hormone: possible influence of balancing selection.

Follicle-stimulating hormone (FSH) is essential for human reproduction. The unique functions of this hormone are provided by the FSH receptor-binding beta-subunit encoded by the FSHB gene. Resequencing and genotyping of FSHB in three European, two Asian and one African population, as well as in the great apes (chimpanzee, gorilla, orangutan), revealed low diversity and significant excess of polymorphisms with intermediate frequency alleles. Statistical tests for FSHB showed deviations from neutrality in all populations suggesting a possible effect of balancing selection. Two core haplotypes were identified (carried by 76-96.6% of each population's sample), the sequences of which are clearly separated from each other. As fertility most directly affects an organism's fitness, the carriers of these haplotypes have apparently had more success in human history to contribute to the next generation. There is a preliminary observation suggesting that the second most frequent FSHB haplotype may be associated with rapid conception success in females. Interestingly, the same haplotype is related to an ancestral FSHB variant shared with the ancestor of the great apes. The determination of the functional consequence of the two core FSHB variants may have implications for understanding and regulating human fertility, as well as in assisting infertility treatments.

Adaptive selection of mitochondrial complex I subunits during primate radiation.

Mammalian oxidative phosphorylation (OXPHOS) complexes I, III, IV and V are assembled from both mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) encoded subunits, with complex I encompassing 39 nDNA and seven mtDNA subunits. Yet the sequence variation of the mtDNA genes is more than ten fold greater than that of the nDNA encoded genes of the OXPHOS complexes and the mtDNA proteins have been found to be influenced by positive (adaptive) selection. To maintain a functional complex I, nDNA and mtDNA subunits must interact, implying that certain nDNA complex I genes may also have been influenced by positive selection. To determine if positive selection has influenced nDNA complex I genes, we analyzed the DNA sequences of all of the nDNA and mtDNA encoded complex I subunits from orangutan, gorilla, chimpanzee, human and all available vertebrate sequences. This revealed that three nDNA complex I genes (NDUFC2, NDUFA1, and NDUFA4) had significantly increased amino acid substitution rates by both PAML and Z-test, suggesting that they have been subjected to adaptive selection during primate radiation. Since all three of these subunits reside in the membrane domain of complex I along with the mtDNA subunits, we compared amino acid changes in these three nDNA genes with those of the mtDNA genes across species. Changes in the nDNA NDUFC2 cysteine 39 were found to correlate with those in the mtDNA ND5 cysteine 330. Therefore, adaptive selection has influenced some nDNA complex I genes and nDNA and mtDNA complex I genes may have co-evolved.

Expression profiling in primates reveals a rapid evolution of human transcription factors.

Although it has been hypothesized for thirty years that many human adaptations are likely to be due to changes in gene regulation, almost nothing is known about the modes of natural selection acting on regulation in primates. Here we identify a set of genes for which expression is evolving under natural selection. We use a new multi-species complementary DNA array to compare steady-state messenger RNA levels in liver tissues within and between humans, chimpanzees, orangutans and rhesus macaques. Using estimates from a linear mixed model, we identify a set of genes for which expression levels have remained constant across the entire phylogeny (approximately 70 million years), and are therefore likely to be under stabilizing selection. Among the top candidates are five genes with expression levels that have previously been shown to be altered in liver carcinoma. We also find a number of genes with similar expression levels among non-human primates but significantly elevated or reduced expression in the human lineage, features that point to the action of directional selection. Among the gene set with a human-specific increase in expression, there is an excess of transcription factors; the same is not true for genes with increased expression in chimpanzee.

Genomic structure and paralogous regions of the inversion breakpoint occurring between human chromosome 3p12.3 and orangutan chromosome 2.

Intrachromosomal duplications play a significant role in human genome pathology and evolution. To better understand the molecular basis of evolutionary chromosome rearrangements, we performed molecular cytogenetic and sequence analyses of the breakpoint region that distinguishes human chromosome 3p12.3 and orangutan chromosome 2. FISH with region-specific BAC clones demonstrated that the breakpoint-flanking sequences are duplicated intrachromosomally on orangutan 2 and human 3q21 as well as at many pericentromeric and subtelomeric sites throughout the genomes. Breakage and rearrangement of the human 3p12.3-homologous region in the orangutan lineage were associated with a partial loss of duplicated sequences in the breakpoint region. Consistent with our FISH mapping results, computational analysis of the human chromosome 3 genomic sequence revealed three 3p12.3-paralogous sequence blocks on human chromosome 3q21 and smaller blocks on the short arm end 3p26-->p25. This is consistent with the view that sequences from an ancestral site at 3q21 were duplicated at 3p12.3 in a common ancestor of orangutan and humans. Our results show that evolutionary chromosome rearrangements are associated with microduplications and microdeletions, contributing to the DNA differences between closely related species.

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.

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.

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.

Progressive inactivation of the haploid expressed gene for the sperm-specific endozepine-like peptide (ELP) through primate evolution.

The endozepine-like peptide (ELP) is a novel intracellular molecule which is expressed in high amounts at both mRNA and protein levels very specifically in late haploid male germ cells. It is closely related to the ubiquitous acyl-CoA binding protein, is highly conserved, shares a similar ability to bind mid-long chain acyl-CoA, and is thus likely to be involved in mature sperm metabolism. While it has been characterized from diverse mammals, it has so far not been possible to identify an equivalent molecule in the primate testis. Using a PCR approach, combined with cDNA cloning and Northern hybridization, testicular transcripts and/or genomic DNA were analysed for different primate species, including human. In the marmoset and cynomolgus macaque normally structured transcripts appear to be expressed, though at a low level. In the human testis, two rare transcripts were characterized, hELP1 and hELP2, the products of independent duplicated genes. Both transcripts were longer than in non-mammalian species, included frame-shift mutations and substantial sequence insertions, preventing the translation of a sensible protein. Genomic PCR analysis of three anthropoid species, chimpanzee, gorilla and orangutan, showed the presence of a similarly mutated hELP1 gene. Only in the gorilla was a hELP2 gene identified, apparently lacking the frame-shift mutation, and thus potentially able to give rise to a functional ELP protein. Taken together, these results show that during primate evolution there has been a progressive inactivation of the ELP gene, initially with a down-regulation in lower primates, and subsequently with inactivating mutations in the open reading frame. At some time during simian evolution prior to these mutations there has been a gene duplication, though this second gene has also become inactivated in humans. In its pattern of evolution the ELP gene shows similarities with the MDC/fertilin family, whose members are also considered essential components of haploid sperm in non-primates, but which are progressively inactivated in anthropoids and humans. We should like to speculate that the established subfertility of the human male may not be a recent event, but the consequence of a longer evolutionary process whereby primates have traded off absolute fertility against social or sexual advantages.

Evolutionary divergence of the oncogenes GLI, HST and INT2.

Almost a quarter of a century ago, the banding patterns of human and other higher primate chromosomes were compared, creating a barrage of speculation. Consequently, a number of approaches have been used to understand human descent. Chromosome modifications are believed to be important in the origin of species, and pericentric inversions account for the majority of evolutionary chromosomal alterations seen in Hominoidea. A comparative mapping fluorescence in situ hybridization technique, using locus-specific DNA probes as phylogenotic markers, was used to decipher the pericentric inversions of human chromosomes 11 and 12. Human-derived (Homo sapiens, HSA) DNA probes for GLI, HST and INT2 protooncogenes were used to identify their homologous locations in the chromosomes of chimpanzee (Pan troglodytes, PTR), gorilla (Gorilla gorilla, GGO) and orangutan (Pongo pygmaeus, PPY). The INT2 and HST loci mapping results confirm the earlier putative claim that a pericentric inversion took place in HSA chromosome 11 and its equivalent PTR and GGO chromosomes. In addition, these data provide additional information regarding the orangutan's position on the evolutionary tree of Pongidae and Hominidae. GLI mapping reveals that a pericentric inversion occurred in the HSA chromosome 12 equivalent in PTR and GGO, but was not seen in HSA or PPY. These pericentric inversions in PTR and GGO may have occurred at a period when both PTR and GGO had branched off from the Hominoidae trunk. The use of loci-specific probes to decipher pericentric inversions has proved to be a formidable approach in characterizing chromosome rearrangements and providing further evidence on human descent.

Unique genomic sequences in human chromosome 16p are conserved in the great apes.

In humans, acute myelomonocytic leukemia (AMML) with abnormal bone marrow eosinophilia is diagnosed by the presence of a pericentric inversion in chromosome 16, involving breakpoints p13;q23 [i.e., inv(16)(p13;q23)]. A pericentric inversion involves breaks that have occurred on the p and q arms and the segment in between is rotated 180 degrees and reattaches. The recent development of a "human micro-coatasome" painting probe for 16p contains unique DNA sequences that fluorescently label only the short arm of chromosome 16, which facilitates the identification of such inversions and represents an ideal tool for analyzing the "divergence/convergence" of the equivalent human chromosome 16 (PTR 18, GGO 17 and PPY 19) in the great apes, chimpanzee, gorilla and orangutan. When the probe is used on the type of pericentric inversion characteristic of AMML, signals are observed on the proximal portions (the regions closest to the centromere) of the long and short arms of chromosome 16. The probe hybridized to only the short arm of all three ape chromosomes and signals were not observed on the long arms, suggesting that a pericentric inversion similar to that seen in AMML has not occurred in any of these great apes.

Polymorphism, monomorphism, and sequences in conserved microsatellites in primate species.

Dimeric short tandem repeats are a source of highly polymorphic markers in the mammalian genome. Genetic variation at these hypervariable loci is extensively used for linkage analysis, for the identification of individuals, and may be useful for interpopulation and interspecies studies. In this paper, we analyze the variability and the sequences of a segment including three microsatellites, first described in man, in several species of primates (chimpanzee, orangutan, gibbon, and macaque) using the heterologous primers (man primers). This region is located on the human chromosome 6p, near the tumor necrosis factor genes, in the major histocompatibility complex. The fact that these primers work in all species studied indicates that they are conserved throughout the different lineages of the two superfamilies, the Hominoidea and the Cercopithecidea, represented by the macaques. However, the intervening sequence displays intraspecific and interspecific variability. The sites of base substitutions and the insertion/deletion events are not evenly distributed within this region. The data suggest that it is necessary to have a minimal number of repeats to increase the rate of mutation sufficiently to allow the development of polymorphism. In some species, the microsatellites present single base variations which reduce the number of contiguous repeats, thus apparently slowing the rate of additional slippage events. Species with such variations or a low number of repeats are monomorphic. These microsatellite sequences are informative in the comparison of closely related species and reflect the phylogeny of the Old World monkeys, apes, and man.

Gorilla and orangutan c-myc nucleotide sequences: inference on hominoid phylogeny.

The nucleotide sequences of the gorilla and orangutan myc loci have been determined by the dideoxy nucleotide method. As previously observed in the human and chimpanzee sequences, an open reading frame (ORF) of 188 codons overlapping exon 1 could be deduced from the gorilla sequence. However, no such ORF appeared in the orangutan sequence. The two sequences were aligned with those of human and chimpanzee as hominoids and of gibbon and marmoset as outgroups of hominoids. The branching order in the evolution of primates was inferred from these data by different methods: maximum parsimony and neighbor-joining. Our results support the view that the gorilla lineage branched off before the human and chimpanzee diverged and strengthen the hypothesis that chimpanzee and gorilla are more related to human than is orangutan.

Pvu II polymorphism in the primate homologue of the mouse B144 (LST-1). A novel marker gene within the tumor necrosis factor region.

A Pvu II RFLP was mapped within the LST-1 gene, the human homologue of the mouse B144 sequence, establishing LST-1 as a new marker gene within the TNF region. We investigated the distribution of this Pvu II RFLP in 274 unrelated individuals, 132 additional HLA-DR7-positive individuals, 86 homozygous lymphoblastoid cell lines, and in four families. Seventeen of 274 individuals (6.2%) were heterozygous for the Pvu II restriction site (ADB1 = lack and ADB2 = presence of the Pvu II restriction site). In our study population the polymorphism has a much wider distribution than that previously reported in an analysis of selected haplotypes. Besides a strong association of ADB1 with HLA-B14, -DR7, we found a further association with HLA-B35. These results were also validated by family segregation studies and analyses of homozygous cell lines. In addition, five of 17 individuals carrying the ADB1 allele had HLA types other than B14 or B35, emphasizing that the presence of ADB1 is not limited to the HLA-B14, DR7 haplotype. LST-1 and its polymorphism may be used as an additional marker of the TNF region, where genes responsible for autoimmune diseases are suspected to be localized.

Evolutionary divergence of human chromosome 9 as revealed by the position of the ABL protooncogene in higher primates.

Attempts to solve the fundamental questions regarding the descent of man are dogged by superstitions and unexamined orthodoxies. The origin of humans, established a decade ago based upon cytological analysis of ape chromosomes, continues to be called into question. Although molecular methods have provided a framework for tracing the paths of human evolution, conclusive evidence remains elusive. We have used a single ABL gene probe derived from human chromosome 9 to assess the direction of change in the equivalent ape chromosomes. This approach has resulted in a few surprises which again challenge the prevailing view of early primate evolution based solely on chromosome banding patterns. The ABL proto-oncogene is present on human chromosome 9 at band q34. Similar DNA sequences presumed to represent an ABL gene, are present on chromosome 11 in chimpanzee (Pan troglodytes) but at a different relative location, indicating that the mechanism of the origin of human chromosome 9 is far more complex than has previously been suggested. Nevertheless, in gorilla (Gorilla gorilla) and orangutan (Pongo pygmaeus), the equivalent to human chromosome band 9 q34 is apparently located on chromosome 13 at a putative telomeric position and no discernible differences could be established. Despite the presence of the ABL protooncogene on human equivalent ape chromosomes, molecular systematics will continue to generate enigmas in the evolutionary context until the entire genome is sequenced.

Emergence of the keratinocyte growth factor multigene family during the great ape radiation.

The structural gene for human keratinocyte growth factor (KGF), a member of the fibroblast growth factor family, consists of three coding exons and two introns typical of other fibroblast growth factor loci. A portion of the KGF gene, located on chromosome 15, is amplified to approximately 16 copies in the human genome, and these highly related copies (which consist of exon 2, exon 3, the intron between them, and a 3' noncoding segment of the KGF transcript) are dispersed to multiple human chromosomes. The KGF-like sequences are transcriptionally active, differentially regulated in various tissues, and composed of three distinct classes of coding sequences that are 5% divergent from each other and from the authentic KGF sequence. Multiple copies of KGF-like genes were also discovered in the genomic DNAs of chimpanzee and gorilla but were not found in lesser apes (gibbon), Old World monkeys (African green monkey and macaques), mice, or chickens. The pattern of evolutionary occurrence suggests that a primordial KGF gene was amplified and chromosomally dispersed subsequent to the divergence of orangutan from African apes but before the trichotomous divergence of human, chimpanzee, and gorilla 5-8 million years ago. The appearance of a transcriptionally active and chromosomally dispersed multigene KGF family may have implications in the evolution of the great apes and humans.

A molecular phylogeny of the hominoid primates as indicated by two-dimensional protein electrophoresis.

A molecular phylogeny for the hominoid primates was constructed by using genetic distances from a survey of 383 radiolabeled fibroblast polypeptides resolved by two-dimensional electrophoresis (2DE). An internally consistent matrix of Nei genetic distances was generated on the basis of variants in electrophoretic position. The derived phylogenetic tree indicated a branching sequence, from oldest to most recent, of cercopithecoids (Macaca fascicularis), gibbon-siamang, orangutan, gorilla, and human-chimpanzee. A cladistic analysis of 240 electrophoretic characters that varied between ape species produced an identical tree. Genetic distance measures obtained by 2DE are largely consistent with those generated by other molecular procedures. In addition, the 2DE data set appears to resolve the human-chimpanzee-gorilla trichotomy in favor of a more recent association of chimpanzees and humans.