Response of stomatal numbers to CO2 and humidity: control by transpiration rate and abscisic acid

The observation that stomatal density (number mm(-2)) on herbarium leaves had decreased over the last century represents clear evidence that plants have responded to anthropogenic increases in CO2 concentration. The mechanism of the response has proved elusive but here it is shown that density responses to both CO2 concentration and humidity are correlated with changes in whole-plant transpiration and leaf abscisic acid (ABA) concentration. The transpiration rate of a range of accessions of Arabidopsis thaliana was manipulated by changing CO2 concentration, humidity and by exogenous application of ABA. Stomatal density increased with transpiration and leaf ABA concentration. A common property of signal transduction systems is that they rapidly lose their ability to respond to the co-associated stimulus. Pathways of water movement within the plant are connected and so variations in supply and demand can be signalled throughout the plant directly, modifying stomatal aperture of mature leaves and stomatal density of developing leaves. Furthermore, the system identified here does not conform to the loss of ability to respond. A putative mechanism is proposed for the control of stomatal density by transpiration rate and leaf ABA concentration.

Systemic signalling of environmental cues in Arabidopsis leaves

Light intensity and atmospheric CO2 partial pressure are two environmental signals known to regulate stomatal numbers. It has previously been shown that if a mature Arabidopsis leaf is supplied with either elevated CO2 (750 ppm instead of ambient at 370 ppm) or reduced light levels (50 micromol m-2 s-1 instead of 250 micromol m-2 s-1), the young, developing leaves that are not receiving the treatment grow with a stomatal density as if they were exposed to the treatment. But the signal(s) that it is believed is generated in the mature leaves and transmitted to developing leaves are largely unknown. Photosynthetic rates of treated, mature Arabidopsis leaves increased in elevated CO2 and decreased when shaded, as would be expected. Similarly, the levels of sugars (glucose, fructose, and sucrose) in the treated mature leaves increased in elevated CO2 and decreased with shade treatment. The levels of sugar in developing leaves were also measured and it was found that they mirrored this result even though they were not receiving the shade or elevated CO2 treatment. To investigate the effect of these treatments on global gene expression patterns, transcriptomics analysis was carried out using Affymetrix, 22K, and ATH1 arrays. Total RNA was extracted from the developing leaves after the mature leaves had received either the ambient control treatment, the elevated CO2 treatment, or the shade treatment, or both elevated CO2 and shade treatments for 2, 4, 12, 24, 48, or 96 h. The experiment was replicated four times. Two other experiments were also conducted, one to compare and contrast gene expression in response to plants grown at elevated CO2 and the other to look at the effect of these treatments on the mature leaf. The data were analysed and 915 genes from the untreated, signalled leaves were identified as having expression levels affected by the shade treatment. These genes were then compared with those whose transcript abundance was affected by the shade treatment in the mature treated leaves (1181 genes) and with 220 putative 'stomatal signalling' genes previously identified from studies of the yoda mutant. The results of these experiments and how they relate to environmental signalling are discussed, as well as possible mechanisms for systemic signalling.

An ancestral nuclear protein assembly: crystal structure of the Methanopyrus kandleri histone

Eukaryotic histone proteins condense DNA into compact structures called nucleosomes. Nucleosomes were viewed as a distinguishing feature of eukaryotes prior to identification of histone orthologs in methanogens. Although evolutionarily distinct from methanogens, the methane-producing hyperthermophile Methanopyrus kandleri produces a novel, 154-residue histone (HMk). Amino acid sequence comparisons show that HMk differs from both methanogenic and eukaryotic histones, in that it contains two histone-fold ms within a single chain. The two HMk histone-fold ms, N and C terminal, are 28% identical in amino acid sequence to each other and approximately 21% identical in amino acid sequence to other histone proteins. Here we present the 1.37-A-resolution crystal structure of HMk and report that the HMk monomer structure is homologous to the eukaryotic histone heterodimers. In the crystal, HMk forms a dimer homologous to [H3-H4](2) in the eukaryotic nucleosome. Based on the spatial similarities to structural ms found in the eukaryotic nucleosome that are important for DNA-binding, we infer that the Methanopyrus histone binds DNA in a manner similar to the eukaryotic histone tetramer [H3-H4](2).

A type IB topoisomerase with DNA repair activities

Previously we have characterized type IB DNA topoisomerase V (topo V) in the hyperthermophile Methanopyrus kandleri. The enzyme has a powerful topoisomerase activity and is abundant in M. kandleri. Here we report two characterizations of topo V. First, we found that its N-terminal domain has sequence homology with both eukaryotic type IB topoisomerases and the integrase family of tyrosine recombinases. The C-terminal part of the sequence includes 12 repeats, each repeat consisting of two similar but distinct helix-hairpin-helix motifs; the same arrangement is seen in recombination protein RuvA and mammalian DNA polymerase beta. Second, on the basis of sequence homology between topo V and polymerase beta, we predict and demonstrate that topo V possesses apurinic/apyrimidinic (AP) site-processing activities that are important in base excision DNA repair: (i) it incises the phosphodiester backbone at the AP site, and (ii) at the AP endonuclease cleaved AP site, it removes the 5' 2-deoxyribose 5-phosphate moiety so that a single-nucleotide gap with a 3'-hydroxyl and 5'-phosphate can be filled by a DNA polymerase. Topo V is thus the prototype for a new subfamily of type IB topoisomerases and is the first example of a topoisomerase with associated DNA repair activities.

Isotope fractionation and atmospheric oxygen: implications for phanerozoic O(2) evolution

Models describing the evolution of the partial pressure of atmospheric oxygen over Phanerozoic time are constrained by the mass balances required between the inputs and outputs of carbon and sulfur to the oceans. This constraint has limited the applicability of proposed negative feedback mechanisms for maintaining levels of atmospheric O(2) at biologically permissable levels. Here we describe a modeling approach that incorporates O(2)-dependent carbon and sulfur isotope fractionation using data obtained from laboratory experiments on carbon-13 discrimination by vascular land plants and marine plankton. The model allows us to calculate a Phanerozoic O(2) history that agrees with independent models and with biological and physical constraints and supports the hypothesis of a high atmospheric O(2) content during the Carboniferous (300 million years ago), a time when insect gigantism was widespread.

Horizontal gene transfer among genomes: the complexity hypothesis

Increasingly, studies of genes and genomes are indicating that considerable horizontal transfer has occurred between prokaryotes. Extensive horizontal transfer has occurred for operational genes (those involved in housekeeping), whereas informational genes (those involved in transcription, translation, and related processes) are seldomly horizontally transferred. Through phylogenetic analysis of six complete prokaryotic genomes and the identification of 312 sets of orthologous genes present in all six genomes, we tested two theories describing the temporal flow of horizontal transfer. We show that operational genes have been horizontally transferred continuously since the divergence of the prokaryotes, rather than having been exchanged in one, or a few, massive events that occurred early in the evolution of prokaryotes. In agreement with earlier studies, we found that differences in rates of evolution between operational and informational genes are minimal, suggesting that factors other than rate of evolution are responsible for the observed differences in horizontal transfer. We propose that a major factor in the more frequent horizontal transfer of operational genes is that informational genes are typically members of large, complex systems, whereas operational genes are not, thereby making horizontal transfer of informational gene products less probable (the complexity hypothesis).

Formation of a primitive ectoderm like cell population, EPL cells, from ES cells in response to biologically derived factors

The primitive ectoderm of the mouse embryo arises from the inner cell mass between 4.75 and 5.25 days post coitum, around the time of implantation. Positioned at a pivotal time in development, just prior to formation of the three germ layers of the embryo proper, the primitive ectoderm responds directly to the signals generated during gastrulation. We have identified a conditioned medium, MEDII, which caused the homogeneous conversion of ES cells to a morphologically distinct cell population, termed early primitive ectoderm-like (EPL) cells. EPL cells expressed the pluripotent cell markers Oct4, SSEA1 and alkaline phosphatase. However, the formation of EPL cells was accompanied by alterations in Fgf5, Gbx2 and Rex1 expression, a loss in chimaera forming ability, changes in factor responsiveness and modified differentiation capabilities, all consistent with the identification of EPL cells as equivalent to the primitive ectoderm population of the 5.5 to 6.0 days post coitum embryo. EPL cell formation could be reversed in the presence of LIF and withdrawal of MEDII, which suggested that EPL cell formation was not a terminal differentiation event but reflected the ability of pluripotent cells to adopt distinct cell states in response to specific factors. Partial purification of MEDII revealed the presence of two separable biological activities, both of which were required for the induction and maintenance of EPL cells. We show here the first demonstration of uniform differentiation of ES cells in response to biological factors. The formation of primitive ectoderm, both in vivo and in vitro, appears to be an obligatory step in the differentiation of the inner cell mass or ES cells into cell lineages of the embryonic germ layers. EPL cells potentially represent a model for the development of lineage specific differentiation protocols and analysis of gastrulation at a molecular level. An understanding of the active components of MEDII may provide a route for the identification of factors which induce primitive ectoderm formation in vivo.

Optimally recovering rate variation information from genomes and sequences: pattern filtering

Nucleotide substitution rates vary at different positions within genes and genomes, but rates are difficult to estimate, because they are masked by the stochastic nature of substitutions. In this paper, a linear method, pattern filtering, is described which can optimally separate the signals (related to substitution rates or to other measures of sequence change) from stochastic noise. Pattern filtering promises to be useful in both genomic and molecular evolution studies. In an example using mitochondrial genomes, it is shown that pattern filtering can reveal coding and non-coding regions without the need for prior identification of reading frames or other knowledge of the sequence and promises to be an important tool for genomic analysis. In a second example, it is shown that pattern filtering allows one to classify sites on the basis of an estimator of substitution rates. Using elongation factor EF-1 alpha sequences, it is shown that the fastest sites favor archaea as the sister taxon of eukaryotes, whereas the slower sites support the eocyte prokaryotes as the sister taxon of eukaryotes, suggesting that the former result is an artifact of "long branch attraction."

Genomic evidence for two functionally distinct gene classes

Analyses of complete genomes indicate that a massive prokaryotic gene transfer (or transfers) preceded the formation of the eukaryotic cell. In comparisons of the entire set of Methanococcus jannaschii genes with their orthologs from Escherichia coli, Synechocystis 6803, and the yeast Saccharomyces cerevisiae, it is shown that prokaryotic genomes consist of two different groups of genes. The deeper, diverging informational lineage codes for genes which function in translation, transcription, and replication, and also includes GTPases, vacuolar ATPase homologs, and most tRNA synthetases. The more recently diverging operational lineage codes for amino acid synthesis, the biosynthesis of cofactors, the cell envelope, energy metabolism, intermediary metabolism, fatty acid and phospholipid biosynthesis, nucleotide biosynthesis, and regulatory functions. In eukaryotes, the informational genes are most closely related to those of Methanococcus, whereas the majority of operational genes are most closely related to those of Escherichia, but some are closest to Methanococcus or to Synechocystis.

Evidence for an early prokaryotic origin of histones H2A and H4 prior to the emergence of eukaryotes

Histones have been identified recently in many prokaryotes. These histones, unlike their eukaryotic homologs, are of a single uniform type that is thought to resemble the archetypal ancestor of the eukaryotic histone family. In this paper we report the finding, the cloning and the phylogenetic analysis of the sequence of a prokaryotic histone from the hyperthermophile Methanopyrus kandleri . Unlike previously described prokaryotic histones, the Methanopyrus sequence has a novel structure consisting of two tandemly repeated histone fold motifs in a single polypeptide. Sequence analyses indicate that the N-terminal repeat is most closely related to eukaryotic H2A and H4 histones, whereas the C-terminal repeat resembles that found in prokaryotic histones. These results imply an early divergence within the histone gene family prior to the emergence of eukaryotes and may represent an evolutionary step leading to eukaryotic histones.

Evidence for a clade of nematodes, arthropods and other moulting animals

The arthropods constitute the most diverse animal group, but, despite their rich fossil record and a century of study, their phylogenetic relationships remain unclear. Taxa previously proposed to be sister groups to the arthropods include Annelida, Onychophora, Tardigrada and others, but hypotheses of phylogenetic relationships have been conflicting. For example, onychophorans, like arthropods, moult periodically, have an arthropod arrangement of haemocoel, and have been related to arthropods in morphological and mitochondrial DNA sequence analyses. Like annelids, they possess segmental nephridia and muscles that are a combination of smooth and obliquely striated fibres. Our phylogenetic analysis of 18S ribosomal DNA sequences indicates a close relationship between arthropods, nematodes and all other moulting phyla. The results suggest that ecdysis (moulting) arose once and support the idea of a new clade, Ecdysozoa, containing moulting animals: arthropods, tardigrades, onychophorans, nematodes, nematomorphs, kinorhynchs and priapulids. No support is found for a clade of segmented animals, the Articulata, uniting annelids with arthropods. The hypothesis that nematodes are related to arthropods has important implications for developmental genetic studies using as model systems the nematode Caenorhabditis elegans and the arthropod Drosophila melanogaster, which are generally held to be phylogenetically distant from each other.

Phylogenetic inference: how much evolutionary history is knowable?

In order to reconstruct phylogenetic trees from extremely dissimilar sequences it is necessary to estimate accurately the extent of sequence divergence. In this paper a new method of sequence analysis, Markov triple analysis, is developed for determining the relative frequencies of nucleotide substitutions within the three branches of a three-taxon dendrogram. Assuming that nucleotide sites are independently and identically distributed and assuming a Markov model for nucleotide (or protein) evolution, it is shown that the unique Markov matrices can be reconstructed given only the joint probability distribution relating three taxa. (In the much simpler case involving only two taxa and two character states, Markov matrices can also be reconstructed, provided symmetry assumptions are placed on the elements of the matrices.) The method is illustrated using sequence data from the combined first and second codon positions derived from complete human, mouse, and cow mitochondrial sequences.

Inference for phylogenies under a hybrid parsimony method: evolutionary-symmetric transversion parsimony

A new method is proposed for inferring topology for evolutionary trees. Existing methods have complementary strengths and weaknesses. Maximum and transversion parsimony are powerful methods, but they lack statistical consistency, that is, they do not always infer the correct tree as the sequence length becomes very large. Evolutionary parsimony overcomes this deficiency, but it may lack sufficient power when sequence length is small (less than 1000 aligned nucleotides; Sinsheimer, Lake, and Little, 1996, Biometrics 52, 193-210). Our proposed method, evolutionary-symmetric transversion parsimony, is a hybrid that retains the consistency of evolutionary parsimony, while increasing power by incorporating a modified form of transversion parsimony within a statistical model. The method requires choice of a parameter gamma that represents the prior probability that symmetric transversion parsimony yields consistent results. Properties of the method are assessed for a variety of choices of gamma in a large simulation study. In general, inference under the evolutionary-symmetric transversion parsimony has more discriminating power than inference under evolutionary parsimony and is better calibrated than inference under symmetric transversion parsimony. The results are quite robust to the choice of gamma, indicating a value of 0.90 as a reasonable overall choice when the true value of gamma ranges between 0.85 to 1.00. Our method is, like evolutionary parsimony and maximum parsimony, computationally straightforward. The same statistical approach can be applied to combine evolutionary parsimony with other inconsistent methods, such as maximum parsimony, but at the expense of more difficult computations.

Phylogeny of Methanopyrus kandleri based on methyl coenzyme M reductase operons

The mcrBDCGA operon that encodes methyl coenzyme M reductase (MR) in the hyperthermophile Methanopyrus kandleri was cloned and sequenced. The results of a phylogenetic analysis of the nine MR sequences now available support the position that M. kandleri is a separate methanogen lineage. As in other methanogens, the M. kandleri mcr operon is located immediately upstream of the mtrE gene, the promoter-proximal gene in an operon that encodes the N5-methyltetrahydromethanopterin:coenzyme M methyltransferase that catalyzes the step preceding the MR-catalyzed reaction in methanogenesis. In contrast to other methanogens and hyperthermophilic members of the Archaea, CG dinucleotides and CG-containing codons occur frequently in M. kandleri DNA. The MR subunit-encoding genes are preceded by sequences consistent with ribosome binding sites, indicating that mRNA-rRNA base pairing can still direct translation initiation in cells growing at temperatures above 100 degrees C.

Bayesian hypothesis testing of four-taxon topologies using molecular sequence data

The reconstruction of phylogenetic trees from molecular sequences presents unusual problems for statistical inference. For example, three possible alternatives must be considered for four taxa when inferring the correct unrooted tree (referred to as a topology). In our view, classical hypothesis testing is poorly suited to this triangular set of alternative hypotheses. In this article, we develop Bayesian inference to determine the posterior probability that a four-taxon topology is correct given the sequence data and the evolutionary parsimony algorithm for phylogenetic reconstruction. We assess the frequency properties of our models in a large simulation study. Bayesian inference under the principles of evolutionary parsimony is shown to be well calibrated with reasonable discriminating power for a wide range of realistic conditions, including conditions that violate the assumptions of evolutionary parsimony.

The phylogeny of Methanopyrus kandleri

The phylogenetic position of Methanopyrus kandleri has been difficult to determine because reconstructions of phylogenetic trees based on rRNA sequences have been ambiguous. The most probable trees determined by most algorithms place the genus Methanopyrus at the base of a group that includes the halobacteria and the methanogens and their relatives, although occasionally some algorithms place this genus near the eocytes (the hyperthermophilic, sulfur-metabolizing prokaryotes), suggesting that it may belong to this lineage. In order to resolve the phylogeny of the genus Methanopyrus, we determined the sequence of an informative region of elongation factor 1-alpha that contains an 11-amino-acid insertion in eocytes and eukaryotes which is replaced by a 4-amino-acid insertion in methanogens, halobacteria, and eubacteria. On the basis of the results of our elongation factor 1-alpha gene analysis, we concluded that the genus Methanopyrus diverged from the eocyte branch before the eukaryotic and eocyte lineages separated and therefore is not an eocyte.

Calculating the probability of multitaxon evolutionary trees: bootstrappers Gambit

The reconstruction of multitaxon trees from molecular sequences is confounded by the variety of algorithms and criteria used to evaluate trees, making it difficult to compare the results of different analyses. A global method of multitaxon phylogenetic reconstruction described here, Bootstrappers Gambit, can be used with any four-taxon algorithm, including distance, maximum likelihood, and parsimony methods. It incorporates a Bayesian-Jeffreys'-bootstrap analysis to provide a uniform probability-based criterion for comparing the results from diverse algorithms. To examine the usefulness of the method, the origin of the eukaryotes has been investigated by the analysis of ribosomal small subunit RNA sequences. Three common algorithms (paralinear distances, Jukes-Cantor distances, and Kimura distances) support the eocyte topology, whereas one (maximum parsimony) supports the archaebacterial topology, suggesting that the eocyte prokaryotes are the closest prokaryotic relatives of the eukaryotes.

DNA enzymology above 100 degrees C. Topoisomerase V unlinks circular DNA at 80-122 degrees C

The widespread application of polymerase chain reaction and related techniques in biology and medicine has led to a heightened interest in thermophilic enzymes of DNA metabolism. Some of these enzymes are stable for hours at 100 degrees C, but no enzymatic activity on duplex DNA at temperatures above 100 degrees C has so far been demonstrated. Recently, we isolated topoisomerase V from the hyperthermophile Methanopyrus kandleri, which grows up to 110 degrees C. This novel enzyme is similar to eukaryotic topoisomerase I and acts on duplex DNA regions. We now show that topoisomerase V catalyzes the unlinking of double-stranded circular DNA at temperatures up to 122 degrees C. In this in vitro system, maximal DNA unlinking occurs at 108 degrees C and corresponds to complementary strands being linked at most once. These results further imply that in the presence of sufficient positive supercoiling DNA can exist as a double helix even at 122 degrees C.

Evidence from 18S ribosomal DNA that the lophophorates are protostome animals

The suspension-feeding metazoan subkingdom Lophophorata exhibits characteristics of both deuterostomes and protostomes. Because the morphology and embryology of lophophorates are phylogenetically ambiguous, their origin is a major unsolved problem of metazoan phylogenetics. The complete 18S ribosomal DNA sequences of all three lophophorate phyla were obtained and analyzed to clarify the phylogenetic relationships of this subkingdom. Sequence analyses show that lophophorates are protostomes closely related to mollusks and annelids. This conclusion deviates from the commonly held view of deuterostome affinity.

A reverse gyrase with an unusual structure. A type I DNA topoisomerase from the hyperthermophile Methanopyrus kandleri is a two-subunit protein

Reverse gyrase, an ATP-dependent topoisomerase that positively supercoils DNA, has been purified to near-homogeneity from the hyperthermophile Methanopyrus kandleri. It migrates on SDS-polyacrylamide gel electrophoresis as two principal bands with apparent molecular masses of 150 and 50 kDa. Both proteins remain associated throughout all chromatographic steps. Transfer of a radioactive phosphate from DNA to the 50-kDa protein and gel retardation experiments indicate that this protein forms the covalent complex with DNA. A blot overlay assay identifies the 150-kDa protein as the potential ATPase. This is the first evidence that a reverse gyrase can be a topoisomerase consisting of two protomers. In analogy with the DNA gyrase A subunit (DNA breakage and reunion activity) and the B subunit (ATPase), the 50- and 150-kDa components of Mka reverse gyrase have been designated the A and B subunits, respectively. Methanopyrus reverse gyrase changes DNA linking number in steps of one and its A subunit covalently binds to the 5'-DNA phosphoryl group. It nicks DNA at sites that predominantly have a cytosine at the -4-position. The same rule was derived previously for monomeric reverse gyrase from sulfur-metabolizing hyperthermophiles and for topoisomerase I from mesophilic bacteria. Based on these results, Mka reverse gyrase is classified as belonging to group A of type I topoisomerases. The structural diversity of type I group A topoisomerases parallels the diversity of type II enzymes and suggests the evolution of an essential function by gene fusion.

Evolution of RNA editing in kinetoplastid protozoa

The editing of RNA in trypanosomatid mitochondria involves the insertion and occasional deletion of uridine residues within coding regions of maxicircle messenger RNA transcripts. The extent to which the transcripts of homologous genes undergo editing differs in different species. In some, entire genes are edited (pan-editing), whereas in others, editing is limited to the 5' termini of editing domains (5' editing). Here we investigate which type of editing is ancestral and which is derived, by analysing RNA editing in the different lineages, using a kinetoplastid phylogeny reconstructed from nuclear small subunit ribosomal RNA sequences. We conclude that the ancestral cryptogenes were pan-edited, and we hypothesize that the 5'-edited homologues were generated by several independent events from partially edited RNAs, in which case editing may be a more primitive mechanism than previously thought.

Reconstructing evolutionary trees from DNA and protein sequences: paralinear distances

The reconstruction of phylogenetic trees from DNA and protein sequences is confounded by unequal rate effects. These effects can group rapidly evolving taxa with other rapidly evolving taxa, whether or not they are genealogically related. All algorithms are sensitive to these effects whenever the assumptions on which they are based are not met. The algorithm presented here, called paralinear distances, is valid for a much broader class of substitution processes than previous algorithms and is accordingly less affected by unequal rate effects. It may be used with all nucleic acid, protein, or other sequences, provided that their evolution may be modeled as a succession of Markov processes. The properties of the method have been proven both analytically and by computer simulations. Like all other methods, paralinear distances can fail when sequences are misaligned or when site-to-site sequence variation of rates is extensive. To examine the usefulness of paralinear distances, the "origin of the eukaryotes" has been investigated by the analysis of elongation factor Tu sequences with a variety of sequence alignments. It has been found that the order in which sequences are pairwise aligned strongly determines the topology which is reconstructed by paralinear distances (as it does for all other reconstruction methods tested). When the parts of the alignment that are unaffected by alignment order are analyzed, paralinear distances strongly select the eocyte topology. This provides evidence that the eocyte prokaryotes are the closest prokaryotic relatives of the eukaryotes.

Purification and characterization of DNA topoisomerase V. An enzyme from the hyperthermophilic prokaryote Methanopyrus kandleri that resembles eukaryotic topoisomerase I

DNA topoisomerase V is a novel prokaryotic enzyme related to eukaryotic topoisomerase I. The enzyme is a type I DNA topoisomerase and is recognized by polyclonal antibody against human topoisomerase I. We describe its purification from the hyperthermophilic methanogen Methanopyrus kandleri. The enzyme has high activity in crude extracts and is present in at least 1,500 copies/cell. Topoisomerase V migrates as a 110-kDa polypeptide in SDS-polyacrylamide gel electrophoresis and as a 142-kDa globular protein in gel filtration. It is active up to at least 100 degrees C on both positively and negatively supercoiled DNA and is not inhibited by single-stranded DNA. The enzyme works from 1 to 650 mM NaCl and up to 3.1 M potassium glutamate. It acts processively at low ionic strength and distributively at high NaCl or KCl concentration. Magnesium is not required and does not stimulate the enzymatic activity. Under DNA denaturing conditions, topoisomerase V catalyzes an unlinking reaction which results in substantial reduction in the linking number of closed circular DNA. The driving force for this process is DNA melting. Camptothecin is not nearly as good an inhibitor for topoisomerase V as it is for eukaryotic topoisomerase I. The unique occurrence of two major type I topoisomerases (reverse gyrase and topoisomerase V) in M. kandleri may shed new light on the evolution of this family of enzymes and supports the concept of a distant but significant relationship between some hyperthermophilic organisms and eukaryotes.

DNA topoisomerase V is a relative of eukaryotic topoisomerase I from a hyperthermophilic prokaryote

The DNA topoisomerases are ubiquitous enzymes that fulfil vital roles in the replication, transcription and recombination of DNA by carrying out DNA-strand passage reactions. Here we characterize a prokaryotic counterpart to the eukaryotic topoisomerase I in the hyperthermophilic methanogen Methanopyrus kandleri. The new enzyme, called topoisomerase V, has the following properties in common with eukaryotic topoisomerase I, which distinguish it from all other known prokaryotic topoisomerases: (1) its activity is Mg(2+)-independent; (2) it relaxes both negatively and positively supercoiled DNA; (3) it makes a covalent complex with the 3' end of the broken DNA strand; and (4) it is recognized by antibody raised against human topoisomerase I. Eukaryotic-like enzymes have been discovered in some hyperthermophilic prokaryotes, namely the eocytes and the extremely thermophilic archaebacteria, and hyperthermophilic homologues of eukaryotic DNA polymerase-alpha, transcription factor IIB and DNA ligase have all been reported. Thus our findings support the idea that some essential parts of the eukaryotic transcription-translation and replication machineries were in place before the emergence of eukaryotes, and that the closest living relatives of eukaryotes may be hyperthermophiles.

Rough genes with Deformed homeobox substitutions exhibit rough regulatory specificity during Drosophila eye development

In certain cases, homeobox genes with different in vivo roles encode proteins with similar in vitro DNA binding specificities. To test the role of the homeobox in the regulatory specificity of such genes, rough homeobox sequences were changed in part or entirely to those of the Deformed gene, and the modified rough genes tested for their ability to rescue the rough mutant phenotype. Surprisingly, the chimaeric genes retained levels of rough regulatory specificity but acquired no novel functions. These results suggest that factors other than the DNA binding specificity of the homeodomain play crucial roles in determining the target, and thus the regulatory specificity, of such proteins.

Evidence that eukaryotes and eocyte prokaryotes are immediate relatives

The phylogenetic origin of eukaryotes has been unclear because eukaryotic nuclear genes have diverged substantially from prokaryotic ones. The genes coding for elongation factor EF-1 alpha were compared among various organisms. The EF-1 alpha sequences of eukaryotes contained an 11-amino acid segment that was also found in eocytes (extremely thermophilic, sulfur-metabolizing bacteria) but that was absent in all other bacteria. The related (paralogous) genes encoding elongation factor EF-2 and initiation factor IF-2 also lacked the 11-amino acid insert. These data imply that the eocytes are the closest surviving relatives (sister taxon) of the eukaryotes.

Mapping the three-dimensional locations of ribosomal RNA and proteins

Seven regions of 16S rRNA have been located on the surface of the 30S ribosomal subunit by DNA hybridization electron microscopy in our laboratory. In addition, we have recently mapped the three-dimensional locations of an additional seven small ribosomal proteins by immunoelectron microscopy. The information from the direct mapping of the sites on rRNA has been incorporated into a model for the tertiary structure of 16S rRNA, accounting for approximately 40% of the total 16S rRNA. A novel structure, the platform ring, is proposed for a region of rRNA within the central domain. This structure rings the edges of the platform and includes regions 655-751 and 769-810. Another region, the recognition complex, consists of nucleotides 500-545, and occupies a region on the exterior surface of the subunit, near the EF-Tu binding site. In addition, 19 of the 21 small subunit ribosomal proteins have been mapped by immunoelectron microscopy in our laboratory. In order to evaluate the reliability of our model for the three-dimensional distribution of 16S rRNA, we have predicted which sites of rRNA are adjacent to ribosomal proteins and compared these predictions with r-protein protection studies of others. Good correlation between the model, the locations of rRNA sites, the locations of ribosomal proteins, and regions of rRNA protected by ribosomal proteins, provides independent support for this model.

Tracing origins with molecular sequences: metazoan and eukaryotic beginnings

Milestones in the evolution of the eukaryotic cell are being discovered through the analysis of molecular sequences. As sequence data become increasingly plentiful, our ability to reconstruct the most distant evolutionary branchings of evolutionary trees is limited only by the mathematics of phylogenetic reconstruction. Analysis of ribosomal RNAs agrees with traditional analyses of morphological and developmental characters that all multicellular animals probably arose from a common ancestor, but highlights one of the major limitations of the various mathematical algorithms used. Refined methods of sequence analysis also suggest a previously unsuspected sister relationship between the eukaryotic nucleus and eocytes, a group of extremely thermophilic, sulfur-metabolizing bacteria, that questions the classical eukaryote/prokaryote division.

DNA-hybridization electron microscopy. Localization of five regions of 16 S rRNA on the surface of 30 S ribosomal subunits

DNA-hybridization electron microscopy has been used to locate five regions of 16 S rRNA on the surface of 30 S ribosomal subunits. Biotinylated DNA probes that are complementary to selected regions of 16 S rRNA were hybridized to activated 30 S ribosomal subunits. These hybridized probes were reacted with avidin and localized by electron microscopy. The specificity of DNA binding was monitored with RNase H, which recognizes RNA-DNA hybrids and cleaves the RNA. Three of the five sequences examined were mapped on the platform. These sequences are 686-703, 714-733 and 787-803. Region 1492-1505 is mapped in the cleft and region 518-533 is at the neck on the side opposite the platform, respectively.

DNA-hybridization electron microscopy tertiary structure of 16 S rRNA

Seven regions of 16 S rRNA have been located on the surface of the 30 S ribosomal subunit by DNA-hybridization electron microscopy. This information has been incorporated into a model for the tertiary structure of 16 S rRNA, accounting for approximately 40% of the total 16 S rRNA. A structure labeled the platform ring is proposed for a region of rRNA within the central domain. This structure rings the edges of the platform and includes regions 655-751 and 769-810. Another region, the recognition complex, consists of nucleotides 500 to 545, and occupies a region on the exterior surface of the subunit near the elongation factor Tu binding site. Ribosomal proteins that have been mapped by immunoelectron microscopy are superimposed onto the model in order to examine possible regions of interaction. Good correlation between the model locations of ribosomal proteins, and regions of rRNA protected by ribosomal proteins provide independent support for this model.

Origin of the Metazoa

The origin of the multicellular animals has been investigated by rate invariant analysis of 18S rRNA sequences. These analyses indicate that (i) the Metazoa is a monophyletic taxon; (ii) the Deuterostomia is a monophyletic taxon; (iii) the Annelida-Mollusca lineage is the sister group of an arthropod subgroup; and (iv) the last common ancestor of the Annelida-Mollusca lineage is most parsimoniously derived from a segmented, hemocoelic ancestor with an open circulatory system.

Origin of the eukaryotic nucleus: eukaryotes and eocytes are genotypically related

The origin of the eukaryotic nucleus is difficult to reconstruct. While eukaryotic organelles (chloroplast, mitochondrion) are eubacterial endosymbionts, the source of nuclear genes has been obscured by multiple nucleotide substitutions. Using evolutionary parsimony, a newly developed rate-invariant treeing algorithm, the eukaryotic rRNA genes are shown to have evolved from the eocytes, a group of extremely thermophilic, sulfur-metabolizing, anucleate cells. The deepest bifurcation yet found separates the reconstructed tree into two taxonomic divisions. These are a proto-eukaryotic group (karyotes) and an essentially bacterial one (parkaryotes). Within the precision of the rooting procedure, the tree is not consistent with either the prokaryotic--eukaryotic or the archaebacterial--eubacterial--eukaryotic groupings. It implies that the last common ancestor of extant life, and the early ancestors of eukaryotes, very likely lacked nuclei, metabolized sulfur, and lived at near boiling temperatures.

Evolution of parasitism: kinetoplastid protozoan history reconstructed from mitochondrial rRNA gene sequences

A phylogenetic tree for the evolution of five representative species from four genera of kinetoplastid protozoa was constructed from comparison of the mitochondrial 9S and 12S rRNA gene sequences and application of both parsimony and evolutionary parsimony algorithms. In the rooted version of the tree, the monogenetic species Crithidia fasciculata is the most deeply rooted, followed by another monogenetic species, Leptomonas sp. The three digenetic species Trypanosoma cruzi, Trypanosoma brucei, and Leishmania tarentolae branch from the Leptomonas line. The substitution rates for the T. brucei and T. cruzi sequences were 3-4 times greater than that of the L. tarentolae sequences. This phylogenetic tree is consistent with our cladistic analysis of the biological evidence including life cycles for these five species. A tentative time scale can be assigned to the nodes of this tree by assuming that the common ancestor of the digenetic parasites predated the separation of South America and Africa and postdated the first fossil appearance of its host (inferred by parsimony analysis). This time scale predicts that the deepest node occurred at 264 +/- 51 million years ago, at a time commensurate with the fossil origins of the Hemiptera insect host. This implies that the ancestral kinetoplastid and its insect host appeared at approximately the same time. The molecular data suggest that these eukaryotic parasites have an evolutionary history that extends back to the origin of their insect host.

Origin of the eukaryotic nucleus determined by rate-invariant analysis of rRNA sequences

The origin of the eukaryotic nucleus is difficult to reconstruct. Eukaryotic organelles (chloroplast, mitochondrion) are eubacterial endosymbionts, but the source of nuclear genes has been obscured by multiple nucleotide substitutions. Using evolutionary parsimony, a newly developed rate-invariant treeing algorithm, the eukaryotic ribosomal rRNA genes are shown to have evolved from the eocytes, a group of extremely thermophilic, sulphur-metabolizing, anucleate cells. The deepest bifurcation yet found separates the reconstructed tree into two taxonomic divisions. These are a proto-eukaryotic group (karyotes) and an essentially bacterial one (parkaryotes). Within the precision of the rooting procedure, the tree is not consistent with either the prokaryotic-eukaryotic or the archaebacterial-eubacterial-eukaryotic groupings. It implies that the last common ancestor of extant life, and the early ancestors of eukaryotes, probably lacked nuclei, metabolized sulphur and lived at near-boiling temperatures.

Determining evolutionary distances from highly diverged nucleic acid sequences: operator metrics

Operator metrics are explicitly designed to measure evolutionary distances from nucleic acid sequences when substitution rates differ greatly among the organisms being compared, or when substitutions have been extensive. Unlike lengths calculated by the distance matrix and parsimony methods, in which substitutions in one branch of a tree can alter the measured length of another branch, lengths determined by operator metrics are not affected by substitutions outside the branch. In the method, lengths (operator metrics) corresponding to each of the branches of an unrooted tree are calculated. The metric length of a branch reconstructs the number of (transversion) differences between sequences at a tip and a node (or between nodes) of a tree. The theory is general and is fundamentally independent of differences in substitution rates among the organisms being compared. Mathematically, the independence has been obtained because the metrics are eigenvectors of fundamental equations which describe the evolution of all unrooted trees. Even under conditions when both the distance matrix method or a simple parismony length method are shown to indicate lengths that are an order of magnitude too large or too small, the operator metrics are accurate. Examples, using data calculated with evolutionary rates and branchings designed to confuse the measurement of branch lengths and to camouflage the topology of the true tree, demonstrate the validity of operator metrics. The method is robust. Operator metric distances are easy to calculate, can be extended to any number of taxa, and provide a statistical estimate of their variances. The utility of the method is demonstrated by using it to analyze the origins and evolution of chloroplasts, mitochondria, and eubacteria.

A rate-independent technique for analysis of nucleic acid sequences: evolutionary parsimony

The method of evolutionary parsimony--or operator invariants--is a technique of nucleic acid sequence analysis related to parsimony analysis and explicitly designed for determining evolutionary relationships among four distantly related taxa. The method is independent of substitution rates because it is derived from consideration of the group properties of substitution operators rather than from an analysis of the probabilities of substitution in branches of a tree. In both parsimony and evolutionary parsimony, three patterns of nucleotide substitution are associated one-to-one with the three topologically linked trees for four taxa. In evolutionary parsimony, the three quantities are operator invariants. These invariants are the remnants of substitutions that have occurred in the interior branch of the tree and are analogous to the substitutions assigned to the central branch by parsimony. The two invariants associated with the incorrect trees must equal zero (statistically), whereas only the correct tree can have a nonzero invariant. The chi 2-test is used to ascertain the nonzero invariant and the statistically favored tree. Examples, obtained using data calculated with evolutionary rates and branchings designed to camouflage the true tree, show that the method accurately predicts the tree, even when substitution rates differ greatly in neighboring peripheral branches (conditions under which parsimony will consistently fail). As the number of substitutions in peripheral branches becomes fewer, the parsimony and the evolutionary-parsimony solutions converge. The method is robust and easy to use.

Elongation factor Tu localized on the exterior surface of the small ribosomal subunit

Protein synthesis elongation factor Tu has been mapped on the surface of the ribosome by immunoelectron microscopy. The EF-Tu binding site is located near the concave portion of the small subunit at the level of the neck and extends 60 to 70 A above and below the neck. The binding site is present on the side of the subunit opposite the platform, i.e. on the exterior subunit surface. This EF-Tu site differs from that found for elongation factor G in that the EF-Tu site is significantly more exposed to the cytoplasm.

Monomers, dimers, and minifilaments of vertebrate skeletal myosin in the presence of sodium pyrophosphate

The self-assembly of myosin in the presence of sodium pyrophosphate was studied in the pH range between 7.0 and 8.5. As evidenced by sedimentation velocity (S0(20,w) = 6.30 S) and light-scattering measurements (molecular weight of 470 000; radius of gyration = 45 nm), myosin existed in a predominantly monomeric form in the presence of 5 mM sodium pyrophosphate at pH 8.5 and above. The concentration-dependent monomer-dimer equilibrium could be easily shifted toward dimeric species at pH 8.0 in the presence of 5 mM sodium pyrophosphate and 5 mM 2-[bis(2-hydroxyethyl)amino]-2-(hydroxymethyl)-1,3-propanediol. The estimated parameters of the dimeric particles were S0(20,w) between 10 and 11 S, molecular weight of 1.1 X 10(6), and radius of gyration = 52 nm. These results are consistent with a head to tail (parallel) arrangement of staggered myosin molecules in the dimer. At lower pH values (7.5), and in the presence of 10 mM sodium pyrophosphate, the monomer-dimer species were in dynamic equilibrium with myosin minifilaments. At pH 7.0, the minifilaments appeared to be the only detectable species present in solutions of myosin in 5 mM sodium pyrophosphate. The molecular parameters of these minifilaments, including sedimentation and viscosity coefficients, molecular weight, radius of gyration, and morphological appearance, were almost indistinguishable from those obtained for myosin minifilaments prepared in 10 mM citrate-tris(hydroxymethyl)aminomethane at pH 8.0 [Reisler, E., Smith, C., & Seegan, G. (1980) J. Mol. Biol. 143, 129-145]. The equilibrium polymerization reactions of myosin in sodium pyrophosphate are discussed in the context of minifilament assembly.