Skip to main content
Advertisement
  • Loading metrics

Heterogeneous Genomic Molecular Clocks in Primates

  • Seong-Ho Kim ,

    Contributed equally to this work with: Seong-Ho Kim, Navin Elango

    Affiliation School of Biology, Georgia Institute of Technology, Atlanta, Georgia, United States of America

  • Navin Elango ,

    Contributed equally to this work with: Seong-Ho Kim, Navin Elango

    Affiliation School of Biology, Georgia Institute of Technology, Atlanta, Georgia, United States of America

  • Charles Warden,

    Affiliation School of Biology, Georgia Institute of Technology, Atlanta, Georgia, United States of America

  • Eric Vigoda,

    Affiliation College of Computing, Georgia Institute of Technology, Atlanta, Georgia, United States of America

  • Soojin V Yi

    To whom correspondence should be addressed. E-mail: soojinyi@gatech.edu

    Affiliation School of Biology, Georgia Institute of Technology, Atlanta, Georgia, United States of America

Abstract

Using data from primates, we show that molecular clocks in sites that have been part of a CpG dinucleotide in recent past (CpG sites) and non-CpG sites are of markedly different nature, reflecting differences in their molecular origins. Notably, single nucleotide substitutions at non-CpG sites show clear generation-time dependency, indicating that most of these substitutions occur by errors during DNA replication. On the other hand, substitutions at CpG sites occur relatively constantly over time, as expected from their primary origin due to methylation. Therefore, molecular clocks are heterogeneous even within a genome. Furthermore, we propose that varying frequencies of CpG dinucleotides in different genomic regions may have contributed significantly to conflicting earlier results on rate constancy of mammalian molecular clock. Our conclusion that different regions of genomes follow different molecular clocks should be considered when inferring divergence times using molecular data and in phylogenetic analysis.

Synopsis

The rate at which mutations accumulate in a genome, referred as a “molecular clock,” is an instrumental tool in molecular evolution and phylogenetics. Different types of mutations occur via distinctive molecular pathways. In particular, while most mutations occur from errors in DNA replication, spontaneous deamination of methylated CpG dinucleotides is another important source of mutation in mammalian genomes. Molecular clock studies typically combined all types of mutations together. In this paper, the authors analyze molecular clocks of replication-origin and methylation-origin mutations separately. By utilizing high-quality sequence data from several primate species and fossil calibration, the authors demonstrate that the two types of mutations follow statistically different molecular clocks. Methylation-origin mutations accumulate relatively constantly over time, while replication-origin mutations scale with generation-times. Therefore, the genomic molecular clock, as a whole, is shaped by the molecular origins of mutations that have accumulated over time. The authors' results have direct implications on phylogenetic analyses, estimation of species divergence dates, and studies of the mechanisms and processes of evolution, where molecular clocks are imperative.

Introduction

Organisms with longer generation-time tend to exhibit slower molecular clock than those with shorter generation-time, an effect known as “generation-time effect” [15]. However, the extent (or even the existence) of generation-time effect is of significant debate [3,6,7]. An opposing theory posits that molecular evolution occurs relatively constantly over time: in other words, molecular clocks are time dependent [6,8]. Here we show that molecular evolution follows both generation-time–dependent and time-dependent molecular clocks, depending upon the molecular origins of the mutations considered.

A generation-time–dependent molecular clock implies that the majority of single nucleotide substitutions in germlines arise from errors during DNA replication [3,9]. However, some mutations may occur independently from DNA replication. This is especially pertinent for transitions at CpG dinucleotides (henceforth, CpG substitutions). CpG substitutions are the most frequent single nucleotide substitutions in vertebrate genomes, accounting for more than a quarter of all substitutions between the genomes of human and chimpanzee [10,11]. Naturally, they play critical roles in several key genetic mechanisms and disease [1216].

CpG dinucleotides are hypermutable because the cytosines in CpG dinucleotides are targets of DNA methylation in vertebrate genomes [17]. Methylated cytosine rapidly mutates to thymine via spontaneous deamination, causing a C to T (G to A in the complementary strand) transition [17,18]. While DNA replication occurs in a specialized stage of the cell cycle, methylation is not confined to replicating DNA: germline cells are methylated early in their development and stay methylated until global demethylation occurs after fertilization [19,20]. Therefore, methylation-origin mutations will accumulate in a rate proportional to the total amount of time germ cells are methylated between generations. In other words, the molecular clock at CpG dinucleotides should be relatively constant over time.

Indeed, statistical inferences using approximately 2 Mbp of sequence data have suggested that CpG substitutions follow relatively constant molecular clock in mammals [21]. In addition, a recent analysis of male mutation bias in humans and chimpanzees have shown that CpG dinucleotides exhibit much lower male mutation bias than other sites [22]. Since male-mutation bias is caused by the more frequent DNA replications in male germlines compared to female germlines [14], the finding that there is lower male mutation bias in CpG dinucleotides is consistent with the idea that CpG substitutions follow a relatively time-dependent molecular clock.

In this paper, we sought to directly compare genomic molecular clocks of CpG dinucleotides and other sites. To achieve this goal, we focused on catarrhines, specifically two hominoid species (human and chimpanzee) and two Old World monkeys (rhesus macaque and baboon). These four species are chosen because they satisfy two criteria. First, because these species are closely related, we can identify sites that have been part of a CpG dinucleotide in recent past (CpG sites) and other sites with high confidence [23]. Second, hominoids and Old World monkeys have markedly differently generation times. According to Gage [24], average generation times in Old World monkeys is 11.4 years, while in chimpanzees and humans, they are 22 and 28 years, respectively. As a consequence of the difference in generation times, evolutionary rates of replication-dependent substitutions are slower in hominoids than in Old World monkeys [2,4,25].

Utilizing genomic data from these species, we demonstrate that indeed CpG substitutions exhibit a relatively time-dependent molecular clock, in contrast to generation-time–dependent genomic molecular clock. Furthermore, we propose that heterogeneous molecular clocks among different genomic regions may have contributed to conflicting earlier results on the degree of generation-time effect in mammals.

Results/Discussion

Slower Molecular Evolution of Hominoid Genomes than Old World Monkey Genomes

We first reevaluated the difference in evolutionary rates between hominoids and Old World monkeys. We analyzed approximately 28 Mbp of genomic sequence alignments to compare rates in human (a hominoid) and baboon (an Old World monkey) using a relative rate test [4,26]. Sequence data from marmoset (a New World monkey) were used as an outgroup. We found that rates in humans are on average 28.4% slower than those in baboons in introns and intergenic regions (Table 1, p < 0.001), confirming earlier results [2,4,27]. Because data used in this analysis account for approximately 1% of the human genome and from several different chromosomes, we can conclude that the canonical genomic molecular clocks in primates exhibit significant generation-time effect.

thumbnail
Table 1.

Hominoid-Rate Slowdown Tested Using Genomic Sequence Data from Human, Baboon, and Marmoset

https://doi.org/10.1371/journal.pgen.0020163.t001

We also constructed a five-species phylogeny of human, chimpanzee, baboon, macaque, and marmoset using data for 1.9 Mbp of sequences orthologous to the human chromosome 7 (hg17.chr7: 115404472–117281897; ENCODE region ENm001). High-quality sequence data are available for all five species analyzed in this study. Figure 1 shows a Neighbor-Joining tree [28] of the five species. Focusing on the ancestral hominoid and ancestral Old World monkey branches, the ratio of the number of substitutions in the Old World monkey branch to the hominoid branch is approximately 1.36, similar to the values estimated from the comparison between the human and baboon genomes. These results confirm the “hominoid rate slowdown” theory proposed more than 40 yeasr ago [9,25].

thumbnail
Figure 1. A Neighbor-Joining Tree of Five Primate Species, Generated Using High-Quality Data from the Encode Region ENm001

The numbers of substitutions per 100 sites in each branch, using the two-parameter correction [58], are shown.

https://doi.org/10.1371/journal.pgen.0020163.g001

Our next goal was to compare the molecular clocks at CpG and non-CpG sites separately. However, because of the difficulty in correcting for multiple hits, we cannot easily analyze substitutions at CpG sites in this phylogenetic setting. Therefore, we proceeded to use data only in catarrhines, where we can accurately infer rates in CpG and non-CpG sites [12,22,23].

Different Molecular Clocks of CpG Sites and Non-CpG Sites

We constructed four-species alignments of two hominoids (human and chimpanzee) and two Old World monkeys (rhesus macaque and baboon) (Figure 2). These species pairs provide a unique opportunity to study time-dependent and generation-time–dependent clocks. Critical to our work, the divergence time between the hominoid pair is similar to that of the Old World monkey pair [27,29,30]. The split between human and chimpanzee is estimated to be 6 to 8 million years ago (Mya), based upon fossil records. In particular, the earliest fossil hominin, Sahelanthropus tchadensis, has been dated to late Miocene, at least 7 Mya [30,31]. The split between rhesus macaque and baboon is calibrated by using an estimate for the split between macaques and papionins. The earliest fossil evidence of papionins is dated to be 6 to 8 Mya [27,29]. Therefore, divergence times of the two species within each pair are similar. In other words, TO/TH ≈ 1 (Figure 2). In contrast to this similarity of within-pair divergence times, evolutionary rates are known to differ between these two groups: as explained in the introduction and demonstrated above, genomic evolutionary rates in hominoids are slower than rates in Old World monkeys.

thumbnail
Figure 2. Phylogeny of the Four Taxa Analyzed in This Study

TO denotes the time since the split between the two Old World monkey species, and TH denotes the time since the split between the two hominoids. Fossil records suggest that TO and TH are very close to each other. X and Y denote the common ancestors of human-chimpanzee and of macaque-baboon, respectively. The genetic divergence between the two hominoid species (KH) is the sum of KHX and KCX. Likewise, KO is the sum of KMY and KBY.

https://doi.org/10.1371/journal.pgen.0020163.g002

We have two contrasting predictions for a time-dependent versus a generation-time–dependent molecular clock. For replication-origin (hence, generation-time–dependent) mutations, the pairwise sequence divergence in the Old World monkey pair (KO = KMY + KBY in Figure 2) should be greater than the pairwise sequence divergence in the hominoid pair (KH = KHX + KCX in Figure 2). On the other hand, a time-dependent molecular clock predicts that KO is similar to KH.

We examined the molecular clocks in CpG and non-CpG sites separately (see Materials and Methods). To directly compare mutations caused by deamination of methylated cytosines to other transitions occurring during replication, we first analyzed only C-to-T (and G-to-A) transitions. A distinctive pattern emerged: KO/KH is 1.03 in CpG sites (95% confidence interval [CI], 0.92 to 1.15), while it is 1.31 in non-CpG sites (95% CI, 1.25 to 1.37). These two types of sites clearly harbor different molecular clocks. Similar trends were discovered when introns and intergenic regions are considered separately, or when repetitive and nonrepetitive sequences are compared separately (Figure 3).

thumbnail
Figure 3. Contrasting Molecular Clocks of Transitions at CpG Sites versus Those at Non-CpG Sites

The y-axis shows the rate difference in the baboon-macaque pair to that in the human-chimpanzee pair. The Old World monkey pair has accumulated significantly more transitions in non-CpG sites, as expected by the generation time effect. In contrast, transitions at CpG sites, which are primarily of methylation origin, show no difference between the two pairs. Data are shown for all sites, repetitive sites (as identified from the RepeatMasker program [57]), and nonrepetitive sites (after removing repetitive sites). Confidence intervals are generated by bootstrapping 10,000 times.

https://doi.org/10.1371/journal.pgen.0020163.g003

We then considered all single nucleotide substitutions that occurred in CpG and non-CpG sites and found the same pattern. The ratio KO/KH in non-CpG sites is 1.18 (95% CI, 1.15 to 1.22). In comparison, in CpG sites, KO/KH is 1.00 (95% CI, 0.89 to 1.11). Again, the results are similar when introns and intergenic regions are considered separately, or when repetitive and nonrepetitive sequences are compared separately.

Because human-chimpanzee (hominoid pair) and rhesus macaque-baboon (Old World monkey pair) are extremely closely related, estimates of pairwise sequence divergence are affected by common ancestral polymorphism [3234]. The common ancestor of the human and chimpanzee is thought to have much larger effective population size than the current human population [35,36]. Rhesus macaque and baboon also harbor comparable levels of genetic diversity to hominoids. For example, Rogers and Kidd [37] reported the nucleotide diversity of Papio hamadryas to be approximately 0.3%. Wall et al. [38] estimated a nucleotide diversity of 0.13% in a noncoding region of rhesus macaques.

Such substantial ancestral polymorphism will effectively reduce the observed rate difference between hominoid and Old World monkey pair: the observed pairwise divergence between rhesus macaque and baboon (KO) is the sum of ancestral diversity (πY, see Figure 2) and the fixed difference between rhesus macaque and baboon (denoted as PO). Likewise, the pairwise divergence between human and chimpanzee, KH = πX + PH. We are interested in the ratio PO/PH while we only have access to KO/KH. When comparing distantly related species, the level of ancestral diversity is negligible relative to the fixed difference. However, between closely related species such as human-chimpanzee and macaque-baboon, ancestral diversity is substantial compared to the fixed difference. For example, πX can be as much as ½ PH [35]. Therefore, KO/KH will underestimate PO/PH.

To address this concern, we used the estimates obtained for CpG and non-CpG sites in hominoids [22] to correct for the effect of ancestral polymorphism. After such corrections, KO/KH for non-CpG sites is 1.18 to 1.26 (Table 2). In contrast, in CpG sites, KO/KH is close to 1.00 even after correcting for the effect of ancestral polymorphism using estimates for CpG sites (Table 2). However, these values should be taken with caution, given the uncertainties associated with ancestral diversity as well as with divergence time estimated from fossil records.

thumbnail
Table 2.

The Ratio of the Pairwise Divergence between Macaque and Baboon (KO) to the Pairwise Divergence between Human and Chimpanzee (KH), Using All Substitutions

https://doi.org/10.1371/journal.pgen.0020163.t002

For completeness, we also analyzed the rate difference for CpG and non-CpG sites using the above three-species alignment (human, baboon, and marmoset). Even though this comparison is less reliable due to the difficulty in correcting for multiple hits (see above), we obtained similar results. We observe that the non-CpG sites (the majority of sites) show substantial rate difference between the human and the baboon genomes. In contrast, CpG sites show little difference in evolutionary rates between hominoid and Old World monkeys (Table 1).

In summary, CpG and non-CpG sites show statistically different molecular clocks in various phylogenetic comparisons, indicating that the difference in two types of molecular clocks is a salient picture of molecular evolution in primate genomes.

Factors that May Affect KO/KH for CpG and Non-CpG Sites

Here we review some of the potential factors that can affect our conclusions. An important assumption in our work is that the divergence time between the hominoid pair is similar to that of the Old World monkey pair. This was mainly based upon fossil records [27,29,30]. However, because fossil records are inherently associated with large variance in dates, let us consider the inference from molecular data.

If we measure the divergence between the Old World monkey pair to that between the hominoid pair in the five species phylogeny shown in Figure 1 (equivalent to KO/KH in Figure 2), it is 1.2. This is different from the ratio obtained from the comparison of the ancestral Old World monkey branch to the hominoid branch, which was 1.36. The discrepancy between these two estimates can be explained by at least two mechanisms, which are not mutually exclusive of each other.

First, as mentioned earlier, estimating evolutionary rates between closely related species, such as human-chimpanzee and macaque-baboon, is significantly affected by ancestral polymorphism [3234]. If we use estimates of the ancestral polymorphism in hominoids [35,36] to correct for the effect of ancestral polymorphism, the ratio of KO/KH increases, close to the value estimated from the ancestral branch. For example, if we assume that the average nucleotide diversities of the ancestral Old World monkey and hominoid populations were 0.4%, the corrected ratio of KO/KH increases to 1.32.

The second possibility is that the actual time in the Old World monkey pair (TO) is slightly shorter than the time in the hominoid pair (TH). Because fossil records provide only the “minimum” divergence time between lineages, the actual divergence time can differ significantly, and the divergence of human and chimpanzee may have occurred before the divergence of macaque and baboon. Therefore, KO/KH will underestimate the true rate difference. According to this possibility, the CpG clock in our data also underestimates the actual rate difference, indicating that some fraction of CpG substitutions follows a generation-time–dependent molecular clock. We believe that this scenario at least partially explains the observed discrepancy, because some substitutions at CpG sites occur during replication. This interpretation is also in accord with the weak but still significant male mutation bias in hominoids [22].

Our study uncovered significant heterogeneity in the degree of generation time effect among different types of single nucleotide substitutions. In particular, when substitutions are divided into transitions and transversions, the latter exhibited less generation-time effect than transitions. In fact, in CpG sites, there were more transversions in the human-chimpanzee pair than in the baboon-macaque pair (58 versus 39). However, the numbers are rather small (since most substitutions at CpG sites are transitions due to methylation), so it is not clear whether this reflects a true underlying pattern. In non-CpG sites, the ratio KO /KH estimated from transitions was 1.31, while the ratio from transversions was 1.14 (the overall ratio was 1.18). Whether this discrepancy reflects differences in molecular mechanisms between transitions and transversions is an interesting question and should be pursued further.

Effect of CpG Dinucleotides on Hominoid Rate Slowdown and Mammalian Molecular Clock

Our findings shed important light on the controversy over mammalian molecular clock. Generation-time effect was clearly demonstrated when closely related species were compared or when noncoding sequences were used [21,27]. However, among relatively distant mammalian species, weak generation-time effect was observed [6,26]. Note that due to sequence availability and alignability, synonymous sites were often used when comparing distantly related species.

We propose that varying proportions of CpG dinucleotides in different data sources can contribute to conflicting conclusions on the nature of genomic molecular clocks. Three observations led to this hypothesis. First, CpG molecular clock runs much faster than clocks at other sites, at least in primates. Assuming that human and chimpanzee diverged 7 Mya [30], we estimate that CpG sites and non-CpG sites undergo single nucleotide substitutions at a rate of 1.03 × 10−8 per site per year and 0.68 × 10−9 per site per year, respectively, from our data. Second, molecular clocks at CpG sites are relatively constant over time. Third, the proportion of CpG dinucleotides is heterogeneous among different genomic regions [39]. In particular, 4-fold degenerate sites are enriched with CpG sites, over 10% [39], while noncoding regions have less than 3% CpG dinucleotides [22,39]. Hence, molecular clocks in regions with relatively abundant CpG sites (such as 4-fold degenerate sites) may be dominated by the rapid and time-dependent CpG clock, while regions relatively devoid of CpG sites (such as noncoding regions) follow generation-time–dependent molecular clock.

To investigate this prediction, we compared results from different studies in Table 3, focusing on two comparisons: between hominoids and Old World monkeys (hominoid rate slowdown), and between primates and rodents. Note that earlier studies on molecular clock did not consider CpG content as a determinant of molecular clock. Therefore, they did not investigate the effect of CpG content on molecular clock. Because some studies used noncoding regions while others used 4-fold degenerate sites, different studies analyzed different data in relation to CpG content (Table 3). We did not include the results from [6] in this table, because they removed a substantial amount of data that did not pass the “homogeneity test,” and the relationship between this test and CpG dinucleotide content is not clear. For example, they discarded 46% of the data in their human-mouse comparison [6].

thumbnail
Table 3.

Rate Differences between Lineages from Various Data Sources

https://doi.org/10.1371/journal.pgen.0020163.t003

We can now compare how the data in Table 3 fit our hypothesis. First, when we compare results from all sites, the rate difference between lineages is greater in noncoding regions than in 4-fold degenerate sites. Moreover, in noncoding regions, the rate difference for CpG sites is lower than for all sites or non-CpG sites. Similarly, in 4-fold degenerate sites, the rate difference in non-CpG sites is higher than in all sites. These trends support our hypothesis.

Since we have reasonable estimates of CpG and non-CpG rates in primates (see above), we can investigate how well our hypothesis fits the data in detail. The number of substitutions in hominoids since the split from Old World monkeys can be approximated as where p is the proportion of CpG sites, kCpG and knon-CpG represent substitution rates per site per year in CpG sites and non-CpG sites, respectively, and T is the time since the split. The observed ratio of Old World monkey branch to hominoid branch can then be expressed as where r represents the ratio of the branch lengths determined by the generation-time–dependent molecular clock. Figure 4 shows this ratio as a function of p, using the rates inferred from our data. In case when r = 1.4, the observed ratios from regions with 12% and 2.5% CpG dinucleotides (analogous to 4-fold degenerate sites and intergenic regions) are 1.12 and 1.29, respectively.

thumbnail
Figure 4. The Proportion of CpG Sites in Data Affects the Degree of Hominoid-Rate Slowdown

We considered a simple model in which all sites can be classified into either CpG sites or non-CpG sites and estimated evolutionary rates in hominoids from the human-chimpanzee comparison. The x-axis is the proportion of CpG sites in the data. The y-axis is the observed degree of hominoid rate slowdown, shown as the ratio of the substitution rate in Old World monkeys to the rate in hominoids, given the “true” ratio (determined by the generation-time effect), depicted as r. While regions relatively devoid of CpG sites will reflect the true generation-time effect, the observed ratio approaches 1 as the data include more CpG sites (i.e., the substitution rate in hominoids and Old World monkeys will be similar). Data points for when data consists of 2.5% and 12% CpG sites for r = 1.3 and 1.4 are shown for convenience.

https://doi.org/10.1371/journal.pgen.0020163.g004

We compared these theoretical expectations to observed values by analyzing rates between hominoids and Old World monkeys in 4-fold degenerate sites, from 41 autosomal genes (Table S2). The proportion of 4-fold degenerate sites that belong to CpG dinucleotides in any of the three species compared in this dataset is 11.0%. This is likely an underestimate of the true proportions of sites that have been part of a CpG dinucleotide, since the divergence time between the three species is rather long. The ratio of the Old World monkey branch to the hominoid branch was 1.09 when all sites were used (Table 3). When we removed CpG-prone sites (sites preceded by C or followed by G, as used in [12,23,40]) from the 4-fold degenerate sites, the aforementioned ratio was increased to 1.27 (Table 3). Recall, when only noncoding sites were used, this ratio was 1.28 (Table 1), which increased to 1.31 when we removed CpG sites. The proportion of sites that belong to CpG dinucleotides in noncoding sites in our data is 2.5%. Therefore, these values are in excellent accord with the above-mentioned model.

It should be noted, however, that the above model ignores other factors that affect regional mutation rate variation, such as GC content and recombination [4,41]. Also, as discussed above, different mutations (such as transitions and transversions) may have different substitution rates between lineages. Hence, partitioning rates into only two categories is likely to be a simplification. Furthermore, identifying sites that have been part of a CpG dinucleotide in the past is a challenging problem [42,43]. Lineage-specific rates are also affected by ancestral generation times and effective population sizes. Further studies are necessary to determine the roles of generation-time–dependent and time-dependent molecular clocks on genome evolution.

Nevertheless, it is clear that the heterogeneity of molecular clocks due to different mutational origins can significantly alter rate differences between taxa. This effect should be taken into account when molecular clocks are used to infer divergence times and to reconstruct phylogenetic history.

Materials and Methods

Noncoding data mining and assembly.

Because accurate identification of CpG sites is critical in our analyses, we used two precautions. First, we analyzed sequences between closely related primates only. Earlier studies have shown that within catarrhines (hominoids and Old World monkeys), we can directly derive rates of CpG substitutions using comparative methods. Specifically, we can confidently determine “CpG sites” (sites for which the ancestral state was part of a CpG) and extract rates of CpG substitutions using parsimony [12,22,23]. Moreover, we can also identify sites that have not been a part of CpG dinucleotides (non-CpG sites), to be used as a control for replication-origin substitutions [12,22,23]. Second, we only used high-quality sequence data, because data obtained from whole genome assemblies include errors in sequencing and assembly that can cause erroneous conclusions regarding rate difference between lineages [34,44].

For the human-baboon-marmoset dataset, we obtained approximately 28 Mbp of high-quality data (BAC-based) from the ENCODE project [45].

For the human-chimpanzee-baboon-macaque (HCBM) dataset, we mined high-quality BAC-based sequences from GenBank. The HCBM dataset consists of BAC-based sequence data orthologous to human Chromosome 7 (hg17.chr7:114505472–117281897; Encode region ENm001). This is obtained by aligning NT_086357.2 [46], NT_165329.1 (chimpanzee), NT_086378.3 (baboon), and NT_165339.1 (macaque) sequences.

We assembled additional orthologous alignments among the four species using the following procedure. First, we searched the GenBank database for sequences from baboon (Papio anubis or P. hamadryas), macaque (Macaca mulatta), and chimpanzee (Pan troglodytes) BAC clones. We obtained sequence data for 377, 276, and 1,641 BACs from baboon, macaque, and chimpanzee, respectively. Next, we identified orthologous BAC clones among these species, using BLAST [47] and other methods as in [48]. We found 25 baboon BAC clones that had both macaque and chimpanzee orthologs. We then localized orthologous human region for each of these 35 orthologous clones using BLAT [49]. We reconfirmed the orthology between baboon, chimpanzee, and macaque BAC clones by ensuring that the regions where these BAC clones independently map to the human genome overlap with each other. We then removed the BAC clones overlapping with ENm001. Finally, we removed sequences from sex chromosomes. As a result, we obtained 16 genomic regions, shown in Table S1.

Analysis of 4-fold degenerate sites.

For primate comparison, all sequence data for the primate 4-fold degenerate site comparisons were downloaded from GenBank [50]. Accession numbers for all genes used in primate comparison are available in Table S2. A portion of the homologous genes in this dataset was also identified via the HOVERGEN database [51]. Sequence data for the human-mouse-dog comparison were downloaded from the Ensembl database [52]. Any genes that underwent recent gene duplications or did not meet the stringent minimum length of 445 nucleotides were removed from the dataset. Sequences were aligned using CLUSTALW [53] via a BioPerl package [54]. After alignment of homologous genes, any genes containing lineages with a negative K4 value were removed from the dataset.

For primate-rodent comparison, known genes from human, mouse, and dog were downloaded from Ensembl [52]. To find orthologous sequences, we used the OrthoMCL algorithm [55], which uses an all-to-all BLASTP results to generate a graph of orthologs and paralogs. We used default parameters except for E-value < 10−10 to ensure orthology. As a result, we constructed 3,494 orthologous gene trios among the three species. The next steps were performed as described in the primate comparison described above.

Sequence curation, data annotation, and statistical analyses.

CpG islands were identified using the algorithm by Takai and Jones [56] with the following conditions: GC content greater than 55%, observed/expected CpG contents greater than 0.65, length 200 or greater. Since the majority of CpG islands are hypomethylated and do not reflect substitutions of methylation origin, we removed them from further analysis.

Repetitive elements were annotated using the RepeatMasker program [57]. Noncoding regions are identified as in Elango et al [34].

The two-parameter model [58] was used to correct for multiple hits. We used a relative rate test [4,26] to test for rate difference between hominoid and Old World monkeys using New World monkey species as outgroup (Table S2). To compare rate difference between human and mouse, we used dog as an outgroup.

For classification and rate estimation of CpG sites and non-CpG sites, we used the method in Meunier et al. [12] to identify CpG and non-CpG sites. Specifically, CpG sites are defined as the middle base of the following patterns: XNG/XCG/XCG/XCG, with X denoting any nucleotide except C to avoid overlapping CpGs. N can occur in any of the four sequences. Sites fitting the complementary pattern (CGY/CGY/CGY/CNY, Y not G) are also considered as CpG sites. As a control, sites expected to have never been part of a CpG dinucleotides since the last common ancestor of the four species (“non-CpG sites”) are defined as sites not preceded by C nor followed by G [12,22]. Sites that do not satisfy either classification are defined as “ambiguous sites” and excluded from the analysis. A simulation study has shown that this classifying scheme can accurately identify CpG sites and non-CpG sites in catarrhines [23]. Substitutions are then inferred using unweighted parsimony using only such sites. Confidence intervals for estimated rates are derived from bootstrapping 10,000 times.

Supporting Information

Table S1. Accession Numbers of Orthologous Baboon, Chimpanzee, and Macaque BACs and Their Locations on the Human Genome (hg 17; NCBI build 35)

* This region has a partial overlap with ENm001. The portion that overlapped with ENm001 was removed.

https://doi.org/10.1371/journal.pgen.0020163.st001

(58 KB DOC)

Table S2. Accession Numbers for Genes Used in Primate Fourfold Degenerate Site Comparison

https://doi.org/10.1371/journal.pgen.0020163.st002

(103 KB DOC)

Accession Numbers

The National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov) accession numbers for human, chimpanzee, baboon, macaque, and marmoset are NT_086357.2, NT_165329.1, NT_086378.3, NT_165339.1, and NT_086504.2, respectively.

Acknowledgments

We thank Ryan Raaum for information on primate generation times, Adam Eyre-Walker for discussions, and several anonymous reviewers for comments on the manuscript.

Author Contributions

EV and SVY conceived and designed the experiments. SHK, NE, CW, EV, and SVY performed the experiments. SHK, NE, CW, EV, and SVY analyzed the data. SVY contributed reagents/materials/analysis tools. EV and SVY wrote the paper.

References

  1. 1. Nei M, Kumar S (2000) Molecular evolution and phylogenetics. New York: Oxford University Press.
  2. 2. Li WH, Ellsworth DL, Krushkal J, Chang BHJ, Hewett-Emmett D (1996) Rates of nucleotide substitution in primates and rodents and the generation-time effect hypothesis. Mol Phylogenet Evol 5: 182–187.
  3. 3. Li WH (1997) Molecular evolution. Sunderland (Massachusetts): Sinauer.
  4. 4. Yi S, Ellsworth DL, Li WH (2002) Slow molecular clocks in Old World monkeys, apes, and humans. Mol Biol Evol 19: 2191–2198.
  5. 5. Laird CD, McConaughy BL, McCarthy BJ (1969) Rate of fixation of nucleotide substitutions in evolution. Nature 224: 149–154.
  6. 6. Kumar S, Subramanian S (2002) Mutation rates in mammalian genomes. Proc Natl Acad Sci U S A 99: 803–808.
  7. 7. Kumar S (2005) Molecular clocks: Four decades of evolution. Nat Rev Genet 6: 654–662.
  8. 8. Easteal S, Collet C (1994) Consistent variation in amino-acid substitution rate, despite uniformity of mutation rate: Protein evolution in mammals is not neutral. Mol Biol Evol 11: 643–647.
  9. 9. Goodman M (1961) The role of immunologic differences in the phyletic development of human behavior. Hum Biol 33: 131–162.
  10. 10. Nachman MW, Crowell SL (2000) Estimate of the mutation rate per nucleotide in humans. Genetics 156: 297–304.
  11. 11. The Chimpanzee Sequencing and Analysis Consortium (2005) Initial sequence of the chimpanzee genome and comparison with the human genome. Nature 437: 69–87.
  12. 12. Meunier J, Khelifi A, Navratil V, Duret L (2005) Homology-dependent methylation in primate repetitive DNA. Proc Natl Acad Sci U S A 102: 5471–5476.
  13. 13. Robertson KD, Wolffe AP (2000) DNA methylation in health and disease. Nat Rev Genet 1: 11–19.
  14. 14. Li WH, Yi S, Makova K (2002) Male-driven evolution. Curr Opin Genet Dev 12: 650–656.
  15. 15. Jones PA, Laird PW (1999) Cancer epigenetics comes of age. Nat Genet 21: 163–167.
  16. 16. Keshet I, Schlesinger Y, Farkash S, Rand E, Hecht M, et al. (2006) Evidence for an instructive mechanism of de novo methylation in cancer cells. Nat Genet 38: 149–157.
  17. 17. Bird A (1980) DNA methylation and the frequency of CpG in animal DNA. Nucleic Acids Res 8: 1499–1504.
  18. 18. Duncan BK, Miller JH (1980) Mutagenic deamination of cytosine residues in DNA. Nature 287: 560–561.
  19. 19. Bird A (2002) DNA methylation patterns and epigenetic memory. Genes Dev 16: 6–21.
  20. 20. Li E (2002) Chromatin modification and epigenetic reprogramming in mammalian development. Nat Rev Genet 3: 662–673.
  21. 21. Hwang DG, Green P (2004) Inaugural Article: Bayesian Markov chain Monte Carlo sequence analysis reveals varying neutral substitution patterns in mammalian evolution. Proc Natl Acad Sci U S A 101: 13994–14001.
  22. 22. Taylor J, Tyekucheva S, Zody M, Chiaromonte F, Makova KD (2006) Strong and weak male mutation bias at different sites in the primate genomes: Insights from the human-chimpanzee comparison. Mol Biol Evol 23: 565–573.
  23. 23. Meunier J, Duret L (2004) Recombination drives the evolution of GC-content in the human genome. Mol Biol Evol 21: 984–990.
  24. 24. Gage TB (1998) The comparative demography of primates: with some comments on the evolution of life histories. Annu Rev Anthropol 27: 197–221.
  25. 25. Goodman M (1962) Evolution of the immunologic species specificity of human serum proteins. Hum Biol 34: 104–150.
  26. 26. Wu CI, Li WH (1985) Evidence for higher rates of nucleotide substitution in rodents than in man. Proc Natl Acad Sci U S A 82: 1741–1745.
  27. 27. Steiper ME, Young NM, Sukrarna TY (2004) Genomic data support the hominoid slowdown and an early Oligocene estimate for the hominoid-cercopithecoid divergence. Proc Nat Acad Sci U S A 101: 17021–17026.
  28. 28. Saitou N, Nei M (1987) The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol 4: 406–425.
  29. 29. Delson E, Tattersall I, Van Couvering JA, Brooks AS (2000) Encyclopedia of human evolution and prehistory. 2nd edition. New York: Garland. pp. 166–171.
  30. 30. Brunet M, Guy F, Pilbeam D, Mackaye HT, Likius A, et al. (2002) A new hominid from the Upper Miocene of Chad, Central Africa. Nature 418: 145–151.
  31. 31. Brunet M, Guy F, Pilbeam D, Lieberman DE, Likius A, et al. (2005) New material of the earliest hominid from the Upper Miocene of Chad. Nature 434: 752–755.
  32. 32. Ebersberger I, Metzler D, Schwartz C, Pääbo S (2002) Genomewide comparison of DNA sequences between human and chimpanzees. Am J Hum Genet 70: 1490–1497.
  33. 33. Makova KD, Li WH (2002) Strong male-driven evolution of DNA sequences in humans and apes. Nature 416: 624–626.
  34. 34. Elango N, Thomas JW, Program NCS, Yi S (2006) Variable molecular clocks in hominoids. Proc Nat Acad Sci U S A 103: 1370–1375.
  35. 35. Chen FC, Li WH (2001) Genomic divergence between humans and other hominoids and the effective population size of the common ancestor of humans and chimpanzees. Am J Hum Genet 68: 444–456.
  36. 36. Wall JD (2003) Estimating ancestral population sizes and divergence times. Genetics 163: 395–404.
  37. 37. Rogers J, Kidd KK (1996) Nucleotide polymorphism, effective population size, and dispersal distances in the yellow baboons (Papio hamadryas cynocephalus) of Mikumi National Park, Tanzania. Am J Primatol 38: 157–168.
  38. 38. Wall JD, Frisse LA, Hudson RR, Di Rienzo A (2003) Comparative linkage-disequilibrium analysis of the beta-globin hotspot in primates. Am J Hum Genet 73: 1330–1340.
  39. 39. Subramanian S, Kumar S (2003) Neutral substitutions occur at a faster rate in exons than in noncoding DNA in primate genomes. Genome Res 13: 838–844.
  40. 40. Keightley PD, Lercher MJ, Eyre-Walker A (2005) Evidence for widespread degradation of gene control regions in hominoid genomes. PLoS Biol 3: e42.. DOI: https://doi.org/10.1371/journal.pbio.0030042 .
  41. 41. Hellmann I, Ebersberger I, Ptak SE, Paabo S, Przeworski M (2003) A neutral explanation for the correlation of diversity with recombination rates in humans. Am J Hum Genet 72: 1527–1535.
  42. 42. Siepel A, Haussler D (2004) Phylogenetic estimation of context-dependent substitution rates by maximum likelihood. Mol Biol Evol 21: 468–488.
  43. 43. Arndt PF, Burge CB, Hwa T (2002) DNA sequence evolution with neighbor-dependent mutation. 6th Annu Int Conf Comp Biol. pp. 32–38.
  44. 44. Taudien S, Ebersberger I, Glöckner G, Platzer M (2006) Should the draft chimpanzee sequence be finished? Trends Genet 22: 122–125.
  45. 45. The ENCODE Project Consortium (2004) The ENCODE (ENCyclopedia Of DNA Elements) Project. Science 306: 636–640.
  46. 46. International Human Genome Sequencing Consortium (2001) Initial sequencing and analysis of the human genome. Nature 409: 860–921.
  47. 47. Altschul DA, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215: 403–410.
  48. 48. Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D (2003) Evolution's cauldron: Duplication, deletion and rearrangement in the mouse and human genomes. Proc Nat Acad Sci U S A 100: 11484–11489.
  49. 49. Kent WJ (2002) BLAT: The BLAST-like alignment tool. Genome Res 12: 656–664.
  50. 50. Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL (2006) GenBank. Nucleic Acids Res 34: D16–D20.
  51. 51. Duret L, Mouchiroud D, Gouy M (1994) HOVERGEN: A database of homologous vertebrate genes. Nucleic Acids Res 22: 2360–2365.
  52. 52. Birney E, Andrews D, Caccamo M, Chen Y, Clarke L, et al. (2006) Ensembl 2006. Nucleic Acids Res 34: D556–D561.
  53. 53. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.
  54. 54. Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, et al. (2002) The Bioperl toolkit: Perl modules for the life sciences. Genome Res 12: 1611–1618.
  55. 55. Li L, Stoeckert CJJ, Roos DS (2003) OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res 13: 2178–2189.
  56. 56. Takai D, Jones PA (2002) Comprehensive analysis of CpG islands in human chromosomes 21 and 22. Proc Natl Acad Sci U S A 99: 3740–3745.
  57. 57. Smit AFA, Hubely R, Green P (2004) RepeatMasker Open-3.0. Available: http://www.repeatmasker.org. Accessed 6 September 2006.
  58. 58. Kimura M (1980) A simple method for estimating evolutionary rate of base substitution through comparative studies of nucleotide sequences. J Mol Evol 16: 111–120.