Phylogeny of mRNA capping enzymes

The m7GpppN cap structure of eukaryotic mRNA is formed cotranscriptionally by the sequential action of three enzymes: RNA triphosphatase, RNA guanylyltransferase, and RNA (guanine-7)-methyltransferase. A multifunctional polypeptide containing all three active sites is encoded by vaccinia virus. In contrast, fungi and Chlorella virus encode monofunctional guanylyltransferase polypeptides that lack triphosphatase and methyltransferase activities. Transguanylylation is a two-stage reaction involving a covalent enzyme-GMP intermediate. The active site is composed of six protein motifs that are conserved in order and spacing among yeast and DNA virus capping enzymes. We performed a structure–function analysis of the six motifs by targeted mutagenesis of Ceg1, the Saccharomyces cerevisiae guanylyltransferase. Essential acidic, basic, and aromatic functional groups were identified. The structural basis for covalent catalysis was illuminated by comparing the mutational results with the crystal structure of the Chlorella virus capping enzyme. The results also allowed us to identify the capping enzyme of Caenorhabditis elegans. The 573-amino acid nematode protein consists of a C-terminal guanylyltransferase domain, which is homologous to Ceg1 and is strictly conserved with respect to all 16 amino acids that are essential for Ceg1 function, and an N-terminal phosphatase domain that bears no resemblance to the vaccinia triphosphatase domain but, instead, has strong similarity to the superfamily of protein phosphatases that act via a covalent phosphocysteine intermediate.

Phylogeny of rice genomes with emphasis on origins of allotetraploid species

The rice genus, Oryza, which comprises 23 species and 9 recognized genome types, represents an enormous gene pool for genetic improvement of rice cultivars. Clarification of phylogenetic relationships of rice genomes is critical for effective utilization of the wild rice germ plasm. By generating and comparing two nuclear gene (Adh1 and Adh2) trees and a chloroplast gene (matK) tree of all rice species, phylogenetic relationships among the rice genomes were inferred. Origins of the allotetraploid species, which constitute more than one-third of rice species diversity, were reconstructed based on the Adh gene phylogenies. Genome types of the maternal parents of allotetraploid species were determined based on the matK gene tree. The phylogenetic reconstruction largely supports the previous recognition of rice genomes. It further revealed that the EE genome species is most closely related to the DD genome progenitor that gave rise to the CCDD genome. Three species of the CCDD genome may have originated through a single hybridization event, and their maternal parent had the CC genome. The BBCC genome species had different origins, and their maternal parents had either a BB or CC genome. An additional genome type, HHKK, was recognized for Oryza schlechteri and Porteresia coarctata, suggesting that P. coarctata is an Oryza species. The AA genome lineage, which contains cultivated rice, is a recently diverged and rapidly radiated lineage within the rice genus.

The phylogeny of closely related species as revealed by the genealogy of a speciation gene, Odysseus

Molecular differentiation between races or closely related species is often incongruent with the reproductive divergence of the taxa of interest. Shared ancient polymorphism and/or introgression during secondary contact may be responsible for the incongruence. At loci contributing to speciation, these two complications should be minimized (1, 2); hence, their variation may more faithfully reflect the history of the species' reproductive differentiation. In this study, we analyzed DNA polymorphism at the Odysseus (OdsH) locus of hybrid sterility between Drosophila mauritiana and Drosophila simulans and were able to verify such a prediction. Interestingly, DNA variation only a short distance away (1.8 kb) appears not to be influenced by the forces that shape the recent evolution of the OdsH coding region. This locus thus may represent a test case of inferring phylogeny of very closely related species.

Gene content phylogeny of herpesviruses

Clusters of orthologous groups [COGs; Tatusov, R. L., Koonin, E. V. & Lipman, D. J. (1997) Science 278, 631–637] were identified for a set of 13 completely sequenced herpesviruses. Each COG represented a family of gene products conserved across several herpes genomes. These families were defined without using an arbitrary threshold criterion based on sequence similarity. The COG technique was modified so that variable stringency in COG construction was possible. High stringencies identify a core set of highly conserved genes. Varying COG stringency reveals differences in the degree of conservation between functional classes of genes. The COG data were used to construct whole-genome phylogenetic trees based on gene content. These trees agree well with trees based on other methods and are robust when tested by bootstrap analysis. The COG data also were used to construct a reciprocal tree that clustered genes with similar phylogenetic profiles. This clustering may give clues to genes with related functions or with related histories of acquisition and loss during herpesvirus evolution.

Diversity and phylogeny of gephyrin: Tissue-specific splice variants, gene structure, and sequence similarities to molybdenum cofactor-synthesizing and cytoskeleton-associated proteins

Gephyrin is essential for both the postsynaptic localization of inhibitory neurotransmitter receptors in the central nervous system and the biosynthesis of the molybdenum cofactor (Moco) in different peripheral organs. Several alternatively spliced gephyrin transcripts have been identified in rat brain that differ in their 5′ coding regions. Here, we describe gephyrin splice variants that are differentially expressed in non-neuronal tissues and different regions of the adult mouse brain. Analysis of the murine gephyrin gene indicates a highly mosaic organization, with eight of its 29 exons corresponding to the alternatively spliced regions identified by cDNA sequencing. The N- and C-terminal domains of gephyrin encoded by exons 3–7 and 16–29, respectively, display sequence similarities to bacterial, invertebrate, and plant proteins involved in Moco biosynthesis, whereas the central exons 8, 13, and 14 encode motifs that may mediate oligomerization and tubulin binding. Our data are consistent with gephyrin having evolved from a Moco biosynthetic protein by insertion of protein interaction sequences.

Phylogeny of genes for secretion NTPases: Identification of the widespread tadA subfamily and development of a diagnostic key for gene classification

Macromolecular transport systems in bacteria currently are classified by function and sequence comparisons into five basic types. In this classification system, type II and type IV secretion systems both possess members of a superfamily of genes for putative NTP hydrolase (NTPase) proteins that are strikingly similar in structure, function, and sequence. These include VirB11, TrbB, TraG, GspE, PilB, PilT, and ComG1. The predicted protein product of tadA, a recently discovered gene required for tenacious adherence of Actinobacillus actinomycetemcomitans, also has significant sequence similarity to members of this superfamily and to several unclassified and uncharacterized gene products of both Archaea and Bacteria. To understand the relationship of tadA and tadA-like genes to those encoding the putative NTPases of type II/IV secretion, we used a phylogenetic approach to obtain a genealogy of 148 NTPase genes and reconstruct a scenario of gene superfamily evolution. In this phylogeny, clear distinctions can be made between type II and type IV families and their constituent subfamilies. In addition, the subgroup containing tadA constitutes a novel and extremely widespread subfamily of the family encompassing all putative NTPases of type IV secretion systems. We report diagnostic amino acid residue positions for each major monophyletic family and subfamily in the phylogenetic tree, and we propose an easy method for precisely classifying and naming putative NTPase genes based on phylogeny. This molecular key-based method can be applied to other gene superfamilies and represents a valuable tool for genome analysis.

Molecular phylogeny of the extinct ground sloth Mylodon darwinii.

DNA was extracted from the remains of 35 ground sloths from various parts of North and South America. Two specimens of Mylodon darwinii, a species that went extinct at the end of the last glaciation, yielded amplifiable DNA. However, of the total DNA extracted, only approximately 1/1000 originated from the sloth, whereas a substantial part of the remainder was of bacterial and fungal origin. In spite of this, > 1100 bp of sloth mitochondrial rDNA sequences could be reconstructed from short amplification products. Phylogenetic analyses using homologous sequences from all extant edentate groups suggest that Mylodon darwinii was more closely related to the two-toed than the three-toed sloths and, thus, that an arboreal life-style has evolved at least twice among sloths. The divergence of Mylodon and the two-toed sloth furthermore allows a date for the radiation of armadillos, anteaters, and sloths to be estimated. This result shows that the edentates differ from other mammalian orders in that they contain lineages that diverged before the end of the Cretaceous Period.

Different population dynamics of human T cell lymphotropic virus type II in intravenous drug users compared with endemically infected tribes

The phylogeny of human T cell lymphotropic virus type II (HTLV-II) was investigated by using strains isolated from Amerindian and Pygmy tribes, in which the virus is maintained primarily through mother-to-child transmission via breast-feeding, and strains from intravenous drug users (IDUs), in which spread is mainly blood-borne via needle sharing. Molecular clock analysis showed that HTLV-II has two different evolutionary rates with the molecular clock for the virus in IDUs ticking 150–350 times faster than the one in endemically infected tribes: 2.7 × 10−4 compared with 1.71/7.31 × 10−7 nucleotide substitutions per site per year in the long terminal repeat region. This dramatic acceleration of the evolutionary rate seems to be related with the mode of transmission. Mathematical models showed the correlation of these two molecular clocks with an endemic spread of HTLV-II in infected tribes compared with the epidemic spread in IDUs. We also noted a sharp increase in the population size of the virus among IDUs during the last decades probably caused by the worldwide increase in intravenous drug use.

Recurrent invasion and extinction of a selfish gene

Homing endonuclease genes show super-Mendelian inheritance, which allows them to spread in populations even when they are of no benefit to the host organism. To test the idea that regular horizontal transmission is necessary for the long-term persistence of these genes, we surveyed 20 species of yeasts for the ω-homing endonuclease gene and associated group I intron. The status of ω could be categorized into three states (functional, nonfunctional, or absent), and status was not clustered on the host phylogeny. Moreover, the phylogeny of ω differed significantly from that of the host, strong evidence of horizontal transmission. Further analyses indicate that horizontal transmission is more common than transposition, and that it occurs preferentially between closely related species. Parsimony analysis and coalescent theory suggest that there have been 15 horizontal transmission events in the ancestry of our yeast species, through simulations indicate that this value is probably an underestimate. Overall, the data support a cyclical model of invasion, degeneration, and loss, followed by reinvasion, and each of these transitions is estimated to occur about once every 2 million years. The data are thus consistent with the idea that frequent horizontal transmission is necessary for the long-term persistence of homing endonuclease genes, and further, that this requirement limits these genes to organisms with easily accessible germ lines. The data also show that mitochondrial DNA sequences are transferred intact between yeast species; if other genes do not show such high levels of horizontal transmission, it would be due to lack of selection, rather than lack of opportunity.

The heme-copper oxidases of Thermus thermophilus catalyze the reduction of nitric oxide: Evolutionary implications

We show that the heme-copper terminal oxidases of Thermus thermophilus (called ba3 and caa3) are able to catalyze the reduction of nitric oxide (NO) to nitrous oxide (N2O) under reducing anaerobic conditions. The rate of NO consumption and N2O production were found to be linearly dependent on enzyme concentration, and activity was abolished by enzyme denaturation. Thus, contrary to the eukaryotic enzyme, both T. thermophilus oxidases display a NO reductase activity (3.0 ± 0.7 mol NO/mol ba3 × min and 32 ± 8 mol NO/mol caa3 × min at [NO] ≈ 50 μM and 20°C) that, though considerably lower than that of bona fide NO reductases (300–4,500 mol NO/mol enzyme × min), is definitely significant. We also show that for ba3 oxidase, NO reduction is associated to oxidation of cytochrome b at a rate compatible with turnover, suggesting a mechanism consistent with the stoichiometry of the overall reaction. We propose that the NO reductase activity of T. thermophilus oxidases may depend on a peculiar CuB+ coordination, which may be revealed by the forthcoming three-dimensional structure. These findings support the hypothesis of a common phylogeny of aerobic respiration and bacterial denitrification, which was proposed on the basis of structural similarities between the Pseudomonas stutzeri NO reductase and the cbb3 terminal oxidases. Our findings represent functional evidence in support of this hypothesis.

Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.

Relationships of cereal crops and other grasses

The grass family includes some 10,000 species, and it encompasses tremendous morphological, physiological, ecological, and genetic diversity. The phylogeny of the family is becoming increasingly well understood. There were two major radiations of grasses, an early diversification leading to the subfamilies Pooideae, Bambusoideae, and Oryzoideae, and a later one leading to Panicoideae, Chloridoideae, Centothecoideae, and Arundinoideae. The phylogeny can be used to determine the direction of changes in genome arrangement and genome size.

Retroposon analysis of major cetacean lineages: The monophyly of toothed whales and the paraphyly of river dolphins

SINE (short interspersed element) insertion analysis elucidates contentious aspects in the phylogeny of toothed whales and dolphins (Odontoceti), especially river dolphins. Here, we characterize 25 informative SINEs inserted into unique genomic loci during evolution of odontocetes to construct a cladogram, and determine a total of 2.8 kb per taxon of the flanking sequences of these SINE loci to estimate divergence times among lineages. We demonstrate that: (i) Odontocetes are monophyletic; (ii) Ganges River dolphins, beaked whales, and ocean dolphins diverged (in this order) after sperm whales; (iii) three other river dolphin taxa, namely the Amazon, La Plata, and Yangtze river dolphins, form a monophyletic group with Yangtze River dolphins being the most basal; and (iv) the rapid radiation of extant cetacean lineages occurred some 28–33 million years B.P., in strong accord with the fossil record. The combination of SINE and flanking sequence analysis suggests a topology and set of divergence times for odontocete relationships, offering alternative explanations for several long-standing problems in cetacean evolution.

Evolution and divergence of sodium channel genes in vertebrates

Invertebrate species possess one or two Na+ channel genes, yet there are 10 in mammals. When did this explosive growth come about during vertebrate evolution? All mammalian Na+ channel genes reside on four chromosomes. It has been suggested that this came about by multiple duplications of an ancestral chromosome with a single Na+ channel gene followed by tandem duplications of Na+ channel genes on some of these chromosomes. Because a large-scale expansion of the vertebrate genome likely occurred before the divergence of teleosts and tetrapods, we tested this hypothesis by cloning Na+ channel genes in a teleost fish. Using an approach designed to clone all of the Na+ channel genes in a genome, we found six Na+ channel genes. Phylogenetic comparisons show that each teleost gene is orthologous to a Na+ channel gene or gene cluster on a different mammalian chromosome, supporting the hypothesis that four Na+ channel genes were present in the ancestors of teleosts and tetrapods. Further duplications occurred independently in the teleost and tetrapod lineages, with a greater number of duplications in tetrapods. This pattern has implications for the evolution of function and specialization of Na+ channel genes in vertebrates. Sodium channel genes also are linked to homeobox (Hox) gene clusters in mammals. Using our phylogeny of Na+ channel genes to independently test between two models of Hox gene evolution, we support the hypothesis that Hox gene clusters evolved as (AB) (CD) rather than {D[A(BC)]}.

Population dynamics of flaviviruses revealed by molecular phylogenies.

The phylogeny of 123 complete envelope gene sequences was reconstructed in order to understand the evolution of tick- and mosquito-borne flaviviruses. An analysis of phylogenetic tree structure reveals a continual and asymmetric branching process in the tick-borne flaviviruses, compared with an explosive radiation in the last 200 years in viruses transmitted by mosquitoes. The distinction between these two viral groups probably reflects differences in modes of dispersal, propagation, and changes in the size of host populations. The most serious implication of this work is that growing human populations are being exposed to an expanding range of increasingly diverse viral strains.