Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “divergence time”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8Linked to original sources

Site stripping based on likelihood ratio reduction is a useful tool to evaluate the impact of non-clock-like behavior on viral phylogenetic reconstructions.

The site stripping for clock detection procedure was implemented in the recently developed maximum likelihood framework for estimating evolutionary rates and divergence times in measurably evolving populations. The method was used to investigate the effect of rate variability on estimating divergence times in non-clock-like trees for human immunodeficiency viruses and hepatitis C viruses. We validate our approach by comparing dated coalescent nodes in molecular phylogenies with known dates of transmission. Our method was able to rapidly recover clock-like behavior and to indicate the presence and direction of a bias when estimates of divergence times using the unstripped data were flawed.

Algorithms↗

Recent African origin of modern humans revealed by complete sequences of hominoid mitochondrial DNAs.

We analyzed the complete mitochondrial DNA (mtDNA) sequences of three humans (African, European, and Japanese), three African apes (common and pygmy chimpanzees, and gorilla), and one orangutan in an attempt to estimate most accurately the substitution rates and divergence times of hominoid mtDNAs. Nonsynonymous substitutions and substitutions in RNA genes have accumulated with an approximately clock-like regularity. From these substitutions and under the assumption that the orangutan and African apes diverged 13 million years ago, we obtained a divergence time for humans and chimpanzees of 4.9 million years. This divergence time permitted calibration of the synonymous substitution rate (3.89 x 10(-8)/site per year). To obtain the substitution rate in the displacement (D)-loop region, we compared the three human mtDNAs and measured the relative abundance of substitutions in the D-loop region and at synonymous sites. The estimated substitution rate in the D-loop region was 7.00 x 10(-8)/site per year. Using both synonymous and D-loop substitutions, we inferred the age of the last common ancestor of the human mtDNAs as 143,000 +/- 18,000 years. The shallow ancestry of human mtDNAs, together with the observation that the African sequence is the most diverged among humans, strongly supports the recent African origin of modern humans, Homo sapiens sapiens.

Africa↗

The eastern Asian and eastern and western North American floristic disjunction: congruent phylogenetic patterns in seven diverse genera.

One of the most remarkable examples of intercontinental disjunction of the North Temperate Flora involves eastern Asia and eastern and western North America. Although there has been considerable interest in this phytogeographic pattern for over 150 years (e.g., Gray, 1859; Li, 1952; Graham, 1972; Boufford and Spongberg, 1983; Wu, 1983; Tiffney, 1985a, 1985b), relationships among taxa displaying the disjunction remain obscure. Understanding phylogenetic relationships is, however, a prerequisite for historical biogeographic analyses of this distributional pattern. To understand better the relationships of taxa displaying this intercontinental disjunction, phylogenetic analyses were conducted using a variety of DNA data sets for species of four genera (Cornus, Boykinia, Tiarella, and Trautvetteria) that occur in eastern Asia, eastern North America, and western North America. An area cladogram was constructed for each of the four genera, all of which show a similar pattern of relationship: the eastern Asian species are sister to all North American species. An identical phylogenetic pattern is also found in three other taxa exhibiting this disjunction (Aralia sect. Aralia, Calycanthus, and Adiantum pedatum). The congruent phylogenetic pattern found in these seven diverse genera raises the possibility of a common origin of the eastern Asia, eastern and western North America disjunction. The data are in agreement with the long-standing hypothesis that this well-known floristic disjunction represents the fragmentation of a once continuous Mixed Mesophytic forest community and suggest that the disjunction may have involved only two major vicariance events: an initial split between Eurasia and North America, followed by the isolation of floras between eastern and western North America. However, congruence between phylogenies and geographic distributions does not necessarily indicate an identical phytogeographic history. Taxa exhibiting the same phylogenetic pattern may have originated at different geological times. Analysis of divergence times using the molecular clock indicates that species of Cornus, Boykinia, and Calycanthus may have diverged at different geological times, suggesting that the floristic disjunction involving eastern Asia and North America may not be simple; it may have involved multiple historical events at very different geological times in different genera.

Asia, Eastern↗

Evolutionary models for insertions and deletions in a probabilistic modeling framework.

BACKGROUND: Probabilistic models for sequence comparison (such as hidden Markov models and pair hidden Markov models for proteins and mRNAs, or their context-free grammar counterparts for structural RNAs) often assume a fixed degree of divergence. Ideally we would like these models to be conditional on evolutionary divergence time. Probabilistic models of substitution events are well established, but there has not been a completely satisfactory theoretical framework for modeling insertion and deletion events. RESULTS: I have developed a method for extending standard Markov substitution models to include gap characters, and another method for the evolution of state transition probabilities in a probabilistic model. These methods use instantaneous rate matrices in a way that is more general than those used for substitution processes, and are sufficient to provide time-dependent models for standard linear and affine gap penalties, respectively. Given a probabilistic model, we can make all of its emission probabilities (including gap characters) and all its transition probabilities conditional on a chosen divergence time. To do this, we only need to know the parameters of the model at one particular divergence time instance, as well as the parameters of the model at the two extremes of zero and infinite divergence. I have implemented these methods in a new generation of the RNA genefinder QRNA (eQRNA). CONCLUSION: These methods can be applied to incorporate evolutionary models of insertions and deletions into any hidden Markov model or stochastic context-free grammar, in a pair or profile form, for sequence modeling.

Algorithms↗

The "inverse relationship between evolutionary rate and age of mammalian genes" is an artifact of increased genetic distance with rate of evolution and time of divergence.

It has recently been claimed that older genes tend to evolve more slowly than newer ones (Alba and Castresana 2005). By simulation of genes of equal age, we show that the inverse correlation between age and rate is an artifact caused by our inability to detect homology when evolutionary distances are large. Since evolutionary distance increases with time of divergence and rate of evolution, homologs of fast-evolving genes are frequently undetected in distantly related taxa and are, hence, misclassified as "new." This misclassification causes the mean genetic distance of 'new' genes to be overestimated and the mean genetic distance of "old" genes to be underestimated.

Animals↗

The origin and evolution of porcine reproductive and respiratory syndrome viruses.

Porcine reproductive and respiratory syndrome viruses (PRRSV) are divided into North American and European types, which show about 40% difference in their amino acid sequences. The divergence time of these two types has been estimated to be about 1980 from epidemiological data. This suggested that PRRSV have evolved at a higher evolutionary rate (order of 10(-2)/site/year) compared with other RNA viruses of 10(-3) to 10(-5)/site/year. Here, to test the evolutionary history of PRRSV speculated by the epidemiological background, we estimated the divergence time and evolutionary rate of PRRSV with molecular evolutionary analysis. Estimated divergence time (1972-1988) corresponded well to that estimated by the epidemiological data, and the evolutionary rate (4.71-9.8) x 10(-2) of PRRSV was indeed the highest among RNA viruses so far reported. Furthermore, we inferred important sites for the adaptation in order to examine how PRRSV have adapted to swine since they emerged. The adaptive sites were located not only in the epitopes related to immunity but also in the transmembrane regions including a signal peptide. In particular, the adaptive sites in the transmembrane regions were considered to affect compatibility to the host cell membrane. We conclude that PRRSV were transmitted from another host species to swine in about 1980 and have adapted to swine by altering the transmembrane regions.

Amino Acid Sequence↗

On morphological clocks and paleophylogeography: towards a timescale for Sorex hybrid zones.

Phylogeography--the study of within-species phylogenetic and geographic divergence--has been primarily the domain of molecular evolutionists because molecular markers record population structure on smaller scales than do traditional morphological traits. But when geometric morphometrics are combined with distance-based phylogenetics molar shape divergence appears to record population-level phylogeny, a fact that allows extant and fossil populations to be combined in a single phylogeographic study. The European Sorex araneus complex--a genetically complicated group composed of multiple karyotypic races and species--illustrates the principle. The phylogeographic patterns revealed by molar shape broadly agree with scenarios based on molecular data and circumstantial evidence. Importantly, the inclusion of fossil samples of known age allows minimum divergence times to be inferred. Some races of S. araneus may have diverged more than 120,000 years ago, but others may have diverged less than 14,000. Supporting evidence that molar shape can be used to reconstruct phylogeographic relationships comes from strong correlations between molar shape distances and both phylogenetic divergence time and cytochrome b sequence divergence in datasets where these variables are known independently (fossil carnivorans from a well-constrained stratigraphic setting and shrew species of the genus Sorex, respectively). However, molar shape may have a 'saturation point' beyond which it is not applicable.

Animals↗

A genomic timescale of prokaryote evolution: insights into the origin of methanogenesis, phototrophy, and the colonization of land.

BACKGROUND: The timescale of prokaryote evolution has been difficult to reconstruct because of a limited fossil record and complexities associated with molecular clocks and deep divergences. However, the relatively large number of genome sequences currently available has provided a better opportunity to control for potential biases such as horizontal gene transfer and rate differences among lineages. We assembled a data set of sequences from 32 proteins (approximately 7600 amino acids) common to 72 species and estimated phylogenetic relationships and divergence times with a local clock method. RESULTS: Our phylogenetic results support most of the currently recognized higher-level groupings of prokaryotes. Of particular interest is a well-supported group of three major lineages of eubacteria (Actinobacteria, Deinococcus, and Cyanobacteria) that we call Terrabacteria and associate with an early colonization of land. Divergence time estimates for the major groups of eubacteria are between 2.5-3.2 billion years ago (Ga) while those for archaebacteria are mostly between 3.1-4.1 Ga. The time estimates suggest a Hadean origin of life (prior to 4.1 Ga), an early origin of methanogenesis (3.8-4.1 Ga), an origin of anaerobic methanotrophy after 3.1 Ga, an origin of phototrophy prior to 3.2 Ga, an early colonization of land 2.8-3.1 Ga, and an origin of aerobic methanotrophy 2.5-2.8 Ga. CONCLUSIONS: Our early time estimates for methanogenesis support the consideration of methane, in addition to carbon dioxide, as a greenhouse gas responsible for the early warming of the Earths' surface. Our divergence times for the origin of anaerobic methanotrophy are compatible with highly depleted carbon isotopic values found in rocks dated 2.8-2.6 Ga. An early origin of phototrophy is consistent with the earliest bacterial mats and structures identified as stromatolites, but a 2.6 Ga origin of cyanobacteria suggests that those Archean structures, if biologically produced, were made by anoxygenic photosynthesizers. The resistance to desiccation of Terrabacteria and their elaboration of photoprotective compounds suggests that the common ancestor of this group inhabited land. If true, then oxygenic photosynthesis may owe its origin to terrestrial adaptations.

Bacteria↗

Assessing concordance of fossil calibration points in molecular clock studies: an example using turtles.

Although still controversial, estimation of divergence times using molecular data has emerged as a powerful tool to examine the tempo and mode of evolutionary change. Two primary obstacles in improving the accuracy of molecular dating are heterogeneity in DNA substitution rates and accuracy of the fossil record as calibration points. Recent methodological advances have provided powerful methods that estimate relative divergence times in the face of heterogeneity of nucleotide substitution rates among lineages. However, relatively little attention has focused on the accuracy of fossil calibration points that allow one to translate relative divergence times into absolute time. We present a new cross-validation method that identifies inconsistent fossils when multiple fossil calibrations are available for a clade and apply our method to a molecular phylogeny of living turtles with fossil calibration times for 17 of the 22 internal nodes in the tree. Our cross-validation procedure identified seven inconsistent fossils. Using the consistent fossils as calibration points, we found that despite their overall antiquity as a lineage, the most species-rich clades of turtles diversified well within the Cenozoic. Many of the truly ancient lineages of turtles are currently represented by a few, often endangered species that deserve high priority as conservation targets.

Animals↗

Mammalian genes as molecular clocks?

In analyzing the silent nucleotide substitutions in some mammalian mitochondrial mRNA coding genes, we had found that the frequency of each of the four nucleotides in rat, mouse, and cow, but not in humans, is the same in the silent third codon position (Lanave C, Preparata G, Saccone C, Serio G (1984) J Mol Evol 20:86-93). Because our findings for these three species were compatible with a stationary Markov process for the evolution of nucleotide sequences, we applied such a model to calculate the effective evolutionary silent substitution rate (vs) and the divergence times among the species. In this paper we have analyzed the first and second codon positions in the same mammalian mitochondrial genes. We found that in the first and second codon positions the human mitochondrial genes satisfy the stationarity conditions. This has allowed us to use the stochastic model mentioned above to calculate the divergence times among mouse, rat, cow, and human. Furthermore, we have analyzed the silent substitution rate in one nuclear gene for these four mammals. We found that in this gene the effective silent substitution rate is about 3 times lower than in mitochondrial genes, and that humans are in this case stationary with respect to the other three mammals in the third codon position as well. Application of our Markov model to this latter gene yields divergence times consistent with our previous determinations.

Animals↗

Timing the eastern Asian-eastern North American floristic disjunction: molecular clock corroborates paleontological estimates.

Sequence data of the chloroplast gene rbcL were used to estimate the time of the well-known eastern Asian-eastern North American floristic disjunction. Sequence divergence of rbcL was examined for 22 species of 11 genera (Campsis, Caulophyllum, Cornus, Decumaria, Liriodendron, Menispermum, Mitchella, Pachysandra, Penthorum, Podophyllum, and Phryma) representing a diverse array of flowering plants occurring disjunctly in eastern Asia and eastern North America. Divergence times of putative disjunct species pairs were estimated from synonymous substitutions, using rbcL molecular clocks calibrated for Cornus. Relative rate tests were performed to assess rate constancy of rbcL evolution among lineages. Corrections of estimates of divergence times for each species pair were made based on rate differences of rbcL between Cornus and other species pairs. Results of these analyses indicate that the time of divergence of species pairs examined ranges from 12.56 +/- 4.30 million years to recent (<0.31 million years), with most within the last 10 million years (in the late Miocene and Pliocene). These results suggest that the isolation of most morphologically similar disjunct species in eastern Asia and eastern North America occurred during the global climatic cooling period that took place throughout the late Tertiary and Quaternary. This estimate is closely correlated with paleontological evidence and in agreement with the hypothesis that considers the eastern Asian-eastern North American floristic disjunction to be the result of the range restriction of a once more or less continuously distributed mixed mesophytic forest of the Northern Hemisphere that occurred during the late Tertiary and Quaternary. This implies that in most taxa the disjunction may have resulted from vicariance events. However, long-distance dispersal may explain the disjunct distribution of taxa with low divergence, such as Menispermum.

Asia↗

Phylogenetics of Perissodactyla and tests of the molecular clock.

Two mitochondrial genes, the protein-coding cytochrome c oxidase subunit II (COII) gene and a portion of the 12S rRNA gene, were used for phylogenetic investigation of the mammalian order Perissodactyla. The primary objective of the study was to utilize the extensive fossil record of perissodactyls for calibrating molecular clocks and comparing estimates of divergence times using both genes and two fossil calibration points. Secondary objectives included clarification of previously unresolved relationships within Tapiridae and comparison of the results of separate and combined analyses of two genes. Analyses included several perissodactyl lineages representing all three families (Tapiridae, Equidae, and Rhinocerotidae), most extant genera, all four species of tapirs, two to four species of rhinoceros, and two species of Equus. The application of a relatively recent fossil calibration point and a relatively ancient calibration point produced greatly different estimates of evolutionary rates and divergence times for both genes, even though a relative rates test did not find significant rate differences among taxa. A likelihood-ratio test, however, rejected a molecular clock for both genes. Neither calibration point produced estimates of divergence times consistent with paleontological evidence over a range of perissodactyl radiations. The combined analysis of both genes produces a well-resolved phylogeny with Perissodactyla that conforms to traditional views of interfamilial relationships and supports monophyly of neotropical tapirs. Combining the data sets increases support for most nodes but decreases the support for a neotropical tapir clade because the COII and 12S rRNA data sets are in conflict for tapir relationships.

Animals↗

Polymorphism and divergence for island-model species.

Estimates of the scaled selection coefficient, gamma of Sawyer and Hartl, are shown to be remarkably robust to population subdivision. Estimates of mutation parameters and divergence times, in contrast, are very sensitive to subdivision. These results follow from an analysis of natural selection and genetic drift in the island model of subdivision in the limit of a very large number of subpopulations, or demes. In particular, a diffusion process is shown to hold for the average allele frequency among demes in which the level of subdivision sets the timescale of drift and selection and determines the dynamic equilibrium of allele frequencies among demes. This provides a framework for inference about mutation, selection, divergence, and migration when data are available from a number of unlinked nucleotide sites. The effects of subdivision on parameter estimates depend on the distribution of samples among demes. If samples are taken singly from different demes, the only effect of subdivision is in the rescaling of mutation and divergence-time parameters. If multiple samples are taken from one or more demes, high levels of within-deme relatedness lead to low levels of intraspecies polymorphism and increase the number of fixed differences between samples from two species. If subdivision is ignored, mutation parameters are underestimated and the species divergence time is overestimated, sometimes quite drastically. Estimates of the strength of selection are much less strongly affected and always in a conservative direction.

Biological Evolution↗

Codon equilibrium II: Its use in estimating silent-substitution rates.

We study the equilibrium in the use of synonymous codons by eukaryotic organisms and find five equations involving substitution rates that we believe embody the important implications of equilibrium for the process of silent substitution. We then combine these five equations with additional criteria to determine sets of substitution rates applicable to eukaryotic organisms. One method employs the equilibrium equations and a principle of maximum entropy to find the most uniform set of rates consistent with equilibrium. In a second method we combine the equilibrium equations with data on the man-mouse divergence to determine that set of rates that is most neutral yet consistent with both types of data (i.e., equilibrium and divergence data). Simulations show this second method to be quite reliable in spite of significant saturation in the substitution process. We find that when divergence data are included in the calculation of rates, even though these rates are chosen to be as neutral as possible, the strength of selection inferred from the nonuniformity of the rates is approximately doubled. Both sets of rates are applied to estimate the human-mouse divergence time based on several independent subsets of the divergence data consisting of the quartet, C- or T-ending duet, and A- or G-ending duet codon sets. Both rate sets produce patterns of divergence times that are shortest for the quartet data, intermediate for the CT-ending duets, and longest for the AG-ending duets. This indicates that rates of transitions in the duet-codon sets are significantly higher than those in the quartet-codon sets; this effect is especially marked for A----G, the rate of which in duets must be about double that in quartets.

Animals↗

Estimating effective population size from samples of sequences: inefficiency of pairwise and segregating sites as compared to phylogenetic estimates.

It is known that under neutral mutation at a known mutation rate a sample of nucleotide sequences, within which there is assumed to be no recombination, allows estimation of the effective size of an isolated population. This paper investigates the case of very long sequences, where each pair of sequences allows a precise estimate of the divergence time of those two gene copies. The average divergence time of all pairs of copies estimates twice the effective population number and an estimate can also be derived from the number of segregating sites. One can alternatively estimate the genealogy of the copies. This paper shows how a maximum likelihood estimate of the effective population number can be derived from such a genealogical tree. The pairwise and the segregating sites estimates are shown to be much less efficient than this maximum likelihood estimate, and this is verified by computer simulation. The result implies that there is much to gain by explicitly taking the tree structure of these genealogies into account.

Likelihood Functions↗

Cytonuclear disequilibria in wild populations of rabbit (Oryctolagus cuniculus L.) suggest unequal allele turnover rates at the b locus (IGKC1).

DNA sequence comparisons suggest that evolutionary rates at the rabbit IGKC1 locus can differ among allelic lineages. Here we address the question of whether population turnover rates can vary among IGKC1 alleles. We studied the distribution of sixteen IGKC1 (or b-locus) allotypes in areas comprising the aboriginal species range (Iberian peninsula). Rabbits in this area belong to one of two distantly related mitochondrial lineages (mtDNA types) A and B. In the more recent distribution area of the species, all rabbits belong to the mtDNA type B lineage, and IGKC1 alleles b4 and b5 comprise over 90% of the gene pool. These two alleles are also predominant in areas of mtDNA type B prevalence within the Iberian range. However, in areas of mtDNA type A prevalence, the b4 and b5 allotypes are rare or absent; they apparently have been replaced by serologically related, but distinct, 'endemic' variants. The cytonuclear disequilibria were highly significant, also within the subsample consisting of populations from Spain. These observations suggest that allelic persistence times for the predominant IGKC1 lineages could be shorter than the divergence time of the major mtDNA lineages A and B. In contrast, the relative gene frequencies of the IGKC1 allele b9 were similar among the type A and type B rabbits; it was present in most populations at low frequency. In consequence, persistence times of the b9 allele appear to be longer than the divergence time of lineages A and B. The data reported here are in agreement with the DNA sequence data, providing further proof that the molecular clock can run at different rates among allelic lineages at the rabbit IGKC1 locus.

Alleles↗

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species.

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein-coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).

Amino Acid Sequence↗

Single-copy DNA distance between two congeneric sea urchin species exhibiting radically different modes of development.

We have investigated the differences between nuclear genomes of two purportedly congeneric species of sea urchin that differ radically in early development. Heliocidaris tuberculata develops by means of a typical pluteus larva, whereas H. erythrogramma develops directly from an egg that is 100-fold the volume of the H. tuberculata egg. Reassociation kinetic analysis shows that the kinetic components of the genomic DNA from the two species are essentially the same. No single repeat component explains the 30% difference between the H. erythrogramma and H. tuberculata genomes. Reciprocal hybridization of tracer-labeled single-copy DNA fractions between these species indicates that approximately 50% of the single-copy DNA is sufficiently similar to form hybrids at standard hybridization criterion. Thermal denaturation profiles of the hybridized single-copy DNA sequence yields median (T50H) values of 13.8 degrees-16.5 degrees C. This result suggests a divergence time of 10-13 Mya, which is comparable to divergence times between congeneric sea urchin species in other genera that do not differ significantly in development. Radical differences in early developmental processes can evolve rapidly between closely related forms.

Animals↗