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

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

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.     

Evolutionary conservation of 5' upstream sequence of nine genes between human and great apes

Nucleotide sequences of nine 5' upstream gene regions for human, chimpanzee, gorilla, and orangutan were determined. We estimated nucleotide differences (d) for each region between human and great apes. The overall d was 0.027 (ranged from 0.004 to 0.052). Rates of nucleotide substitution were estimated by using d and divergence times of human, chimpanzee, gorilla, and orangutan. The overall rate of nucleotide substitution between human and other hominoids was estimated to be 0.52-0.85 x 10(-9). This rate in 5' upstream regions was lower than that of synonymous sites, suggesting that 5' upstream regions have evolved under some functional constraints. Because lower rates have been reported for coding sequences in primates compared to rodents, we also estimated the rate (1.17-1.76 x 10(-9)) of nucleotide substitutions for the corresponding 5' upstream regions in rodents (mouse/rat comparison). Thus the primate rate was lower than rodent rate also for the 5' upstream regions.     

Molecular systematics and paleobiogeography of the South American sigmodontine rodents [published erratum appears in Mol Biol Evol 1998 Feb;15(2):224]

The murid rodent subfamily Sigmodontinae contains 79 genera which are distributed throughout the New World. The time of arrival of the first sigmodontines in South America and the estimated divergence time(s) of the different lineages of South American sigmodontines have been controversial due to the lack of a good fossil record and the immense number of extant species. The "early-arrival hypothesis" states that the sigmodontines must have arrived in South America no later than the early Miocene, at least 20 MYA, in order to account for their vast present-day diversity, whereas the "late-arrival hypothesis" includes the sigmodontines as part of the Plio-Pleistocene Great American Interchange, which occurred approximately 3.5 MYA. The phylogenetic relationships among 33 of these genera were reconstructed using mitochondrial DNA (mtDNA) sequence data from the ND3, ND4L, arginine tRNA, and ND4 genes, which we show to be evolving at the same rate. A molecular clock was calibrated for these genes using published fossil dates, and the genetic distances were estimated from the DNA sequences in this study. The molecular clock was used to estimate the dates of the South American sigmodontine origin and the main sigmodontine radiation in order to evaluate the "early-" and "late-arrival" scenarios. We estimate the time of the sigmodontine invasion of South America as between approximately 5 and 9 MYA, supporting neither of the scenarios but suggesting two possible models in which the invading lineage was either (1) ancestral to the oryzomyines, akodonts, and phyllotines or (2) ancestral to the akodonts and phyllotines and accompanied by the oryzomyines. The sigmodontine invasion of South America provides an example of the advantage afforded to a lineage by the fortuitous invasion of a previously unexploited habitat, in this case an entire continent.      

Ancient Origin of the Null Allele se 428 of the Human ABO-Secretor Locus (FUT2)

In human populations, a null allele having several nucleotide differences from the wild-type allele is segregating at the FUT2 locus (the ABO-Secretor locus) encoding α(1,2)fucosyltransferase. To estimate the age of the most recent common ancestor (MRCA) of these two alleles, we sequenced FUT2 homologues from chimpanzee, gorilla, orangutan, and green monkey. Since we did not detect acceleration or any heterogeneity in the substitution rate at this locus among these species, the age of the MRCA was estimated to be around 3 MYA, assuming the divergence time of human and chimpanzee to be 5 MYA. We developed a simple test to examine whether or not the old age of the MRCA of the FUT2 is consistent with that expected for two divergent neutral alleles sampled from a random mating population. An application of the test to the data at FUT2 indicated that the age of the MRCA is too old to be explained by the simple neutral assumptions, although our test depends on accurate estimation of the divergence time of human and chimpanzee in units of twice the human population size. Various possibilities including balancing selection are discussed to explain this old age of the MRCA. Received: 9 May 1999 / Accepted: 20 September 1999     

A mitogenomic timescale for birds detects variable phylogenetic rates of molecular evolution and refutes the standard molecular clock

Current understanding of the diversification of birds is hindered by their incomplete fossil record and uncertainty in phylogenetic relationships and phylogenetic rates of molecular evolution. Here we performed the first comprehensive analysis of mitogenomic data of 48 vertebrates, including 35 birds, to derive a Bayesian timescale for avian evolution and to estimate rates of DNA evolution. Our approach used multiple fossil time constraints scattered throughout the phylogenetic tree and accounts for uncertainties in time constraints, branch lengths, and heterogeneity of rates of DNA evolution. We estimated that the major vertebrate lineages originated in the Permian; the 95% credible intervals of our estimated ages of the origin of archosaurs (258 MYA), the amniote-amphibian split (356 MYA), and the archosaur-lizard divergence (278 MYA) bracket estimates from the fossil record. The origin of modern orders of birds was estimated to have occurred throughout the Cretaceous beginning about 139 MYA, arguing against a cataclysmic extinction of lineages at the Cretaceous/Tertiary boundary. We identified fossils that are useful as time constraints within vertebrates. Our timescale reveals that rates of molecular evolution vary across genes and among taxa through time, thereby refuting the widely used mitogenomic or cytochrome b molecular clock in birds. Moreover, the 5-Myr divergence time assumed between 2 genera of geese (Branta and Anser) to originally calibrate the standard mitochondrial clock rate of 0.01 substitutions per site per lineage per Myr (s/s/l/Myr) in birds was shown to be underestimated by about 9.5 Myr. Phylogenetic rates in birds vary between 0.0009 and 0.012 s/s/l/Myr, indicating that many phylogenetic splits among avian taxa also have been underestimated and need to be revised. We found no support for the hypothesis that the molecular clock in birds "ticks" according to a constant rate of substitution per unit of mass-specific metabolic energy rather than per unit of time, as recently suggested. Our analysis advances knowledge of rates of DNA evolution across birds and other vertebrates and will, therefore, aid comparative biology studies that seek to infer the origin and timing of major adaptive shifts in vertebrates.     

A study of the phylogeny of Brassica rapa,B. nigra,Raphanus sativus,and their related genera using noncoding regions of chloroplast DNA

There are two evolutionary lineages in the genus Brassica: the rapa/oleracea lineage and the nigra lineage. Using nuclear DNA sequences such as the intergenic spacer between 5S rRNA genes and the internal transcribed spacer between 18S and 25S rRNA genes, we and others had previously demonstrated that Raphanus sativus is closely related to the nigra lineage. In the present study, we sequenced the chloroplast noncoding region between trnT and trnF and that between trnD and trnT in seven species and showed that R. sativus is more closely related to the rapa/oleracea lineage than to the nigra lineage. The conflicting results from nuclear DNA and chloroplast DNA support the hypothesis that Raphanus was derived from a hybridization between the rapa/oleracea and the nigra lineages. We estimated the date of this hybridization event to be 60% of the divergence time between the two Brassica lineages. In addition, the pattern and rate of nucleotide substitution were studied. There were more transversions than transitions in these noncoding regions, which have a high AT content. Furthermore, the proportion of transversions among the substitutions at a site increases with increasing A + T content of its two adjacent nucleotides. An influence of immediate 5(') pyrimidine on substitution pattern is also observed when both adjacent bases in the two DNA strands are A or T. The rate of nucleotide substitution in the trnL group I intron is only about one third of the rate in the nearby intergenic spacers in the trnT-trnF fragment. The rate of nucleotide substitution in the rapa/oleracea lineage is at least 1.5 times that in the nigra lineage.     

Estimation of primate speciation dates using local molecular clocks

Protein-coding genes of the mitochondrial genomes from 31 mammalian species were analyzed to estimate the speciation dates within primates and also between rats and mice. Three calibration points were used based on paleontological data: one at 20-25 MYA for the hominoid/cercopithecoid divergence, one at 53-57 MYA for the cetacean/artiodactyl divergence, and the third at 110-130 MYA for the metatherian/eutherian divergence. Both the nucleotide and the amino acid sequences were analyzed, producing conflicting results. The global molecular clock was clearly violated for both the nucleotide and the amino acid data. Models of local clocks were implemented using maximum likelihood, allowing different evolutionary rates for some lineages while assuming rate constancy in others. Surprisingly, the highly divergent third codon positions appeared to contain phylogenetic information and produced more sensible estimates of primate divergence dates than did the amino acid sequences. Estimated dates varied considerably depending on the data type, the calibration point, and the substitution model but differed little among the four tree topologies used. We conclude that the calibration derived from the primate fossil record is too recent to be reliable; we also point out a number of problems in date estimation when the molecular clock does not hold. Despite these obstacles, we derived estimates of primate divergence dates that were well supported by the data and were generally consistent with the paleontological record. Estimation of the mouse-rat divergence date, however, was problematic.     

Investigation of loss and gain of introns in the compact genomes of pufferfishes (Fugu and Tetraodon)

We have investigated intron evolution in the compact genomes of 2 closely related species of pufferfishes, Fugu rubripes and Tetraodon nigroviridis, that diverged about 32 million years ago (MYA). Analysis of 148,028 aligned intron positions in 13,547 gene pairs using human as an outgroup identified 57 and 24 intron losses in Tetraodon and fugu lineages, respectively, and no gain in either lineage. For comparison, we analyzed 144,545 intron positions in 12,866 orthologous pairs of genes in human and mouse that diverged about 61 MYA using dog as an outgroup and identified 51 intron losses in mouse and 3 losses in human and no gain. The rate of intron loss in Tetraodon is higher than that in fugu, mouse, and human but lower than the previous estimates for other eukaryotes. The introns lost in pufferfishes and mammals are significantly shorter than the mean size of introns in the genome. One intron deleted in fugu and another in Tetraodon have left behind 6 and 3 nucleotides, respectively, suggesting that they were lost due to genomic deletions. Such losses of introns are likely to be the result of a higher rate of DNA deletions experienced by the genomes of pufferfishes compared with mammals. The shorter generation time of Tetraodon compared with fugu, and the rich diversity and higher activity of transposable elements in pufferfishes compared with mammals, may be responsible for the higher rate of intron loss in Tetraodon. Our findings indicate that overall very little intron turnover has occurred in pufferfishes and mammals during recent evolution and that intron gain is an extremely rare event in vertebrate evolution.     

Extensive mitochondrial DNA transfer in a rapidly evolving rodent has been mediated by independent insertion events and by duplications

Mitochondrial DNA translocations to the nucleus (numt pseudogenes) are pervasive among eukaryotes, but copy number within the nuclear genome varies widely among taxa. As an increasing number of genomes are sequenced in their entirety, the origins, transfer mechanisms and insertion sites of numts are slowly being characterized. We investigated mitochondrial transfers within a genetically diverse rodent lineage and here report 15 numts totaling 21.8 kb that are harbored within the nuclear genome of the vole Microtus rossiaemeridionalis. The 15 numts total 21.8 kb and range from 0.39 to over 3.0 kb in length. Phylogenetic analyses revealed that these numts resulted from three independent insertions to the nucleus, two of which were followed by subsequent nuclear duplication events. The dates of the two translocations that led to subsequent duplications were estimated at 1.97 and 1.19 MYA, which coincide with the origin and radiation of the genus Microtus. Numt sequence data from five Microtus species were used to estimate an average rate of nucleotide substitution as 2.6x10(-8) subs/site/yr. This substitution rate is higher than in many other mammals, but is concordant with the elevated rate of mtDNA substitution in this lineage. Our data suggest that numt translocation in Microtus is more extensive than in either Mus or in Rattus, consistent with the elevated rate of speciation, karyotypic rearrangement, and mitochondrial DNA evolution in Microtus.     

Reexamination of the African hominoid trichotomy with additional sequences from the primate beta-globin gene cluster.

Additional DNA sequence information from a range of primates, including 13.7 kb from pygmy chimpanzee (Pan paniscus), was added to data sets of beta-globin gene cluster sequence alignments that span the gamma 1, gamma 2, and psi eta loci and their flanking and intergenic regions. This enlarged body of data was used to address the issue of whether the ancestral separations of gorilla, chimpanzee, and human lineages resulted from only one trichotomous branching or from two dichotomous branching events. The degree of divergence, corrected for superimposed substitutions, seen in the beta-globin gene cluster between human alleles is about a third to a half that observed between two species of chimpanzee and about a fourth that between human and chimpanzee. The divergence either between chimpanzee and gorilla or between human and gorilla is slightly greater than that between human and chimpanzee, suggesting that the ancestral separations resulted from two closely spaced dichotomous branchings. Maximum parsimony analysis further strengthened the evidence that humans and chimpanzees share the longest common ancestry. Support for this human-chimpanzee clade is statistically significant at P = 0.002 over a human-gorilla clade or a chimpanzee-gorilla clade. An analysis of expected and observed homoplasy revealed that the number of sequence changes uniquely shared by human and chimpanzee lineages is too large to be attributed to homoplasy. Molecular clock calculations that accommodated lineage variations in rates of molecular evolution yielded hominoid branching times that ranged from 17-19 million years ago (MYA) for the separation of gibbon from the other hominoids to 5-7 MYA for the separation of chimpanzees from humans. Based on the relatively late dates and mounting corroborative evidence from unlinked nuclear genes and mitochondrial DNA for the close sister grouping of humans and chimpanzees, a cladistic classification would place all apes and humans in the same family. Within this family, gibbons would be placed in one subfamily and all other extant hominoids in another subfamily. The later subfamily would be divided into a tribe for orangutans and another tribe for gorillas, chimpanzees, and humans. Finally, gorillas would be placed in one subtribe with chimpanzees and humans in another, although this last division is not as strongly supported as the other divisions.     

Genomic divergences among cattle,dog and human estimated from large-scale alignments of genomic sequences

Background

Approximately 11 Mb of finished high quality genomic sequences were sampled from cattle, dog and human to estimate genomic divergences and their regional variation among these lineages.

Results

Optimal three-way multi-species global sequence alignments for 84 cattle clones or loci (each >50 kb of genomic sequence) were constructed using the human and dog genome assemblies as references. Genomic divergences and substitution rates were examined for each clone and for various sequence classes under different functional constraints. Analysis of these alignments revealed that the overall genomic divergences are relatively constant (0.32–0.37 change/site) for pairwise comparisons among cattle, dog and human; however substitution rates vary across genomic regions and among different sequence classes. A neutral mutation rate (2.0–2.2 × 10(-9) change/site/year) was derived from ancestral repetitive sequences, whereas the substitution rate in coding sequences (1.1 × 10(-9) change/site/year) was approximately half of the overall rate (1.9–2.0 × 10(-9) change/site/year). Relative rate tests also indicated that cattle have a significantly faster rate of substitution as compared to dog and that this difference is about 6%.

Conclusion

This analysis provides a large-scale and unbiased assessment of genomic divergences and regional variation of substitution rates among cattle, dog and human. It is expected that these data will serve as a baseline for future mammalian molecular evolution studies.

     

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.     

Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for pinus

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.     

Neuroglobins from the zebrafish Danio rerio and the pufferfish Tetraodon nigroviridis

Neuroglobin is a recently discovered respiratory, porphyrin-containing protein that is expressed in the brain of mouse and man. Here we show that neuroglobin is also present in the teleost fish. Complete cDNA sequences are reported from the pufferfish Tetraodon nigroviridis and the zebrafish Danio rerio. In addition, the neuroglobin gene of T. nigroviridis was sequenced, demonstrating the conservation of the B12.2, E11.0 and G7.0 introns plus the presence of an additional intron in the 5' noncoding region. The fish neuroglobins each comprise 159 amino acids and are 84.3% identical. Phylogenetic analyses show a basal position of the neuroglobins within the metazoan globin tree. An enhanced amino acid substitution rate was estimated for the fish neuroglobins ( approximately 0.93 x 10(-9) amino acid substitutions per site and year) compared with their mammalian proteins ( approximately 0.39 x 10(-9) replacements per site and year).     

Molecular evolution of the psi eta-globin gene locus: gibbon phylogeny and the hominoid slowdown

An 8.4-kb genomic region spanning both the psi eta-globin gene locus and flanking DNA was sequenced from the common gibbon (Hylobates lar). In addition, sequencing of the entire orthologous region from galago (Galago crassicaudatus) was completed. The gibbon and galago sequences, along with published orthologous sequences from 10 other species, were aligned. These noncoding nucleotide sequences represented four human alleles, four apes (chimpanzee, gorilla, organgutan, and gibbon), an Old World monkey (rhesus monkey), two New World monkeys (spider and owl monkeys), tarsier, two strepsirhines (galago and lemur), and goat. Divergence and maximum parsimony analyses of the psi eta genomic region first groups humans and chimpanzees and then, at progressively more ancient branch points, successively joins gorillas, orangutans, gibbons, Old World monkeys, New World monkeys, tarsiers, and strepsirhines (the lemuriform-lorisiform branch of primates). This cladistic pattern supports the taxonomic grouping of all extant hominoids into family Hominidae, the division of Hominidae into subfamilies Hylobatinae (gibbons) and Homininae, the division of Homininae into tribes Pongini (orangutans) and Hominini, and the division of Hominini into subtribes Gorillina (gorillas) and Hominina (chimpanzees and humans). The additional gibbon and galago sequence data provide further support for the occurrence of a graded evolutionary-rate slowdown in the descent of simian primates, with the slowing rate being more pronounced in the great-ape and human lineages than in the gibbon or monkey lineages. A comparison of global versus local molecular clocks reveals that local clock predictions, when focused on a specific number of species within a narrow time frame, provide a more accurate estimate of divergence dates than do those of global clocks.     

The evolutionary history of human and chimpanzee Y-chromosome gene loss

Recent studies have suggested that gene gain and loss may contribute significantly to the divergence between humans and chimpanzees. Initial comparisons of the human and chimpanzee Y-chromosomes indicate that chimpanzees have a disproportionate loss of Y-chromosome genes, which may have implications for the adaptive evolution of sex-specific as well as reproductive traits, especially because one of the genes lost in chimpanzees is critically involved in spermatogenesis in humans. Here we have characterized Y-chromosome sequences in gorilla, bonobo, and several chimpanzee subspecies for 7 chimpanzee gene-disruptive mutations. Our analyses show that 6 of these gene-disruptive mutations predate chimpanzee-bonobo divergence at approximately 1.8 MYA, which indicates significant Y-chromosome change in the chimpanzee lineage relatively early in the evolutionary divergence of humans and chimpanzees.     

真核生物DNA非编码区的组分分析

在全基因组水平上,用直方图、混沌表示灰度图、距离差异度和信息熵差异度四种方法,研究了拟南芥、线虫、果蝇的DNA内含子、基因间隔区DNA、外显子三种区域的核苷酸短序列组分及组分复杂度.结果表明：a.不同基因组之间,不管基因数目多少,用4种方法得到的外显子部分其组分复杂度都比较接近,而非编码区部分的组分复杂度却很大.这一点定量地说明了物种之间的复杂程度,主要不体现在编码区部分,而体现在非编码区部分.b.同一基因组中,内含子的核苷酸短序列组分复杂度都是相似的,外显子和intergenic DNA部分的组分复杂度也是相似的.c.内含子和intergenic DNA在转录、剪切、二级结构等方面有很大的不同,但它们在核苷酸短序列组分上的差异却很小,说明内含子和intergenic DNA在转录、剪切、二级结构上的不同并不通过核苷酸短序列组分来进行限制.