Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup

Comparative genetic analysis between human and chimpanzee may detect genetic divergences responsible for human-specific characteristics. Previous studies have identified a series of genes that potentially underwent Darwinian positive selection during human evolution. However, without a closely related species as outgroup, it is difficult to identify human-lineage-specific changes, which is critical in delineating the biological uniqueness of humans. In this study, we conducted phylogeny-based analyses of 2633 human brain-expressed genes using rhesus macaque as the outgroup. We identified 47 candidate genes showing strong evidence of positive selection in the human lineage. Genes with maximal expression in the brain showed a higher evolutionary rate in human than in chimpanzee. We observed that many immune-defense-related genes were under strong positive selection, and this trend was more prominent in chimpanzee than in human. We also demonstrated that rhesus macaque performed much better than mouse as an outgroup in identifying lineage-specific selection in humans.     

A comparative analysis of regulatory regions of the transthyretin gene in the mouse, human, and chimpanzee genomes

A full genome analysis of differences between the gene expression in the human and chimpanzee brains revealed that the gene for transthyretin, the carrier of thyroid hormones, is differently transcribed in the cerebella of these species. A 7-kbp DNA fragment of chimpanzee was sequenced to identify possible regulatory sequences responsible for the differences in expression. One hundred and thirteen substitutions were found in the chimpanzee sequence in comparison with the human sequence. About 40% of the substitutions were revealed within the repeating elements of the genome; their location and sizes did not differ from those in the corresponding fragments of the human genome, and the nucleotide sequences had a high degree of identity. A comparison of nucleotide sequences of the transthyretin region of human, chimpanzee, and mouse genes revealed substantial differences in the distribution of G + C content along the examined fragment in the human (chimpanzee) and mouse genes and allowed us to localize three sequence tracts with a higher degree of identity in the three species. One of these tracts is located in the promoter region of the gene, and the other two probably determine the specificity of transthyretin gene expression in the liver and brain. One of the conserved tracts of the chimpanzee genome was found to have a single and a triple nucleotide substitution. The triple substitution distinguishes chimpanzees from humans and mice, which have identical sequences of this site. It is likely that these substitutions are responsible for the differences in the expression levels of the transthyretin gene in the human and chimpanzee brains.     

A Comparative Analysis of Regulatory Regions of the Transthyretin Gene in the Mouse, Human, and Chimpanzee Genomes

A full genome analysis of differences between the gene expression in the human and chimpanzee brains revealed that the gene for transthyretin, the carrier of thyroid hormones, is differently transcribed in the cerebella of these species. A 7-kbp DNA fragment of chimpanzee was sequenced to identify possible regulatory sequences responsible for the differences in expression. One hundred and thirteen substitutions were found in the chimpanzee sequence in comparison with the human sequence. About 40% of the substitutions were revealed within the repeating elements of the genome; their location and sizes did not differ from those in the corresponding fragments of the human genome, and the nucleotide sequences had a high degree of identity. A comparison of nucleotide sequences of the transthyretin region of human, chimpanzee, and mouse genes revealed substantial differences in the distribution of G + C content along the examined fragment in the human (chimpanzee) and mouse genes and allowed us to localize three sequence tracts with a higher degree of identity in the three species. One of these tracts was located in the promoter region of the gene, and the other two probably determine the specificity of transthyretin gene expression in the liver and brain. One of the conserved tracts of the chimpanzee genome was found to have a single and a triple nucleotide substitution. The triple substitution distinguishes chimpanzees from humans and mice, which have identical sequences of this site. It is likely that these substitutions are responsible for the differences in the expression levels of the transthyretin gene in the human and chimpanzee brains.     

Low rates of expression profile divergence in highly expressed genes and tissue-specific genes during mammalian evolution

Evolutionary rates provide important information about the pattern and mechanism of evolution. Although the rate of gene sequence evolution has been well studied, the rate of gene expression evolution is poorly understood. In particular, it is unclear whether the gene expression level and tissue specificity influence the divergence of expression profiles between orthologous genes. Here we address this question using a microarray data set comprising the expression signals of 10,607 pairs of orthologous human and mouse genes from over 60 tissues per species. We show that the level of gene expression and the degree of tissue specificity are generally conserved between the human and mouse orthologs. The rate of gene expression profile change during evolution is negatively correlated with the level of gene expression, measured by either the average or the highest level among all tissues examined. This is analogous to the observation that the rate of gene (or protein) sequence evolution is negatively correlated with the gene expression level. The impacts of the degree of tissue specificity on the evolutionary rate of gene sequence and that of expression profile, however, are opposite. Highly tissue-specific genes tend to evolve rapidly at the gene sequence level but slowly at the expression profile level. Thus, different forces and selective constraints must underlie the evolution of gene sequence and that of gene expression.     

Comparative analysis of testis protein evolution in rodents

Genes expressed in testes are critical to male reproductive success, affecting spermatogenesis, sperm competition, and sperm-egg interaction. Comparing the evolution of testis proteins at different taxonomic levels can reveal which genes and functional classes are targets of natural and sexual selection and whether the same genes are targets among taxa. Here we examine the evolution of testis-expressed proteins at different levels of divergence among three rodents, mouse (Mus musculus), rat (Rattus norvegicus), and deer mouse (Peromyscus maniculatus), to identify rapidly evolving genes. Comparison of expressed sequence tags (ESTs) from testes suggests that proteins with testis-specific expression evolve more rapidly on average than proteins with maximal expression in other tissues. Genes with the highest rates of evolution have a variety of functional roles including signal transduction, DNA binding, and egg-sperm interaction. Most of these rapidly evolving genes have not been identified previously as targets of selection in comparisons among more divergent mammals. To determine if these genes are evolving rapidly among closely related species, we sequenced 11 of these genes in six Peromyscus species and found evidence for positive selection in five of them. Together, these results demonstrate rapid evolution of functionally diverse testis-expressed proteins in rodents, including the identification of amino acids under lineage-specific selection in Peromyscus. Evidence for positive selection among closely related species suggests that changes in these proteins may have consequences for reproductive isolation.     

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25–50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events.     

Identification of two novel primate-specific genes in DSCR.

We recently helped to complete the sequence of human chromosome 21 at a very high level of accuracy. Using this sequence we identified two novel genes, designated DSCR9 and DSCR10, in the so-called Down Syndrome Critical Region (DSCR) by computational gene prediction and subsequent cDNA cloning. Both DSCR9 and DSCR10 are expressed preferentially in testis and encode functionally unknown proteins with 149 and 87 amino acid residues, respectively. Zoo blot analysis suggested that both genes are exclusive to primate genomes such as chimpanzee, gorilla, orangutan, crab-eating monkey and African green monkey but are not present in other non-primate mammals including mouse, dog, cat, and chicken. Comparative genomic sequence analysis of DSCR9 and DSCR10 with the corresponding mouse syntenic region confirmed the lack of these genes in the mouse. These results strongly suggest that DSCR9 and DSCR10 have emerged as a new class of gene in the primate lineage during evolution.     

Detection and genetic characterization of enteroviruses circulating among wild populations of chimpanzees in Cameroon: relationship with human and simian enteroviruses

Enteroviruses (EVs), members of the family Picornaviridae, are a genetically and antigenically diverse range of viruses causing acute infections in humans and several Old World monkey (OWM) species. Despite their known wide distribution in primates, nothing is currently known about the occurrence, frequency, and genetic diversity of enteroviruses infecting apes. To investigate this, 27 chimpanzee and 27 gorilla fecal samples collected from undisturbed jungle areas with minimal human contact in Cameroon were screened for EVs. Four chimpanzee samples were positive, but none of the gorilla samples were positive. Genetic characterization of the VP1, VP4, and partial VP2 genes, the 5' untranslated region, and partial 3Dpol sequences enabled chimpanzee-derived EVs to be identified as (i) the species A type, EV76, (ii) a new species D type assigned as EV111, along with a human isolate from the Democratic Republic of Congo previously described by the International Committee on the Taxonomy of Viruses, and (iii) a new species B type (assigned as EV110) most closely related to, although a distinct type from, the SA5 isolate recovered from a vervet monkey. The identification of EVs infecting chimpanzees related to those circulating in human and OWM populations provides evidence for cross-species transmission of EVs between primates. However, the direction of transfer and the existence of primate sources of zoonotic enterovirus infections in humans require further investigation of population exposure and more extensive characterization of EVs circulating in wild ape populations.     

Comparison of chimpanzee and human leukocyte Ig-like receptor genes reveals framework and rapidly evolving genes.

The leukocyte receptor complex (LRC) on human chromosome 19 contains related Ig superfamily killer cell Ig-like receptor (KIR) and leukocyte Ig-like receptor (LIR) genes. Previously, we discovered much difference in the KIR genes between humans and chimpanzees, primate species estimated to have approximately 98.8% genomic sequence similarity. Here, the common chimpanzee LIR genes are identified, characterized, and compared with their human counterparts. From screening a chimpanzee splenocyte cDNA library, clones corresponding to nine different chimpanzee LIRs were isolated and sequenced. Analysis of genomic DNA from 48 unrelated chimpanzees showed 42 to have all nine LIR genes, and six animals to lack just one of the genes. In structural diversity and functional type, the chimpanzee LIRs cover the range of human LIRs. Although both species have the same number of inhibitory LIRs, humans have more activating receptors, a trend also seen for KIRs. Four chimpanzee LIRs are clearly orthologs of human LIRs. Five other chimpanzee LIRs have paralogous relationships with clusters of human LIRs and have undergone much recombination. Like the human genes, chimpanzee LIR genes appear to be organized into two duplicated blocks, each block containing two orthologous genes. This organization provides a conserved framework within which there are clusters of faster evolving genes. Human and chimpanzee KIR genes have an analogous arrangement. Whereas both KIR and LIR genes can exhibit greater interspecies differences than the genome average, within each species the LIR gene family is more conserved than the KIR gene family.     

Genomic divergence between human and chimpanzee estimated from large-scale alignments of genomic sequences

To study the genomic divergence between human and chimpanzee, large-scale genomic sequence alignments were performed. The genomic sequences of human and chimpanzee were first masked with the RepeatMasker and the repeats were excluded before alignments. The repeats were then reinserted into the alignments of nonrepetitive segments and entire sequences were aligned again. A total of 2.3 million base pairs (Mb) of genomic sequences, including repeats, were aligned and the average nucleotide divergence was estimated to be 1.22%. The Jukes-Cantor (JC) distances (nucleotide divergences) in nonrepetitive (1.44 Mb) and repetitive sequences (0.86 Mb) are 1.14% and 1.34%, respectively, suggesting a slightly higher average rate in repetitive sequences. Annotated coding and noncoding regions of homologous chimpanzee genes were also retrieved from GenBank and compared. The average synonymous and nonsynonymous divergences in 88 coding genes are 1.48% and 0.55%, respectively. The JC distances in intron, 5' flanking, 3' flanking, promoter, and pseudogene regions are 1.47%, 1.41%, 1.68%, 0.75%, and 1.39%, respectively. It is not clear why the genetic distances in most of these regions are somewhat higher than those in genomic sequences. One possible explanation is that some of the genes may be located in regions with higher mutation rates.     

Evolution of X-Degenerate Y Chromosome Genes in Greater Apes: Conservation of Gene Content in Human and Gorilla,But Not Chimpanzee

Compared with the X chromosome, the mammalian Y chromosome is considerably diminished in size and has lost most of its ancestral genes during evolution. Interestingly, for the X-degenerate region on the Y chromosome, human has retained all 16 genes, while chimpanzee has lost 4 of the 16 genes since the divergence of the two species. To uncover the evolutionary forces governing ape Y chromosome degeneration, we determined the complete sequences of the coding exons and splice sites for 16 gorilla Y chromosome genes of the X-degenerate region. We discovered that all studied reading frames and splice sites were intact, and thus, this genomic region experienced no gene loss in the gorilla lineage. Higher nucleotide divergence was observed in the chimpanzee than the human lineage, particularly for genes with disruptive mutations, suggesting a lack of functional constraints for these genes in chimpanzee. Surprisingly, our results indicate that the human and gorilla orthologues of the genes disrupted in chimpanzee evolve under relaxed functional constraints and might not be essential. Taking mating patterns and effective population sizes of ape species into account, we conclude that genetic hitchhiking associated with positive selection due to sperm competition might explain the rapid decline in the Y chromosome gene number in chimpanzee. As we found no evidence of positive selection acting on the X-degenerate genes, such selection likely targets other genes on the chimpanzee Y chromosome. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Hotspots of Biased Nucleotide Substitutions in Human Genes

Genes that have experienced accelerated evolutionary rates on the human lineage during recent evolution are candidates for involvement in human-specific adaptations. To determine the forces that cause increased evolutionary rates in certain genes, we analyzed alignments of 10,238 human genes to their orthologues in chimpanzee and macaque. Using a likelihood ratio test, we identified protein-coding sequences with an accelerated rate of base substitutions along the human lineage. Exons evolving at a fast rate in humans have a significant tendency to contain clusters of AT-to-GC (weak-to-strong) biased substitutions. This pattern is also observed in noncoding sequence flanking rapidly evolving exons. Accelerated exons occur in regions with elevated male recombination rates and exhibit an excess of nonsynonymous substitutions relative to the genomic average. We next analyzed genes with significantly elevated ratios of nonsynonymous to synonymous rates of base substitution (d_N/d_S) along the human lineage, and those with an excess of amino acid replacement substitutions relative to human polymorphism. These genes also show evidence of clusters of weak-to-strong biased substitutions. These findings indicate that a recombination-associated process, such as biased gene conversion (BGC), is driving fixation of GC alleles in the human genome. This process can lead to accelerated evolution in coding sequences and excess amino acid replacement substitutions, thereby generating significant results for tests of positive selection.     

Slow molecular clocks in Old World monkeys,apes, and humans

Two longstanding issues on the molecular clock hypothesis are studied in this article. First, is there a global molecular clock in mammals? Although many authors have observed unequal rates of nucleotide substitution among mammalian lineages, some authors have proposed a global clock for all eutherians, i.e., a single global rate of 2.2 x 10(-9) substitutions per nucleotide site per year. We reexamine this issue using noncoding, nonrepetitive DNA from Old World monkeys (OWMs), chimpanzee, and human. First, using the minimal date of 6 MYA for the human-chimpanzee divergence and more than 2.5 million base pairs of genomic sequences from human and chimpanzee, we estimate a maximal rate of 0.99 x 10(-9) for noncoding, nonrepetitive genomic regions for these two species. This estimate is less than half of the proposed global rate and much smaller than the commonly used rate (3.5 x 10(-9)) for eutherians. Further, using a minimal date of 23 MYA for the human-OWM divergence, we estimate a maximal rate of 1.5 x 10(-9) for both introns and fourfold degenerate sites in humans and OWMs. In addition, with the New World monkey (NWM) lineage as an outgroup, we estimate that the rate of substitution in introns is 30% higher in the OWM lineage than in the human lineage. Clearly, there is no global molecular clock in eutherians. Second, although many studies have indicated considerable variation in the mutation rate among regions of the mammalian genome, a recent study proposed a uniform rate. Using new and existing intron sequence data from higher primates, we find significant rate variation among genomic regions and a positive correlation between the rate of substitution and the GC content, refuting the claim of a uniform rate.     

Genomic organization, chromosome location, and expression analysis of mouse beta-synuclein, a candidate for involvement in neurodegeneration

The synuclein family of proteins is a group of primarily brain-expressed polypeptides that show a high degree of amino acid conservation. alpha-Synuclein is the best known of the synuclein family, as it is a major component of the Lewy body, a cytoplasmic inclusion characteristic of Parkinson's disease as well as a variety of related neurodegenerative disorders. With the discovery that mutations in alpha-synuclein can cause Parkinson's disease, a potential role for the other synuclein family members in neurodegenerative disease is being considered. beta-Synuclein in particular may deserve special attention, as it is co-expressed with alpha-synuclein at presynaptic nerve terminals, is subject to phosphorylation by Ca(2+) calmodulin protein kinase II, appears important for neural plasticity, and forms aggregates in the brains of patients with Parkinson's disease and a related disorder. To facilitate study of beta-synuclein, we have cloned the mouse beta-synuclein gene (Sncb) and determined its genomic organization, size, and intron-exon structure. Using an interspecific backcross mapping panel from The Jackson Laboratory, we were then able to localize Sncb to chromosome 13 at the MGD 35.0 cM position. Like the human beta-synuclein gene, Sncb appears to consist of six exons separated by five introns. Unlike the human beta-synuclein gene, the mouse ortholog possesses a variant GC 5' splice donor sequence at the exon 4 - intron 4 boundary in a highly conserved splice junction consensus. Northern blot analysis and Western blot analysis both indicate that Sncb is highly expressed in the brain. Knowledge of the genomic organization and expression pattern of Sncb will allow functional studies of its potential role in neurodegeneration to commence in the mouse.     

Conservative evolution in duplicated genes of the primate Class I ADH cluster

Humans have seven alcohol dehydrogenase genes (ADH) falling into five classes. Three out of the seven genes (ADH1A, ADH1B and ADH1C) belonging to Class I are expressed primarily in liver and code the main enzymes catalyzing ethanol oxidization. The three genes are tandemly arrayed within the ADH cluster on chromosome 4 and have very high nucleotide similarity to each other (exons: >90%; introns: >70%), suggesting the genes have been generated by duplication event(s). One explanation for maintaining similarity of such clustered genes is homogenization via gene conversion(s). Alternatively, recency of the duplications or some other functional constraints might explain the high similarities among the genes. To test for gene conversion, we sequenced introns 2, 3, and 8 of all three Class I genes (total>15.0 kb) for five non-human primates--four great apes and one Old World Monkey (OWM)--and compared them with those of humans. The phylogenetic analysis shows each intron sequence clusters strongly within each gene, giving no evidence for gene conversion(s). Several lines of evidence indicate that the first split was between ADH1C and the gene that gave rise to ADH1A and ADH1B. We also analyzed cDNA sequences of the three genes that have been previously reported in mouse and Catarrhines (OWMs, chimpanzee, and humans) and found that the synonymous and non-synonymous substitution (dN/dS) ratios in all pairs are less than 1 representing purifying selection. This suggests that purifying selection is more important than gene conversion(s) in maintaining the overall sequence similarity among the Class I genes. We speculate that the highly conserved sequences on the three duplicated genes in primates have been achieved essentially by maintaining stability of the hetero-dimer formation that might have been related to dietary adaptation in primate evolution.     

Functional impact of transposable elements using bioinformatic analysis and a comparative genomic approach

A dual coding event, which is the translation of different isoforms from a single gene, is one of the special patterns among the alternative splicing events. This is an important mechanism for the regulation of protein diversity in human and mouse genomes. Although the regulation for dual coding events has been characterized in a few genes, the individual mechanism remains unclear. Numerous studies have described the exonization of transposable elements, which is the splicing mediated insertion of transposable element sequence fragments into mature mRNAs. Therefore, in this study, we investigated the number of transposable element (TE)-derived dual coding genes in human, chimpanzee and mouse genomes. TE fusion exons appeared in the dual coding regions of 309 human genes. Functional protein domain alterations by TE-derived dual coding events were observed in 129 human genes. Comparative TE-derived dual coding events were also analyzed in chimpanzee and mouse orthologs. Seventy chimpanzee orthologs had TE-derived dual coding events, but mouse orthologs did not have any TE-derived dual coding events. Taken together, our analyses listed the number of TE-derived dual coding genes which could be investigated by experimental analysis and suggested that TE-derived dual coding events were major sources for the functional diversity of human genes, but not mouse genes.     

Widespread Genomic Signatures of Natural Selection in Hominid Evolution

Selection acting on genomic functional elements can be detected by its indirect effects on population diversity at linked neutral sites. To illuminate the selective forces that shaped hominid evolution, we analyzed the genomic distributions of human polymorphisms and sequence differences among five primate species relative to the locations of conserved sequence features. Neutral sequence diversity in human and ancestral hominid populations is substantially reduced near such features, resulting in a surprisingly large genome average diversity reduction due to selection of 19–26% on the autosomes and 12–40% on the X chromosome. The overall trends are broadly consistent with “background selection” or hitchhiking in ancestral populations acting to remove deleterious variants. Average selection is much stronger on exonic (both protein-coding and untranslated) conserved features than non-exonic features. Long term selection, rather than complex speciation scenarios, explains the large intragenomic variation in human/chimpanzee divergence. Our analyses reveal a dominant role for selection in shaping genomic diversity and divergence patterns, clarify hominid evolution, and provide a baseline for investigating specific selective events.     

Growth hormone locus expands and diverges after the separation of New and Old World Monkeys

While most mammals including the prosimians have a single copy of the growth hormone (GH) gene, anthropoids possess a cluster of GH-related genes. Throughout the evolution of the main anthropoid groups [New World Monkeys (NWM), Old World Monkeys (OWM), and apes], two features stand out of the GH loci. The first is the appearance of chorionic somatommamotropin hormone (CSH) genes within the OWM lineage and the second is the expansion of the loci intergenic regions in the OWM and apes. In relation with this loci expansion, the NWM possess intergenic regions of homogeneous lengths (3.5 kb). In contrast, heterogeneous lengths (6 and 13 kb) have been reported for species of the OWM. At the present, none of the OWM genomic GH loci organizations have been described. Here, we report the genomic organization of the GH locus in the rhesus monkey, this locus has six GH-related genes separated by five intergenic regions. The 5' end gene (GH-1) encodes for the pituitary GH and is followed by CSH-1, GH-2, CSH-2, CSH-3 and CSH-4 genes. The five intergenic regions have heterogeneous lengths and also present more or less the same Alu distribution as the human GH locus. To analyze the events that contributed to the extension of the intergenic regions of the GH locus and the emergence of the regulatory elements, the five GH locus intergenic regions of the spider monkey (NWM) were sequenced. The results of comparing the loci from both species suggest that the long intergenic regions (13 kb) of the rhesus GH locus share a common ancestor with the 3.5 kb intergenic regions of the spider monkey. However, the observed increased length of the former is due to an insertion (approximately 8.7 kb) at their 3' end. Interestingly in this insert, we discovered a DNA element resembling the enhancer of the CSH genes of the human GH locus. On the other hand, we observed that the short intergenic regions (6 kb) increased by a different recombination event.