On the molecular evolutionary clock

Phylo-Epigenetics in Phylogeny Analyses and Evolution

Long-standing, continuous blurring and controversies in the field of phylogenetic interspecies relations, associated with insufficient explanations for dynamics and variability of speeds of evolution in mammals, hint at a crucial missing link. It has been suggested that transgenerational epigenetic inheritance and the concealed mechanisms behind play a distinct role in mammalian evolution. Here, a comprehensive sequence alignment approach in hominid species, i.e., Homo sapiens, Homo neanderthalensis, Denisovan human, Pan troglodytes, Pan paniscus, Gorilla gorilla, and Pongo pygmaeus, comprising conserved CpG islands of housekeeping genes, uncover evidence for a distinct variability of CpG dinucleotides. Applying solely these evolutionary consistent and inconsistent CpG sites in a classic phylogenetic analysis, calibrated by the divergence time point of the common chimpanzee (P. troglodytes) and the bonobo or pygmy chimpanzee (P. paniscus), a "phylo-epigenetic" tree has been generated, which precisely recapitulates branch points and branch lengths, i.e., divergence events and relations, as they have been broadly suggested in the current literature, based on comprehensive molecular phylogenomics and fossil records of many decades. It is suggested here that CpG dinucleotide changes at CpG islands are of superior importance for evolutionary developments. These changes are successfully inherited through the germ line, determining emerging methylation profiles, and they are a central component of transgenerational epigenetic inheritance. It is hidden in the DNA, what will happen on it later.

Confidence intervals in molecular dating by maximum likelihood

Molecular dating has been widely used to infer the times of past evolutionary events using molecular sequences. This paper describes three bootstrap methods to infer confidence intervals under a penalized likelihood framework. The basic idea is to use data pseudoreplicates to infer uncertainty in the branch lengths of a phylogeny reconstructed with molecular sequences. The three specific bootstrap methods are nonparametric (direct tree bootstrapping), semiparametric (rate smoothing), and parametric (Poisson simulation). Our extensive simulation study showed that the three methods perform generally well under a simple strict clock model of molecular evolution; however, the results were less positive with data simulated using an uncorrelated or a correlated relaxed clock model. Several factors impacted, possibly in interaction, the performance of the confidence intervals. Increasing the number of calibration points had a positive effect, as well as increasing the sequence length or the number of sequences although both latter effects depended on the model of evolution. A case study is presented with a molecular phylogeny of the Felidae (Mammalia: Carnivora). A comparison was made with a Bayesian analysis: the results were very close in terms of confidence intervals and there was no marked tendency for an approach to produce younger or older bounds compared to the other.

The performance of outgroup-free rooting under evolutionary radiations

Tree rooting implies a temporal dimension to phylogenies. Only after defining the position of the root node is that the ancestral-descendant relationship between branches can be fully deduced. Rooting has been usually carried out by employing evolutionarily close outgroup lineages, which is a drawback when these lineages are unavailable or unknown. Alternatively, outgroup-free rooting methods were proposed, which rely on the constancy of evolutionary rates to varying degrees. In this work we analyzed the performance of two of these methods, the midpoint rooting (MPR) and the minimal ancestor deviation (MAD), in rooting topologies evolved under challenging scenarios of fast evolutionary radiations derived from empirical data, characterized by short internal branches near the crown node. Considering all branch length combinations investigated, both methods exhibited average success rates below 50%, although MAD slightly outperformed MPR. Moreover, tree balance significantly impacted the relative performance of the methods. We found that, in four-taxa unrooted trees, the outcome of whether both methodologies will correctly root the tree can be roughly predicted by two simple dimensionless metrics: the coefficient of variation of the external branch lengths, and the ratio between the internal branch length to the total sum of branch lengths, which were employed to devise a general linear model that allowed calculating the probability of correct placing the root node for any four-taxa tree. We predicted that the performance of both outgroup-free rooting methods on loci representing the placental mammal radiation ranged between 50% and 75%.

Spectrum of Protein Location in Proteomes Captures Evolutionary Relationship Between Species

The native subcellular location (also referred to as localization or cellular compartment) of a protein is the one in which it acts most frequently; it is one aspect of protein function. Do ten eukaryotic model organisms differ in their location spectrum, i.e., the fraction of its proteome in each of seven major cellular compartments? As experimental annotations of locations remain biased and incomplete, we need prediction methods to answer this question. After systematic bias corrections, the complete but faulty prediction methods appeared to be more appropriate to compare location spectra between species than the incomplete more accurate experimental data. This work compared the location spectra for ten eukaryotes: Homo sapiens (human), Gorilla gorilla (gorilla), Pan troglodytes (chimpanzee), Mus musculus (mouse), Rattus norvegicus (rat), Drosophila melanogaster (fruit/vinegar fly), Anopheles gambiae (African malaria mosquito), Caenorhabitis elegans (nematode), Saccharomyces cerevisiae (baker’s yeast), and Schizosaccharomyces pombe (fission yeast). The two largest classes were predicted to be the nucleus and the cytoplasm together accounting for 47–62% of all proteins, while 7–21% of the proteins were predicted in the plasma membrane and 4–15% to be secreted. Overall, the predicted location spectra were largely similar. However, in detail, the differences sufficed to plot trees (UPGMA) and 2D (PCA) maps relating the ten organisms using a simple Euclidean distance in seven states (location classes). The relations based on the simple predicted location spectra captured aspects of cross-species comparisons usually revealed only by much more detailed evolutionary comparisons. Most interestingly, known phylogenetic relations were reproduced better by paralog-only than by ortholog-only trees.

Probing periodontal microbial dark matter using metataxonomics and metagenomics

Our view of the periodontal microbial community has been shaped by a century or more of cultivation-based and microscopic investigations. While these studies firmly established the infection-mediated etiology of periodontal diseases, it was apparent from the very early days that periodontal microbiology suffered from what Staley and Konopka described as the "great plate count anomaly", in that these culturable bacteria were only a minor part of what was visible under the microscope. For nearly a century, much effort has been devoted to finding the right tools to investigate this uncultivated majority, also known as "microbial dark matter". The discovery that DNA was an effective tool to "see" microbial dark matter was a significant breakthrough in environmental microbiology, and oral microbiologists were among the earliest to capitalize on these advances. By identifying the order in which nucleotides are arranged in a stretch of DNA (DNA sequencing) and creating a repository of these sequences, sequence databases were created. Computational tools that used probability-driven analysis of these sequences enabled the discovery of new and unsuspected species and ascribed novel functions to these species. This review will trace the development of DNA sequencing as a quantitative, open-ended, comprehensive approach to characterize microbial communities in their native environments, and explore how this technology has shifted traditional dogmas on how the oral microbiome promotes health and its role in disease causation and perpetuation.

Epigenetic pacemaker: closed form algebraic solutions

Background

DNA methylation is widely used as a biomarker in crucial medical applications as well as for human age prediction of very high accuracy. This biomarker is based on the methylation status of several hundred CpG sites. In a recent line of publications we have adapted a versatile concept from evolutionary biology - the Universal Pacemaker (UPM) - to the setting of epigenetic aging and denoted it the Epigenetic PaceMaker (EPM). The EPM, as opposed to other epigenetic clocks, is not confined to specific pattern of aging, and the epigenetic age of the individual is inferred independently of other individuals. This allows an explicit modeling of aging trends, in particular non linear relationship between chronological and epigenetic age. In one of these recent works, we have presented an algorithmic improvement based on a two-step conditional expectation maximization (CEM) algorithm to arrive at a critical point on the likelihood surface. The algorithm alternates between a time step and a site step while advancing on the likelihood surface.

Results

Here we introduce non trivial improvements to these steps that are essential for analyzing data sets of realistic magnitude in a manageable time and space. These structural improvements are based on insights from linear algebra and symbolic algebra tools, providing us greater understanding of the degeneracy of the complex problem space. This understanding in turn, leads to the complete elimination of the bottleneck of cumbersome matrix multiplication and inversion, yielding a fast closed form solution in both steps of the CEM.In the experimental results part, we compare the CEM algorithm over several data sets and demonstrate the speedup obtained by the closed form solutions. Our results support the theoretical analysis of this improvement.

Conclusions

These improvements enable us to increase substantially the scale of inputs analyzed by the method, allowing us to apply the new approach to data sets that could not be analyzed before.

An epigenetic pacemaker is detected via a fast conditional expectation maximization algorithm

AIM

DNA methylation has proven to be a remarkably accurate biomarker for human age, allowing the prediction of chronological age to within a couple of years. Recently, we proposed that the Universal PaceMaker (UPM), a flexible paradigm for modeling evolution, could be applied to epigenetic aging. Nevertheless, application to real data was restricted to small datasets for technical limitations.

MATERIALS & METHODS

We partition the set of variables into to two subsets and optimize the likelihood function on each set separately. This yields an extremely efficient Conditional Expectation Maximization algorithm, alternating between the two sets while increasing the overall likelihood.

RESULTS

Using the technique, we could reanalyze datasets of larger magnitude and show significant advantage to the UPM approach.

CONCLUSION

The UPM more faithfully models epigenetic aging than the time linear approach while methylated sites accelerate and decelerate jointly.

Neutral Theory and the Somatic Evolution of Cancer

Kimura's neutral theory argued that positive selection was not responsible for an appreciable fraction of molecular substitutions. Correspondingly, quantitative analysis reveals that the vast majority of substitutions in cancer genomes are not detectably under selection. Insights from the somatic evolution of cancer reveal that beneficial substitutions in cancer constitute a small but important fraction of the molecular variants. The molecular evolution of cancer community will benefit by incorporating the neutral theory of molecular evolution into their understanding and analysis of cancer evolution-and accepting the use of tractable, predictive models, even when there is some evidence that they are not perfect.

DNA DIVERGENCE AMONG HOMINOIDS

We have determined the degree of single-copy DNA divergence among the extant members of the Hominoidea employing the technique of DNA-DNA hybridization. The species studied include humans, two species of chimpanzees, gorillas, two subspecies of orangutans, and two species of gibbons; as an outgroup we have used a member of the Old World monkeys (Cercopithecidae), the baboon. Our methods are different from those previously used and allow us to control for two factors other than base-pair mismatch that can affect the thermal stability of DNA duplexes: the base composition and duplex length. In addition, we have studied more than one individual for most species and thus are able to assess the effect of intraspecific variation on phylogenetic conclusions. The results indicate that the closest extant relatives of humans are the chimpanzees. Gorillas are the next closest, followed by orangutans and gibbons. This result is strongly supported statistically, as there is virtually no overlap in measurements between different taxa. Our conclusions are in agreement with a growing amount of molecular evidence supporting this pattern of relatedness. The data behave as a reasonably good molecular clock, and we do not see an indication of slowdown in molecular evolution in the clade containing humans and African apes, contrary to what has been documented for protein-coding regions. Because of the clocklike nature of the results, we have estimated that the divergence of humans and chimpanzees occurred about 6-8 million years ago. Results from orangutans indicate that the Borneo and Sumatra populations are genetically distinct, about as different as the named species of chimpanzees.

A Statistical Framework to Identify Deviation from Time Linearity in Epigenetic Aging

In multiple studies DNA methylation has proven to be an accurate biomarker of age. To develop these biomarkers, the methylation of multiple CpG sites is typically linearly combined to predict chronological age. By contrast, in this study we apply the Universal PaceMaker (UPM) model to investigate changes in DNA methylation during aging. The UPM was initially developed to study rate acceleration/deceleration in sequence evolution. Rather than identifying which linear combinations of sites predicts age, the UPM models the rates of change of multiple CpG sites, as well as their starting methylation levels, and estimates the age of each individual to optimize the model fit. We refer to the estimated age as the “epigenetic age”, which is in contrast to the known chronological age of each individual. We construct a statistical framework and devise an algorithm to determine whether a genomic pacemaker is in effect (i.e rates of change vary with age). The decision is made by comparing two competing likelihood based models, the molecular clock (MC) and UPM. For the molecular clock model, we use the known chronological age of each individual and fit the methylation rates at multiple sites, and express the problem as a linear least squares and solve it in polynomial time. For the UPM case, the search space is larger as we are fitting both the epigenetic age of each individual as well as the rates for each site, yet we succeed to reduce the problem to the space of individuals and polynomial in the more significant space—the methylated sites. We first tested our algorithm on simulated data to elucidate the factors affecting the identification of the pacemaker model. We find that, provided with enough data, our algorithm is capable of identifying a pacemaker even when a weak signal is present in the data. Based on these results, we applied our method to DNA methylation data from human blood from individuals of various ages. Although the improvement in variance across sites between the UPM and MC was small, the results suggest that the existence of a pacemaker is highly significant. The PaceMaker results also suggest a decay in the rate of change in DNA methylation with age.

DNA methylation is an important component of the epigenetic code that defines and maintains the state of cells. Recently, it has been found that certain sites in the genome undergo methylation changes at different rates during aging. The seminal work of Steve Horvath found that the methylation of a couple hundred CpG sites could be linearly combined to accurately predict the age of an individual in a number of tissues. Such a pattern resembles the Molecular Clock (MC) concept prevailing in molecular evolution, which suggests that there are sites in the genome that change linearly with age. In this work, we adapt the Universal PaceMaker (UPM) model to the setting of DNA methylation changes during aging. UPM relaxes the rate constancy of MC and was found to provide a better statistical explanation for genome evolution across the entire tree of life. This adaptation requires the solution of a complex optimization problem. Nevertheless, in a series of observations we show that the problem can be solved efficiently under the MC model and slightly less efficiently under the UPM model. This allows us to solve problems of non-trivial size. We chose as a proof of concept to analyze DNA methylation data collected from the blood of humans of different ages. Our results show that, similarly to genome evolution, the UPM provided an improvement of about 2% in the fit to the data. The statistical significance of this improvement is very high. Although tested on a small data set, this improvement demonstrates that the UPM more accurately captures age related DNA methylation changes than the MC model.

Characterization of Complete Genome Sequences of a Porcine Endogenous Retrovirus Isolated From China Bama Minipig Reveals an Evolutionary Time Earlier Than That of Isolates From European Minipigs

BACKGROUND

A porcine endogenous retroviruses (PERV) isolate, PERV-A-BM, was isolated from a Guangxi Bama minipig in China.

METHODS

To understand its genetic variation and evolution, the complete PERV-A-BM genome sequences were determined and compared with isolates from different Sus scrofa breeds and porcine cell lines. A total of 69 nucleotide substitutions were found in the full-length genome, including 26 non-synonymous mutations.

RESULTS

Phylogenetic trees based on the complete genome sequence as well as the gag, pol, and env gene sequences from 21 PERV isolates demonstrated that the PERV-A-BM was closely related to the EF133960 isolate from Chinese Wuzhishan miniature pigs inbred in Hainan, China, and distantly related to strains isolated from European-born pigs.

CONCLUSIONS

The estimation of age in the proviral PERV-A-BM integrating into the host genome reveals that the age of PERV-A-BM is at least 8.3 × 10(6) years, an evolutionary time earlier than that of isolates from European-born pigs.

Histories of molecules: Reconciling the past

Molecular data and methods have become centrally important to evolutionary analysis, largely because they have enabled global phylogenetic reconstructions of the relationships between organisms in the tree of life. Often, however, molecular stories conflict dramatically with morphology-based histories of lineages. The evolutionary origin of animal groups provides one such case. In other instances, different molecular analyses have so far proved irreconcilable. The ancient and major divergence of eukaryotes from prokaryotic ancestors is an example of this sort of problem. Efforts to overcome these conflicts highlight the role models play in phylogenetic reconstruction. One crucial model is the molecular clock; another is that of 'simple-to-complex' modification. I will examine animal and eukaryote evolution against a backdrop of increasing methodological sophistication in molecular phylogeny, and conclude with some reflections on the nature of historical science in the molecular era of phylogeny.

On the number of genomic pacemakers: a geometric approach

The universal pacemaker (UPM) model extends the classical molecular clock (MC) model, by allowing each gene, in addition to its individual intrinsic rate as in the MC, to accelerate or decelerate according to the universal pacemaker. Under UPM, the relative evolutionary rates of all genes remain nearly constant whereas the absolute rates can change arbitrarily. It was shown on several taxa groups spanning the entire tree of life that the UPM model describes the evolutionary process better than the MC model. In this work we provide a natural generalization to the UPM model that we denote multiple pacemakers (MPM). Under the MPM model every gene is still affected by a single pacemaker, however the number of pacemakers is not confined to one. Such a model induces a partition over the gene set where all the genes in one part are affected by the same pacemaker and task is to identify the pacemaker partition, or in other words, finding for each gene its associated pacemaker. We devise a novel heuristic procedure, relying on statistical and geometrical tools, to solve the problem and demonstrate by simulation that this approach can cope satisfactorily with considerable noise and realistic problem sizes. We applied this procedure to a set of over 2000 genes in 100 prokaryotes and demonstrated the significant existence of two pacemakers.

An evaluation of fossil tip-dating versus node-age calibrations in tetraodontiform fishes (Teleostei: Percomorphaceae)

Time-calibrated phylogenies based on molecular data provide a framework for comparative studies. Calibration methods to combine fossil information with molecular phylogenies are, however, under active development, often generating disagreement about the best way to incorporate paleontological data into these analyses. This study provides an empirical comparison of the most widely used approach based on node-dating priors for relaxed clocks implemented in the programs BEAST and MrBayes, with two recently proposed improvements: one using a new fossilized birth-death process model for node dating (implemented in the program DPPDiv), and the other using a total-evidence or tip-dating method (implemented in MrBayes and BEAST). These methods are applied herein to tetraodontiform fishes, a diverse group of living and extinct taxa that features one of the most extensive fossil records among teleosts. Previous estimates of time-calibrated phylogenies of tetraodontiforms using node-dating methods reported disparate estimates for their age of origin, ranging from the late Jurassic to the early Paleocene (ca. 150-59Ma). We analyzed a comprehensive dataset with 16 loci and 210 morphological characters, including 131 taxa (95 extant and 36 fossil species) representing all families of fossil and extant tetraodontiforms, under different molecular clock calibration approaches. Results from node-dating methods produced consistently younger ages than the tip-dating approaches. The older ages inferred by tip dating imply an unlikely early-late Jurassic (ca. 185-119Ma) origin for this order and the existence of extended ghost lineages in their fossil record. Node-based methods, by contrast, produce time estimates that are more consistent with the stratigraphic record, suggesting a late Cretaceous (ca. 86-96Ma) origin. We show that the precision of clade age estimates using tip dating increases with the number of fossils analyzed and with the proximity of fossil taxa to the node under assessment. This study suggests that current implementations of tip dating may overestimate ages of divergence in calibrated phylogenies. It also provides a comprehensive phylogenetic framework for tetraodontiform systematics and future comparative studies.

The genome of the chicken DT40 bursal lymphoma cell line

The chicken DT40 cell line is a widely used model system in the study of multiple cellular processes due to the efficiency of homologous gene targeting. The cell line was derived from a bursal lymphoma induced by avian leukosis virus infection. In this study we characterized the genome of the cell line using whole genome shotgun sequencing and single nucleotide polymorphism array hybridization. The results indicate that wild-type DT40 has a relatively normal karyotype, except for whole chromosome copy number gains, and no karyotype variability within stocks. In a comparison to two domestic chicken genomes and the Gallus gallus reference genome, we found no unique mutational processes shaping the DT40 genome except for a mild increase in insertion and deletion events, particularly deletions at tandem repeats. We mapped coding sequence mutations that are unique to the DT40 genome; mutations inactivating the PIK3R1 and ATRX genes likely contributed to the oncogenic transformation. In addition to a known avian leukosis virus integration in the MYC gene, we detected further integration sites that are likely to de-regulate gene expression. The new findings support the hypothesis that DT40 is a typical transformed cell line with a relatively intact genome; therefore, it is well-suited to the role of a model system for DNA repair and related processes. The sequence data generated by this study, including a searchable de novo genome assembly and annotated lists of mutated genes, will support future research using this cell line.

Universal pacemaker of genome evolution in animals and fungi and variation of evolutionary rates in diverse organisms

Gene evolution is traditionally considered within the framework of the molecular clock (MC) model whereby each gene is characterized by an approximately constant rate of evolution. Recent comparative analysis of numerous phylogenies of prokaryotic genes has shown that a different model of evolution, denoted the Universal PaceMaker (UPM), which postulates conservation of relative, rather than absolute evolutionary rates, yields a better fit to the phylogenetic data. Here, we show that the UPM model is a better fit than the MC for genome wide sets of phylogenetic trees from six species of Drosophila and nine species of yeast, with extremely high statistical significance. Unlike the prokaryotic phylogenies that include distant organisms and multiple horizontal gene transfers, these are simple data sets that cover groups of closely related organisms and consist of gene trees with the same topology as the species tree. The results indicate that both lineage-specific and gene-specific rates are important in genome evolution but the lineage-specific contribution is greater. Similar to the MC, the gene evolution rates under the UPM are strongly overdispersed, approximately 2-fold compared with the expectation from sampling error alone. However, we show that neither Drosophila nor yeast genes form distinct clusters in the tree space. Thus, the gene-specific deviations from the UPM, although substantial, are uncorrelated and most likely depend on selective factors that are largely unique to individual genes. Thus, the UPM appears to be a key feature of genome evolution across the history of cellular life.

Molecular phylogenetics before sequences: oligonucleotide catalogs as k-mer spectra

From 1971 to 1985, Carl Woese and colleagues generated oligonucleotide catalogs of 16S/18S rRNAs from more than 400 organisms. Using these incomplete and imperfect data, Carl and his colleagues developed unprecedented insights into the structure, function, and evolution of the large RNA components of the translational apparatus. They recognized a third domain of life, revealed the phylogenetic backbone of bacteria (and its limitations), delineated taxa, and explored the tempo and mode of microbial evolution. For these discoveries to have stood the test of time, oligonucleotide catalogs must carry significant phylogenetic signal; they thus bear re-examination in view of the current interest in alignment-free phylogenetics based on k-mers. Here we consider the aims, successes, and limitations of this early phase of molecular phylogenetics. We computationally generate oligonucleotide sets (e-catalogs) from 16S/18S rRNA sequences, calculate pairwise distances between them based on D₂ statistics, compute distance trees, and compare their performance against alignment-based and k-mer trees. Although the catalogs themselves were superseded by full-length sequences, this stage in the development of computational molecular biology remains instructive for us today.

Stability along with extreme variability in core genome evolution

The shape of the distribution of evolutionary distances between orthologous genes in pairs of closely related genomes is universal throughout the entire range of cellular life forms. The near invariance of this distribution across billions of years of evolution can be accounted for by the Universal Pace Maker (UPM) model of genome evolution that yields a significantly better fit to the phylogenetic data than the Molecular Clock (MC) model. Unlike the MC, the UPM model does not assume constant gene-specific evolutionary rates but rather postulates that, in each evolving lineage, the evolutionary rates of all genes change (approximately) in unison although the pacemakers of different lineages are not necessarily synchronized. Here, we dissect the nearly constant evolutionary rate distribution by comparing the genome-wide relative rates of evolution of individual genes in pairs or triplets of closely related genomes from diverse bacterial and archaeal taxa. We show that, although the gene-specific relative rate is an important feature of genome evolution that explains more than half of the variance of the evolutionary distances, the ranges of relative rate variability are extremely broad even for universal genes. Because of this high variance, the gene-specific rate is a poor predictor of the conservation rank for any gene in any particular lineage.

Universal pacemaker of genome evolution

A fundamental observation of comparative genomics is that the distribution of evolution rates across the complete sets of orthologous genes in pairs of related genomes remains virtually unchanged throughout the evolution of life, from bacteria to mammals. The most straightforward explanation for the conservation of this distribution appears to be that the relative evolution rates of all genes remain nearly constant, or in other words, that evolutionary rates of different genes are strongly correlated within each evolving genome. This correlation could be explained by a model that we denoted Universal PaceMaker (UPM) of genome evolution. The UPM model posits that the rate of evolution changes synchronously across genome-wide sets of genes in all evolving lineages. Alternatively, however, the correlation between the evolutionary rates of genes could be a simple consequence of molecular clock (MC). We sought to differentiate between the MC and UPM models by fitting thousands of phylogenetic trees for bacterial and archaeal genes to supertrees that reflect the dominant trend of vertical descent in the evolution of archaea and bacteria and that were constrained according to the two models. The goodness of fit for the UPM model was better than the fit for the MC model, with overwhelming statistical significance, although similarly to the MC, the UPM is strongly overdispersed. Thus, the results of this analysis reveal a universal, genome-wide pacemaker of evolution that could have been in operation throughout the history of life.

Distinct co-evolution patterns of genes associated to DNA polymerase III DnaE and PolC

BACKGROUND

Bacterial genomes displaying a strong bias between the leading and the lagging strand of DNA replication encode two DNA polymerases III, DnaE and PolC, rather than a single one. Replication is a highly unsymmetrical process, and the presence of two polymerases is therefore not unexpected. Using comparative genomics, we explored whether other processes have evolved in parallel with each polymerase.

RESULTS

Extending previous in silico heuristics for the analysis of gene co-evolution, we analyzed the function of genes clustering with dnaE and polC. Clusters were highly informative. DnaE co-evolves with the ribosome, the transcription machinery, the core of intermediary metabolism enzymes. It is also connected to the energy-saving enzyme necessary for RNA degradation, polynucleotide phosphorylase. Most of the proteins of this co-evolving set belong to the persistent set in bacterial proteomes, that is fairly ubiquitously distributed. In contrast, PolC co-evolves with RNA degradation enzymes that are present only in the A+T-rich Firmicutes clade, suggesting at least two origins for the degradosome.

CONCLUSION

DNA replication involves two machineries, DnaE and PolC. DnaE co-evolves with the core functions of bacterial life. In contrast PolC co-evolves with a set of RNA degradation enzymes that does not derive from the degradosome identified in gamma-Proteobacteria. This suggests that at least two independent RNA degradation pathways existed in the progenote community at the end of the RNA genome world.

Waiting for sequences: Morris Goodman, immunodiffusion experiments, and the origins of molecular anthropology

During the early 1960s, Morris Goodman used a variety of immunological tests to demonstrate the very close genetic relationships among humans, chimpanzees, and gorillas. Molecular anthropologists often point to this early research as a critical step in establishing their new specialty. Based on his molecular results, Goodman challenged the widely accepted taxonomic classification that separated humans from chimpanzees and gorillas in two separate families. His claim that chimpanzees and gorillas should join humans in family Hominidae sparked a well-known conflict with George Gaylord Simpson, Ernst Mayr, and other prominent evolutionary biologists. Less well known, but equally significant, were a series of disagreements between Goodman and other prominent molecular evolutionists concerning both methodological and theoretical issues. These included qualitative versus quantitative data, the role of natural selection, rates of evolution, and the reality of molecular clocks. These controversies continued throughout Goodman's career, even as he moved from immunological techniques to protein and DNA sequence analysis. This episode highlights the diversity of methods used by molecular evolutionists and the conflicting conclusions drawn from the data that these methods generated.

History in the gene: negotiations between molecular and organismal anthropology

In the advertising discourse of human genetic database projects, of genetic ancestry tracing companies, and in popular books on anthropological genetics, what I refer to as the anthropological gene and genome appear as documents of human history, by far surpassing the written record and oral history in scope and accuracy as archives of our past. How did macromolecules become "documents of human evolutionary history"? Historically, molecular anthropology, a term introduced by Emile Zuckerkandl in 1962 to characterize the study of primate phylogeny and human evolution on the molecular level, asserted its claim to the privilege of interpretation regarding hominoid, hominid, and human phylogeny and evolution vis-à-vis other historical sciences such as evolutionary biology, physical anthropology, and paleoanthropology. This process will be discussed on the basis of three key conferences on primate classification and evolution that brought together exponents of the respective fields and that were held in approximately ten-years intervals between the early 1960s and the 1980s. I show how the anthropological gene and genome gained their status as the most fundamental, clean, and direct records of historical information, and how the prioritizing of these epistemic objects was part of a complex involving the objectivity of numbers, logic, and mathematics, the objectivity of machines and instruments, and the objectivity seen to reside in the epistemic objects themselves.

Studying the effect of environmental change on biotic evolution: past genetic contributions, current work and future directions

Evolutionary geneticists currently face a major scientific opportunity when integrating across the rapidly increasing amount of genetic data and existing biological scenarios based on ecology, fossils or climate models. Although genetic data acquisition and analysis have improved tremendously, several limitations remain. Here, we discuss the feedback between history and genetic variation in the face of environmental change with increasing taxonomic and temporal scale, as well as the major challenges that lie ahead. In particular, we focus on recent developments in two promising genetic methods, those of 'phylochronology' and 'molecular clocks'. With the advent of ancient DNA techniques, we can now directly sample the recent past. We illustrate this amazing and largely untapped utility of ancient DNA extracted from accurately dated localities with documented environmental changes. Innovative statistical analyses of these genetic data expose the direct effect of recent environmental change on genetic endurance, or maintenance of genetic variation. The 'molecular clock' (assumption of a linear relationship between genetic distance and evolutionary time) has been used extensively in phylogenetic studies to infer time and correlation between lineage divergence time and concurrent environmental change. Several studies at both population and species scale support a persuasive relationship between particular perturbation events and time of biotic divergence. However, we are still a way from gleaning an overall pattern to this relationship, which is a prerequisite to ultimately understanding the mechanisms by which past environments have shaped the evolutionary trajectory. Current obstacles include as-yet undecided reasons behind the frequent discrepancy between molecular and fossil time estimates, and the frequent lack of consideration of extensive confidence intervals around time estimates. We suggest that use and interpretation of both ancient DNA and molecular clocks is most effective when results are synthesized with palaeontological (fossil) and ecological (life history) information.

Relative age of proviral porcine endogenous retrovirus sequences in Sus scrofa based on the molecular clock hypothesis

Porcine endogenous retroviruses (PERV) are discussed as putative infectious agents in xenotransplantation. PERV classes A, B, and C harbor different envelope proteins. Two different types of long terminal repeat (LTR) structures exist, of which both are present only in PERV-A. One type of LTR contains a distinct repeat structure in U3, while the other is repeatless, conferring a lower level of transcriptional activity. Since the different LTR structures are distributed unequally among the proviruses and, apparently, PERV is the only virus harboring two different LTR structures, we were interested in determining which LTR is the ancestor. Replication-competent viruses can still be found today, suggesting an evolutionary recent origin. Our studies revealed that the age of PERV is at most 7.6 x 10(6) years, whereas the repeatless LTR type evolved approximately 3.4 x 10(6) years ago, being the phylogenetically younger structure. The age determined for PERV correlates with the time of separation between pigs (Suidae, Sus scrofa) and their closest relatives, American-born peccaries (Tayassuidae, Pecari tajacu), 7.4 x 10(6) years ago.

Making a robust biomolecular time scale for phylogenetic studies

The further evolution of informational molecular sequences should depend on the number of viable alternatives possible for the sequences as set by selection, the unrepaired mutation rate, and time. Most biomolecular clocks are based on Kimura's nearly neutral mutation random-drift hypothesis. This clock assumes that informational sequences are in equilibrium, i.e., the nucleotides mutate at a uniform rate and the number of nucleotides unconstrained by selection remains constant. Correcting for deviations from these assumptions should produce a more accurate clock. Informational molecules probably formed from polynucleotides having some other function such as nitrogen or nucleotide storage, thus being initially functionally unselected. At any time the rate of development of functionality in a protein may be expected to be proportional to the number of viable alternatives of sequence in its potentially interacting regions. Assuming the rate of unrepaired mutations is constant, these clocks should exponentially slow as they evolve, each with a different rate toward individual equilibria. Also if the degree of selection changes, its clock rate should change. For a more precise clock two approaches are suggested to estimate these time dependent changes in evolutionary rate. An improved clock could improve estimation of phylogeny and put a time scale on that phylogeny.

Similarity of phylogenetic trees as indicator of protein-protein interaction

Deciphering the network of protein interactions that underlines cellular operations has become one of the main tasks of proteomics and computational biology. Recently, a set of bioinformatics approaches has emerged for the prediction of possible interactions by combining sequence and genomic information. Even though the initial results are very promising, the current methods are still far from perfect. We propose here a new way of discovering possible protein-protein interactions based on the comparison of the evolutionary distances between the sequences of the associated protein families, an idea based on previous observations of correspondence between the phylogenetic trees of associated proteins in systems such as ligands and receptors. Here, we extend the approach to different test sets, including the statistical evaluation of their capacity to predict protein interactions. To demonstrate the possibilities of the system to perform large-scale predictions of interactions, we present the application to a collection of more than 67 000 pairs of E.coli proteins, of which 2742 are predicted to correspond to interacting proteins.

The origins of bioinformatics

Bioinformatics is often described as being in its infancy, but computers emerged as important tools in molecular biology during the early 1960s. A decade before DNA sequencing became feasible, computational biologists focused on the rapidly accumulating data from protein biochemistry. Without the benefits of super computers or computer networks, these scientists laid important conceptual and technical foundations for bioinformatics today.

Expansion and molecular evolution of the interferon-induced 2'-5' oligoadenylate synthetase gene family

The mammalian 2'-5' oligoadenylate synthetases (2'-5'OASs) are enzymes that are crucial in the interferon-induced antiviral response. They catalyze the polymerization of ATP into 2'-5'-linked oligoadenylates which activate a constitutively expressed latent endonuclease, RNaseL, to block viral replication at the level of mRNA degradation. A molecular evolutionary analysis of available OAS sequences suggests that the vertebrate genes are members of a multigene family with its roots in the early history of tetrapods. The modern mammalian 2'-5'OAS genes underwent successive gene duplication events resulting in three size classes of enzymes, containing one, two, or three homologous domains. Expansion of the OAS gene family occurred by whole-gene duplications to increase gene content and by domain couplings to produce the multidomain genes. Evolutionary analyses show that the 2'-5'OAS genes in rodents underwent gene duplications as recently as 11 MYA and predict the existence of additional undiscovered OAS genes in mammals.

A molecular timescale for vertebrate evolution

A timescale is necessary for estimating rates of molecular and morphological change in organisms and for interpreting patterns of macroevolution and biogeography. Traditionally, these times have been obtained from the fossil record, where the earliest representatives of two lineages establish a minimum time of divergence of these lineages. The clock-like accumulation of sequence differences in some genes provides an alternative method by which the mean divergence time can be estimated. Estimates from single genes may have large statistical errors, but multiple genes can be studied to obtain a more reliable estimate of divergence time. However, until recently, the number of genes available for estimation of divergence time has been limited. Here we present divergence-time estimates for mammalian orders and major lineages of vertebrates, from an analysis of 658 nuclear genes. The molecular times agree with most early (Palaeozoic) and late (Cenozoic) fossil-based times, but indicate major gaps in the Mesozoic fossil record. At least five lineages of placental mammals arose more than 100 million years ago, and most of the modern orders seem to have diversified before the Cretaceous/Tertiary extinction of the dinosaurs.

The meaning of near-neutrality at coding and non-coding regions

The nature of weak selection differs between coding and non-coding regions. Coding regions contain genetic information, whereas most non-coding regions do not have any information. Genetic information may be regarded as interaction systems, and the NK model of Kauffman was analysed. This model assumes that each amino acid makes a fitness contribution that depends on the amino acid and on K other amino acids among the N that make the protein. Through simulations, it was found that there are numerous nearly-neutral mutations under this model. Therefore, evolution is rapid in small populations, and slow in large populations. The variance of the evolutionary rate is not quite as large as data indicate under the model, and additional factors, such as environmental change or population-size fluctuation, need to be considered. Weak selection at non-coding regions may come from chromosome organization, and may be regional in character, which differs from that at coding regions. The problem of genetic load is thought to disappear in these circumstances.

Role of random genetic drift in the evolution of interactive systems

Through comparative studies of DNA sequences it has become possible to test the neutral and the selection theories of molecular evolution. The separate estimation of the numbers of synonymous and non-synonymous substitutions is one of the most powerful tools for detecting selection. The patterns on the average and variance of these two types of substitutions of mammalian genes turned out to be in accord with the slightly deleterious or the nearly neutral mutation theory for non-synonymous changes. Interactive systems at the amino acid level were suggested to be responsible for such nearly neutral or very weak selection. An attractive model is the NK model of Kauffman, which assumes that each amino acid makes a fitness contribution that depends upon the amino acid and upon K other amino acids among the N that make the protein. It is known that the fitness landscape is very rugged for K > or = 2. Population genetic analysis of this model suggest that protein evolution obeys the nearly neutral theory and that random genetic drift is important. In other words, evolution becomes rapid in small populations because proportion of near-neutrality increases among new mutations, and proteins as interactive systems evolve by shifting through random genetic drift on the multipeaked fitness landscape.

Molecular evolution and secondary structural conservation in the B-cell lymphoma leukemia 2 (bcl-2) family of proto-oncogene products

The nature of the bcl-2 family of proto-oncogenes was analyzed by sequence alignment, secondary structure prediction, and phylogenetic techniques. Phylogenies were inferred from both the nucleic acid and amino acid sequences of the human, murine, rat, and chicken sequences for BCL-2 and BCL-X, human MCL1, murine A1, the nematode Caenorhabditis elegans and Caenorhabditis briggsiae ced-9 proteins, and the sequences BHRF1 from Epstein-Barr and LMW5-HL from African swine fever viruses. Both sequence alignment and secondary structure prediction techniques supported the conservation of both the overall secondary structure and the carboxy-terminal transmembrane domain in all members of the family. All the treeing methods employed (distance matrix, maximum likelihood, and parsimony) supported a tree in which the proapoptotic proteins BCL-2 and BCL-X represent the most recent additions to the group. All the trees also indicated that the viral proteins BHRF1 and LMW-HL arose from a common ancestor, an ancestor they shared in common with the pro-apoptotic control protein BAX, indicating that this function of BAX evolved only recently. The most ancient branches are represented by the nematode ced-9 protein and by the control genes MCL1 and A1, which in the treeing methods employed represent separate lineages within the most ancient grouping. These results demonstrate the evolution of a highly conserved family of developmental control genes from nematode to man--genes that encode proteins essential for normal development but which are highly conserved in terms of predicted structure and possible cellular localization. The evolutionary analysis also indicates that the family may be even larger than originally predicted and that other members are waiting to be discovered.

T4 phage evolution data in terms of a time-dependent Topal-Fresco mechanism

A physical interpretation of the Topal-Fresco [Nature 263, 285 (1976)] model for spontaneous base substitutions suggests that hydrogen-bonded DNA protons satisfy the criteria for a classical noninteracting isolated system. Accessible states for duplex G-C protons include the keto-amino state and the six complementary enol-imine isomers. Hydrogen-bonded enol and imine protons occupy symmetric double-minima created by the two sets of indistinguishable electron lone pairs and a single proton belonging to each enol-imine end group. These protons will consequently participate in coupled quantum mechanical flip-flop, tunneling back and forth between symmetric energy wells. This results in a quantum mixing of proton energy states where the lowest energy state will be a linear combination of available G-C isomers. The resulting conclusion is that metastable keto-amino G-C protons will populate accessible enol-imine stationary states at rates governed by quantum laws of statistical equilibrium, consistent with achieving the lowest energy condition for duplex G-C protons. Enol-imine G-C stationary states are bound more tightly, of the order of 3 to 12 kcal/mol, which requires a modified mode of Topal-Fresco replication that will inhibit reequilibration of enol and imine G and C template isomers and, thus, promote the formation of complementary mispairs. The model is demonstrated on time-dependent base substitutions expressed by T4 phage DNA systems where data are consistent with model explanations, including the prediction that time-dependent evolutionary transversion sites will exhibit both G-C-to-T-A and G-C-to-C-G transversions at replication, due to proton flip-flop alteration of G template genetic specificity. The observation that A-T sites are resistant to time-dependent evolutionary base substitutions, expressed exclusively at G-C sites, allows codons to be classified as either evolutionary sensitive (16 codons) or evolutionary resistant (8 codons). These criteria provide possible explanations for expansion properties of the CGG fragile X sequences. Enol-imine G-C stationary states appear to have been misdiagnosed as deamination of cytosine and oxidation of guanine to 8-hydroxy-guanine.

Molecular evolution of interleukin-3

Chimpanzee, tamarin, and marmoset interleukin-3 (IL-3) genes were cloned, sequenced, and expressed. Western blot analysis demonstrated that functional genes were isolated. IL-3 sequences were compared with those of mouse, rat, rhesus monkey, gibbon, and man. Multiple alignment of the IL-3 coding regions showed that only a few regions had been conserved during mammalian evolution, which are likely associated with functional domains of the IL-3 protein. Substitution rates for the various lineages were calculated and the numbers of synonymous and nonsynonymous substitutions were estimated separately. Distance matrices of the IL-3 coding regions were used to construct phylogenetic trees which revealed large differences in IL-3 evolution rate as well as a more rapid substitution rate for rodents and a rate slowdown during hominoid evolution. Extremes were rhesus monkey IL-3, which accumulated few synonymous substitutions, and gibbon IL-3, which had almost exclusively synonymous substitutions. In rhesus monkey IL-3, nonsynonymous substitutions outnumbered synonymous substitutions, which could not be readily explained by a random process of substitutions. We assume that during evolution of IL-3, the majority of the amino acid replacements and the impaired interspecies functional cross-reactivity originate from selection mechanisms with the most likely selective force being the structure of the heterodimeric IL.3 cell-surface receptor. Insight into IL-3 architecture and structural analysis of the IL-3 receptor are needed to analyze the unusually fast evolution of IL-3 in more detail.

Evolution of translational elongation factor (EF) sequences: reliability of global phylogenies inferred from EF-1 alpha(Tu) and EF-2(G) proteins

The EF-2 coding genes of the Archaea Pyrococcus woesei and Desulfurococcus mobilis were cloned and sequenced. Global phylogenies were inferred by alternative tree-making methods from available EF-2(G) sequence data and contrasted with phylogenies constructed from the more conserved but shorter EF-1 alpha(Tu) sequences. Both the monophyly (sensu Hennig) of Archaea and their subdivision into the kingdoms Crenarchaeota and Euryarchaeota are consistently inferred by analysis of EF-2(G) sequences, usually at a high bootstrap confidence level. In contrast, EF-1 alpha(Tu) phylogenies tend to be inconsistent with one another and show low bootstrap confidence levels. While evolutionary distance and DNA maximum parsimony analyses of EF-1 alpha(Tu) sequences do show archaeal monophyly, protein parsimony and DNA maximum-likelihood analyses of these data do not. In no case, however, do any of the tree topologies inferred from EF-1 alpha(Tu) sequence analyses receive significant bootstrap support.

DNA-dependent RNA polymerase subunit B as a tool for phylogenetic reconstructions: branching topology of the archaeal domain

The branching topology of the archaeal (archaebacterial) domain was inferred from sequence comparisons of the largest subunit (B) of DNA-dependent RNA polymerases (RNAP). Both the nucleic acid sequences of the genes coding for RNAP subunit B and the amino acid sequences of the derived gene products were used for phylogenetic reconstructions. Individual analysis of the three nucleotide positions of codons revealed significant inequalities with respect to guanosine and cytosine (GC) content and evolutionary rates. Only the nucleotides at the second codon positions were found to be unbiased by varied GC contents and sufficiently conserved for reliable phylogenetic reconstructions. A decision matrix was used for the combination of the results of distance matrix, maximum parsimony, and maximum likelihood methods. For this purpose the original results (sums of squares, steps, and logarithms of likelihoods) were transformed into comparable effective values and analyzed with methods known from the theory of statistical decisions. Phylogenetic invariants and statistical analysis with resampling techniques (bootstrap and jackknife) confirmed the preferred branching topology, which is significantly different from the topology known from phylogenetic trees based on 16S rRNA sequences. The preferred topology reconstructed by this analysis shows a common stem for the Methanococcales and Methanobacteriales and a separation of the thermophilic sulfur archaea from the methanogens and halophiles. The latter coincides with a unique phylogenetic location of a characteristic splitting event replacing the largest RNAP subunit of thermophilic sulfur archaea by two fragments in methanogens and halophiles. This topology is in good agreement with physiological and structural differences between the various archaea and demonstrates RNAP to be a suitable phylogenetic marker molecule.

The DNA-dependent RNA-polymerase of Thermotoga maritima; characterisation of the enzyme and the DNA-sequence of the genes for the large subunits

An improved purification procedure for Thermotoga maritima RNA-polymerase holoenzyme was developed. The enzyme is highly active with poly dAT or T7 phage DNA as template. DNA gyrase was found to be a side product of this RNA-polymerase purification. The genes for the large subunits beta and beta' of RNA-polymerase were cloned and sequenced. The phylogenetic position of T.maritima within the bacterial domain was determined by various methods. It is the lowest bacterial offspring but slightly higher than the chloroplasts.

Molecular analysis of ependymins from the cerebrospinal fluid of the orders Clupeiformes and Salmoniformes: no indication for the existence of an euteleost infradivision

Ependymins represent the predominant protein constituents in the cerebrospinal fluid of many teleost fish and they are synthesized in meningeal fibroblasts. Here, we present the ependymin sequences from the herring (Clupea harengus) and the pike (Esox lucius). A comparison of ependymin homologous sequences from three different orders of teleost fish (Salmoniformes, Cypriniformes, and Clupeiformes) revealed the highest similarity between Clupeiformes and Cypriniformes. This result is unexpected because it does not reflect current systematics, in which Clupeiformes belong to a separate infradivision (Clupeomorpha) than Salmoniformes and Cypriniformes (Euteleostei). Furthermore, in Salmoniformes the evolutionary rate of ependymins seems to be accelerated mainly on the protein level. However, considering these inconstant rates, neither neighbor-joining trees nor DNA parsimony methods gave any indication that a separate euteleost infradivision exists.

Cloning and sequencing of the gene encoding glutamine synthetase I from the archaeum Pyrococcus woesei: anomalous phylogenies inferred from analysis of archaeal and bacterial glutamine synthetase I sequences

The gene glnA encoding glutamine synthetase I (GSI) from the archaeum Pyrococcus woesei was cloned and sequenced with the Sulfolobus solfataricus glnA gene as the probe. An operon reading frame of 448 amino acids was identified within a DNA segment of 1,528 bp. The encoded protein was 49% identical with the GSI of Methanococcus voltae and exhibited conserved regions characteristic of the GSI family. The P. woesei GSI was aligned with available homologs from other archaea (S. solfataricus, M. voltae) and with representative sequences from cyanobacteria, proteobacteria, and gram-positive bacteria. Phylogenetic trees were constructed from both the amino acid and the nucleotide sequence alignments. In accordance with the sequence similarities, archaeal and bacterial sequences did not segregate on a phylogeny. On the basis of sequence signatures, the GSI trees could be subdivided into two ensembles. One encompassed the GSI of cyanobacteria and proteobacteria, but also that of the high-G + C gram-positive bacterium Streptomyces coelicolor (all of which are regulated by the reversible adenylylation of the enzyme subunits); the other embraced the GSI of the three archaea as well as that of the low-G + C gram-positive bacteria (Clostridium acetobutilycum, Bacillus subtilis) and Thermotoga maritima (none of which are regulated by subunit adenylylation). The GSIs of the Thermotoga and the Bacillus-Clostridium lineages shared a direct common ancestor with that of P. woesei and the methanogens and were unrelated to their homologs from cyanobacteria, proteobacteria, and S. coelicolor. The possibility is presented that the GSI gene arose among the archaea and was then laterally transferred from some early methanogen to a Thermotoga-like organism. However, the relationship of the cyanobacterial-proteobacterial GSIs to the Thermotoga GSI and the GSI of low-G+C gram-positive bacteria remains unexplained.

Sorghum phosphoenolpyruvate carboxylase gene family: structure, function and molecular evolution

Although housekeeping functions have been shown for the phosphoenolpyruvate carboxylase (EC 4.1.1.31, PEPC) in plants and in prokaryotes, PEPC is mainly known for its specific role in the primary photosynthetic CO2 fixation in C4 and CAM plants. We have shown that in Sorghum, a monocotyledonous C4 plant, the enzyme is encoded in the nucleus by a small multigene family. Here we report the entire nucleotide sequence (7.5 kb) of the third member (CP21) that completes the structure of the Sorghum PEPC gene family. Nucleotide composition, CpG islands and GC content of the three Sorghum PEPC genes are analysed with respect to their possible implications in the regulation of expression. A study of structure/function and phylogenetic relationships based on the compilation of all PEPC sequences known so far is presented. Data demonstrated that: (1) the different forms of plant PEPC have very similar primary structures, functional and regulatory properties, (2) neither apparent amino acid sequences nor phylogenetic relationships are specific for the C4 and CAM PEPCs and (3) expression of the different genes coding for the Sorghum PEPC isoenzymes is differently regulated (i.e. by light, nitrogen source) in a spatial and temporal manner. These results suggest that the main distinguishing feature between plant PEPCs is to be found at the level of genes expression rather than in their primary structure.

Early evolutionary relationships among known life forms inferred from elongation factor EF-2/EF-G sequences: phylogenetic coherence and structure of the archaeal domain

Phylogenies were inferred from both the gene and the protein sequences of the translational elongation factor termed EF-2 (for Archaea and Eukarya) and EF-G (for Bacteria). All treeing methods used (distance-matrix, maximum likelihood, and parsimony), including evolutionary parsimony, support the archaeal tree and disprove the "eocyte tree" (i.e., the polyphyly and paraphyly of the Archaea). Distance-matrix trees derived from both the amino acid and the DNA sequence alignments (first and second codon positions) showed the Archaea to be a monophyletic-holophyletic grouping whose deepest bifurcation divides a Sulfolobus branch from a branch comprising Methanococcus, Halobacterium, and Thermoplasma. Bootstrapped distance-matrix treeing confirmed the monophyly-holophyly of Archaea in 100% of the samples and supported the bifurcation of Archaea into a Sulfolobus branch and a methanogen-halophile branch in 97% of the samples. Similar phylogenies were inferred by maximum likelihood and by maximum (protein and DNA) parsimony. DNA parsimony trees essentially identical to those inferred from first and second codon positions were derived from alternative DNA data sets comprising either the first or the second position of each codon. Bootstrapped DNA parsimony supported the monophyly-holophyly of Archaea in 100% of the bootstrap samples and confirmed the division of Archaea into a Sulfolobus branch and a methanogen-halophile branch in 93% of the bootstrap samples. Distance-matrix and maximum likelihood treeing under the constraint that branch lengths must be consistent with a molecular clock placed the root of the universal tree between the Bacteria and the bifurcation of Archaea and Eukarya. The results support the division of Archaea into the kingdoms Crenarchaeota (corresponding to the Sulfolobus branch and Euryarchaeota). This division was not confirmed by evolutionary parsimony, which identified Halobacterium rather than Sulfolobus as the deepest offspring within the Archaea.

Revisiting junk DNA

The distribution of functions within genomes of higher organisms relative to processes that lead to the spread of mutations in populations is examined in its general outlines. A number of points are enumerated that collectively put in question the concept of junk DNA: the plausible compatibility of DNA function with rapid substitution rates; the likelihood of superimposed functions along much of eukaryotic DNA; the potential for a merely conditional functionality in sequence repeats; the apparent adoption of macromolecular waste as a strategy for maintaining a function without selective grooming of individual sequence repeats that carry out the function; the likely requirement that any DNA sequence must be "polite" vis-'a-vis (compatible with) functional sequences in its genomic environment; the existence in germ-cell lineages of selective constraints that are not apparent in populations of individuals; and the fact that DNA techtonics - the appearance and disappearance of genomic DNA - are not incompatible with function. It is pointed out that the inverse correlation between functional constraints and rates of substitution cannot be claimed to be pillar of the neutral theory, because it is also predicted from a selectionist viewpoint. The dispensability of functional structures is brought into relation with the concept of reproductive sufficiency the survivability of genotypes in the absence of fitter alleles.

Uniformity in the nonsynonymous substitution rates of embryonic beta-globin genes of several vertebrate species

The nucleotide substitution rate in structural portions of the embryonic beta-globin genes of placental mammals is lower than that for the adult beta-globin genes. This difference occurs entirely within the class of substitutions that result in nonsynonymous (replacement) differences between these genes, and therefore represents a constraint on the structure of the mammalian embryonic beta-globin proteins relative to the adult proteins (Shapiro et al. 1983; Hardison 1984). A similar effect has also been observed in marsupial mammals (Koop and Goodman 1988). In an effort to determine whether the observed rates are evidence of a uniform degree of selective constraint on the embryonic beta-globin genes, analyses were performed that compared replacement substitution rates. The analyses reveal that embryonic beta-globin genes appear to have been fixing replacement substitutions at nearly the same average rate not only in placental and marsupial mammals but in avian and amphibian species as well. In contrast, the adult beta-globin genes from these organisms appear to have a more variable rate of replacement substitution with an especially low rate for birds. In the chicken (Gallus gallus), the adult beta-globin gene replacement substitution rate appears to be lower than the embryonic replacement substitution rate.

Glutamine synthetase gene evolution: a good molecular clock

Glutamine synthetase (EC 6.3.1.2) gene evolution in various animals, plants, and bacteria was evaluated by a general stationary Markov model. The evolutionary process proved to be unexpectedly regular even for a time span as long as that between the divergence of prokaryotes from eukaryotes. This enabled us to draw phylogenetic trees for species whose phylogeny cannot be easily reconstructed from the fossil record. Our calculation of the times of divergence of the various organelle-specific enzymes led us to hypothesize that the pea and bean chloroplast genes for these enzymes originated from the duplication of nuclear genes as a result of the different metabolic needs of the various species. Our data indicate that the duplication of plastid glutamine synthetase genes occurred long after the endosymbiotic events that produced the organelles themselves.