Comparison of RNA sequences from widely divergent species

Summary Tests that take into account the effects of gaps have been applied to 5S ribosomal RNA sequences from the bacteria,E. coli andP. fluorescens, and from KB carcinoma cells. The 5S RNAs from KB andP. fluorescens, when compared to that ofE. coli are shown to be more similar than random sequences of the same composition. Intrasequential analyses of 5S RNAs give some evidence for partial gene duplication or repetitive subsequences, but the proposed duplication of Brownlee, Sanger and Barrell (1968) is not supported by our data.     

Molecular evolution of theSaccharomyces cerevisiae histone gene loci

Summary The core histone genes ofSaccharomyces cerevisiae are arranged as duplicate nonallelic sets of specifically paired genes. The identity of structural organization between the duplicated gene pairs would have its simplest evolutionary origin in the duplication of a complete locus in a single event. In such a case, the time since the duplication of one of the genes should be identical to that since duplication of the gene adjacent to it on the chromosome. A calculation of the evolutionary distances between the coding DNA sequences of the histone genes leads to a duplication paradox: The extents of sequence divergence in the silent component of third-base positions for adjacent pairs of genes are not identical. Estimates of the evolutionary distance between the two H3-H4 noncoding intergene DNA sequences are large; the divergence between the two separate sequences is indistinguishable from the divergence between either of the regions and a randomly generated permutation of itself. These results suggest that the duplication event may have occurred much earlier than previously estimated. The potential age of the duplication, and the attractive simplicity of the duplication of both the H3-H4 and the H2A-H2B gene pairs having taken place in a single event, leads to the hypothesis that modern haploidS. cerevisiae may have evolved by diploidization or fusion of two ancient fungi.     

Tandem duplication of alkaline phosphatase genes and polymorphism in the intergenic sequence in<Emphasis Type="Italic"> Bombyx mori</Emphasis>

Alkaline phosphatases are ubiquitous in organisms from bacteria to human. Two alkaline phosphatase genes, Alp-m and Alp-s, were independently cloned from the silkworm Bombyx mori. They were mapped to a small DNA region and shown to be organized in tandem. Exon-intron structures of the two genes were highly conserved, with the exception of the second intron in Alp-m, which has no counterpart in Alp-s. The similarity between the nucleotide sequences of the exons of the two genes was strikingly high (60–79%), suggesting that Alp-m and Alp-s originated from a duplication of their common ancestor gene. The intergenic sequence between the two Alp genes shows length polymorphism in different B. mori strains, which can be explained by presence/absence of two putative insertion sequences. This structural variation suggests a possible scenario for the divergence of the two Alp genes after the duplication event.Communicated by G. Reuter     

Amino acid sequences of stomach and nonstomach lysozymes of ruminants

Summary Complete amino acid sequences are presented for lysozymesc from camel and goat stomachs and compared to sequences of other lysozymesc. Tree analysis suggests that the rate of amino acid replacement went up as soon as lysozyme was recruited for the stomach function in early ruminants. The two lysozymes from goat stomach are the products of a gene duplication that probably took place before the divergence of cow, goat, and deer about 25 million years ago. Partial sequences of three lysozymes from goat tears indicated that (a) the goat tear family of lysozymes may have diverged from the stomach lysozyme family by an ancient duplication and (b) later duplications are probably responsible for the multiple forms of tear and milk lysozymes in ruminants.     

Evolution of the proteasome components

 A phylogenetic analysis of proteasome subunits revealed two major families (α and β) which originated by an ancient gene duplication prior to the divergence of archaebacteria and eukaryotes. Numerous gene duplications have subsequently occurred in eukaryotes; at least nine of these duplications were shown to have occurred prior to the divergence of animals and fungi. In mammals, two genes encoding proteasome subunits (LMP2 and LMP7) are located in the major histocompatibility complex (MHC) region and play a specific role in generation of peptides for presentation by class I MHC molecules. Phylogenetic analysis of LMP7 and related sequences from mammals and lower vertebrates indicated that this locus arose by gene duplication prior to the divergence of jawed and jawless vertebrates; the time of this duplication was estimated to have been about 600 million years ago. The evolutionary history of the proteasome subunits provides support for a model of the evolution of new gene function postulating that, after gene duplication, the proteins encoded by daughter loci can adapt to specialized functions previously performed by the product of a single generalized ancestral locus. Received: 19 August 1996 / Revised: 24 December 1996     

Exploring the Evolutionary History of the Alcohol Dehydrogenase Gene (<Emphasis Type="Italic">Adh</Emphasis>) Duplication in Species of the Family Tephritidae

In the olive fruit fly Bactrocera oleae and the med fly Ceratitis capitata previous studies have shown the existence of two Adh genes in each species. This observation, in combination with the former finding that various Drosophila species of virilis and repleta group encode two isozymes of ADH which are the result of a gene duplication, challenged us to address a scenario dealing with the evolutionary history of the Adh gene duplication in Tephritidae. In our lab we proceeded to the cloning and sequence analysis of Adh genes from more tephritid species, a prerequisite for further study of this issue. Here we show that phylogenetic trees produced from either the nucleotide or the amino acid sequences of 14 tephritid Adh genes consisted of two main clusters, with Adh sequences of the same type grouping together (i.e., Adh1 sequences form a cluster and Adh2 sequences form a second one), as expected if there was one duplication event before speciation within the family Tephritidae. We used the amount of divergence between the two isozymic forms of Adh of the species carrying both Adh1 and Adh2 genes to obtain an estimate of the age of the duplication event. Interestingly, our data again support the hypothesis that the duplication of an ancestral Adh single gene in the family Tephritidae occurred before the emergence of the genera Bactrocera and Ceratitis, thus suggesting that Adh duplication was based on a prespeciation rather than a postspeciation event that might have involved two independent duplication events, one in each of the two genera.     

A modular domain of NifU,a nitrogen fixation cluster protein,is highly conserved in evolution

HnifU, a gene exhibiting similarity tonifU genes of nitrogen fixation gene clusters, was identified in the course of expressed sequence tag (EST) generation from a human fetal heart cDNA library. Northern blot of human tissues and polymerase chain reaction (PCR) using human genomic DNA verified that the hnifU gene represented a human gene rather than a microbial contaminant of the cDNA library. Conceptual translation of the hnifU cDNA yielded a protein product bearing 77% and 70% amino acid identity to NifU-like hypothetical proteins fromHaemophilus influenzae andSaccharomyces cerevisiae, respectively, and 40–44% identity to the N-terminal regions of NifU proteins from several diazatrophs (i.e., nitrogen-fixing organisms). Pairwise determination of amino acid identities between the NifU-like proteins of nondiazatrophs showed that these NifU-like proteins exhibited higher sequence identity to each other (63–77%) than to the diazatrophic NifU proteins (40–48%). Further, the NifU-like proteins of non-nitrogenfixing organisms were similar only to the N-terminal region of diazatrophic NifU proteins and therefore identified a novel modular domain in these NifU proteins. These findings support the hypothesis that NifU is indeed a modular protein. The high degree of sequence similarity between NifU-like proteins from species as divergent as humans andH. influenzae suggests that these proteins perform some basic cellular function and may be among the most highly conserved proteins. Correspondence to: C.-C. Liew     

Counterselection of GATC sequences in enterobacteriophages by the components of the methyl-directed mismatch repair system

Summary Weak to severe deficit of GATC sequences in the DNA of enterobacteriophages appears to be correlated with their undermethylation during growth indam ⁺ (GATC ade-methylase) bacteria. This observation is corroborated by the sequence analysis showing no evidence for site-specific mutagenicity of 6meAde. The MutH protein of the methyl-directed mismatch repair system recognizes and cleaves the undermethylated GATC sequences in the course of mismatch repair. To enquire whether the MutH function of the methyldirected mismatch repair system participates in counterselection of GATC sequences in enterobacteriophages, we have studied the yield of bacteriophage X174 containing either 0, 1, or 2 GATC sequences, in wild type,dam, andmut (H, L, S, U) Escherichia coli. Following transfection with unmethylated DNA containing two GATC sequences, a net decrease in the yield of infective particles was observed in all bacterialmutH ⁺ dam^– strains, whereas no detectable decrease was observed in bacteria infected by DNA without GATC sequence. This effect of the MutH function is maximum in wild type andmutL andmutS bacteria whereas the effect is not significant inmutU bacteria, suggesting an interaction of the, helicase II with the MutH protein.However, indam ⁺ bacteria, the presence of GATC sequences leads to an increased yield of infective particles. The effect of GATC sequence and its Dam methylation system on phage yield inmutH ^– bacteria reveals that methylated GATC sequences are advantageous to the phage. These results suggest that the methyl-directed mismatch repair system, and in particular its MutH protein, may have participated in severe counterselection of GATC sequences from enterobacteriophages, presumably, by DNA cleavage or by interfering with DNA replication or packaging when GATC sequences are undermethylated. Coevolution of the Dam and MutH proteins could then account for the loss of GATC sequences from DNA of bacteriophages growing indam ⁺ hosts.     

Glutamine synthetase II inRhizobium: Reexamination of the proposed horizontal transfer of DNA from eukaryotes to prokaryotes

Summary We have determined the DNA sequence of aRhizobium meliloti gene that encodes glutamine synthetase II (GSII). The deduced amino acid sequence was compared to that ofBradyrhizobium japonicum GSII and those of various plant and mammalian glutamine synthetases (GS) in order to evaluate a proposal that the gene for this enzyme was recently transferred from plants to their symbiotic bacteria. There is 83.6% identity between theR. meliloti andB. japonicum proteins. The bacterial GSII proteins average 42.5% identity with the plant GS proteins and 41.8% identity with their mammalian counterparts. The plant proteins average 53.7% identity with the mammalian proteins. Thus, the GS proteins are highly conserved and the divergence of these proteins is proportional to the phylogenetic divergence of the organisms from which the sequences were determined. No transfer of genes across large taxonomic gaps is needed to explain the presence of GSII in these bacteria.     

Rampant horizontal transfer and duplication of rubisco genes in eubacteria and plastids

Previous work has shown that molecular phylogenies of plastids, cyanobacteria, and proteobacteria based on the rubisco (ribulose-1,5- bisphosphate carboxylase/oxygenase) genes rbcL and rbcS are incongruent with molecular phylogenies based on other genes and are also incompatible with structural and biochemical information. Although it has been much speculated that this is the consequence of a single horizontal gene transfer (of a proteobacterial or mitochondrial rubisco operon into plastids of rhodophytic and chromophytic algae), neither this hypothesis nor the alternative hypothesis of ancient gene duplication have been examined in detail. We have conducted phylogenetic analyses of all available bacterial rbcL sequences, and representative plastid sequences, in order to explore these alternative hypothesis and fully examine the complexity of rubisco gene evolution. The rbcL phylogeny reveals a surprising number of gene relationships that are fundamentally incongruent with organismal relationships as inferred from multiple lines of other molecular evidence. On the order of six horizontal gene transfers are implied by the form I (L8S8) rbcL phylogeny, two between cyanobacteria and proteobacteria, one between proteobacteria and plastids, and three within proteobacteria. Alternatively, a single ancient duplication of the form I rubisco operon, followed by repeated and pervasive differential loss of one operon or the other, would account for much of this incongruity. In all probability, the rubisco operon has undergone multiple events of both horizontal gene transfer and gene duplication in different lineages.      

Duplicate Genes and the Root of Angiosperms, with an Example Using Phytochrome Sequences

The root of the angiosperm tree has not yet been established. Major morphological and molecular differences between angiosperms and other seed plants have introduced ambiguities and possibly spurious results. Because it is unlikely that extant species more closely related to angiosperms will be discovered, and because relevant fossils will almost certainly not yield molecular data, the use of duplicate genes for rooting purposes may provide the best hope of a solution. Simultaneous analysis of the genes resulting from a gene duplication event along the branch subtending angiosperms would yield an unrooted network, wherein two congruent gene trees should be connected by a single branch. In these circumstances the best rooted species tree is the one that corresponds to the two gene trees when the network is rooted along the connecting branch. In general, this approach can be viewed as choosing among rooted species trees by minimizing hypothesized events such as gene duplication, gene loss, lineage sorting, and lateral transfer. Of those gene families that are potentially relevant to the angiosperm problem, phytochrome genes warrant special attention. Phylogenetic analysis of a sample of complete phytochrome (PHY) sequences implies that an initial duplication event preceded (or occurred early within) the radiation of seed plants and that each of the two resulting copies duplicated again. In one of these cases, leading to thePHYAandPHYClineages, duplication appears to have occurred before the diversification of angiosperms. Duplicate gene trees are congruent in these broad analyses, but the sample of sequences is too limited to provide much insight into the rooting question. Preliminary analyses of partialPHYAandPHYCsequences from several presumably basal angiosperm lineages are promising, but more data are needed to critically evaluate the power of these genes to resolve the angiosperm radiation.     

Evolution of structure and substrate specificity ind-alanine:d-Alanine ligases and related enzymes

Thed-alanine:d-alanine-ligase-related enzymes can have three preferential substrate specificities. Usually, these enzymes synthesized-alanyl-d-alanine. In vancomycin-resistant Gram-positive bacteria, structurally related enzymes synthesized-alanyl-d-lactate or Dalanyl-d-serine. The sequence of internal fragments of eight structurald-alanine:d-alanine ligase genes from enterococci has been determined. Alignment of the deduced amino acid sequences with those of other related enzymes from Gram-negative and Gram-positive bacteria revealed the presence of four distinct sequence patterns in the putative substrate-binding sites, each correlating with specificity to a particular substrate (d-alanine:d-lactate ligases exhibited two patterns). Phylogenetic analysis showed different clusters. The enterococcal subtree was largely superimposable on that derived from 16S rRNA sequences. In lactic acid bacteria, structural divergence due to differences in substrate specificity was observed. Glycopeptide resistance proteins VanA and VanB, the VanC-type ligases, and Dd1A and DdlB from enteric bacteria andHaemophilus influenzae constituted separate clusters. Correspondence to: P. Courvalin     

Phylogenetic analysis of the acetyl-CoA carboxylase and 3-phosphoglycerate kinase loci in wheat and other grasses

We have applied a two-gene system based on the sequences of nuclear genes encoding multi-domain plastid acetyl-CoA carboxylase (ACCase) and plastid 3-phosphoglycerate kinase (PGK) to study grass evolution. Our analysis revealed that these genes are single-copy in most of the grass species studied, allowing the establishment of orthologous relationships between them. These relationships are consistent with the known facts of their evolution: the eukaryotic origin of the plastid ACCase, created by duplication of a gene encoding the cytosolic multi-domain ACCase gene early in grass evolution, and the prokaryotic (endosymbiont) origin of the plastid PGK. The major phylogenetic relationships among grasses deduced from the nucleotide sequence comparisons of ACCase and PGK genes are consistent with each other and with the milestones of grass evolution revealed by other methods. Nucleotide substitution rates were calculated based on multiple pairwise sequence comparisons. On a relative basis, with the divergence of the Pooideae and Panicoideae subfamilies set at 60 million years ago (MYA), events leading to the Triticum/Aegilops complex occurred at the following intervals: divergence of Lolium (Lolium rigidum) at 35 MYA, divergence of Hordeum (Hordeum vulgare) at 11 MYA and divergence of Secale (Secale cereale) at 7 MYA. On the same scale, gene duplication leading to the multi-domain plastid ACCase in grasses occurred at 129 MYA, divergence of grass and dicot plastid PGK genes at 137 MYA, and divergence of grass and dicot cytosolic PGK genes at 155 MYA. The ACCase and PGK genes provide a well-understood two-locus system to study grass phylogeny, evolution and systematics.     

Gene Transfer of a Part of a β-Lactamase Gene?

β-Lactamase is an enzyme which catalyzes the hydrolysis of the β-lactam ring of penicillins and cephalosporins. By similarity analysis of amino acid sequences in a database, the amino acid sequence deduced from the nucleotide sequence of the upstream region of cytochrome c oxidase subunit II from Paracoccus denitrificans was found to have an unusually high score of homology to that of a portion of β-lactamases from Gram-negative bacteria. Furthermore, the nucleotide sequences corresponding only to this region had a very high score of similarity among them. The phylogenetic tree constructed on the basis of the amino acid sequences was in accord with that constituted on the 5S rRNA's. Moreover, the molar G + C contents and the codon usage were similar to those in their respective bacteria. It is suggested, therefore, that the nucleotide sequence in P. denitrificans was positioned by a transfer of a part of a β-lactamase gene formed as a result of gene duplication or it was formed by a deletion of the essential region of the β-lactamase gene, although no β-lactamase gene has been yet detected in P. denitrificans.     

Molecular and genetic characterization of superoxide dismutase in Haemophilus influenzae type b

Oxygen free radicals present a serious potential threat to microbial survival, through their ability to inflict Indiscriminate damage on proteins and DNA. Superoxide dismutase (SOD, EC 1.15.1.1), among other oxygen-metabolizing enzymes, is essential to prevent these toxic molecules from accumulating in the bacterial cytosol during aerobic metabolism. The gene sodA, encoding manganese-containing SOD ([Mn]-SOD), has been cloned from a virulent strain of Haemophilus influenzae type b using degenerate oligonucleotides encoding regions of the gene conserved across different bacterial species. The gene product has been identified as [Mn]-SOD by its similarity at key amino acid residues to known examples of the enzyme, by expression of enzymatically active protein from cloned DNA expressed in Escherichia coli, and by demonstration that an in-frame deletion in the gene abolishes this activity. In contrast to the situation in E. coli, this [Mn]-SOD is the only active SOD detected in H. influenzae. In further contrast to E. coli, [Mn]-SOD gene expression in H. influenzae has been found to be only partially repressed under anaerobic conditions. When expressed in E. coli the gene is regulated by Fur and Fnr, and the promoter region, identified experimentally, has been found to contain nucleotide sequence motifs similar to the Fur- and Fnr-binding sequences of E. coli, suggesting the involvement of analogues of these aerobiosis- responsive activators in H. influenzae gene expression.     

Characterization of Two Alcohol Dehydrogenase (Adh) Loci from the Olive Fruit Fly, Bactrocera (Dacus) oleae and Implications for Adh Duplication in Dipteran Insects

We report the cloning and structural characterization of two Adh loci of the olive fruit fly, Bactrocera oleae. Each of the two genes, named Adh1 and Adh2, consists of three exons and two introns for a total length of 1981 and 988 nucleotides, respectively. Their deduced amino acid sequences of 257 and 258 residues exhibit a 77% identity and display the characteristics of the insect ADH enzymes, which belong to the short-chain dehydrogenases/reductases family. The Adh genes of B. oleae are compared to the two genes of the Mediterranean fly, Ceratitis capitata, the only other species of the Tephritidae family in which the Adh genes have been studied. On the basis of amino acid divergence the four genes form two clusters each containing one gene from each species, as expected if there was one duplication event before speciation. On the basis of nucleotide sequence the four sequences form two clusters each containing the two sequences from the same species, as expected if there was a separate duplication event in each species. To help decide between the two alternatives, we compared at both the amino acid and DNA level the Adh genes of five Drosophila species that are known to carry two such genes and observed that, with only one exception at the amino acid level, conspecific loci cluster together. We conclude that the information we have at present does not allow a firm choice between the hypothesis of a single duplication event that occurred before the split of Bactrocera and Ceratitis from their common ancestor and the hypothesis of two independent duplication events, one in each of the two genera. Received: 30 May 2000 / Accepted: 17 August 2000     

Speciation in <Emphasis Type="Italic">Chlamydia</Emphasis>: Genomewide Phylogenetic Analyses Identified a Reliable Set of Acquired Genes

Horizontal gene transfer (HGT), a process through which genomes acquire sequences from distantly related organisms, is believed to be a major source of genetic diversity in bacteria. A central question concerning the impact of HGT on bacterial genome evolution is the proportion of horizontally transferred sequences within genomes. This issue, however, remains unresolved because the various methods developed to detect potential HGT events identify different sets of genes. The present-day consensus is that phylogenetic analysis of individual genes is still the most objective and accurate approach for determining the occurrence and directionality of HGT. Here we present a genome-scale phylogenetic analysis of protein-encoding genes from five closely related Chlamydia, identifying a reliable set of sequences that have arisen via HGT since the divergence of the Chlamydia lineage. According to our knowledge, this is the first systematic phylogenetic inference-based attempt to establish a reliable set of acquired genes in a bacterial genome. Although Chlamydia are obligate intracellular parasites of higher eukaryotes, and thus suspected to be isolated from HGT more than the free-living species, our results show that their diversification has involved the introduction of foreign sequences into their genome. Furthermore, we also identified a complete set of genes that have undergone deletion, duplication, or rearrangement during this evolutionary period leading to the radiation of Chlamydia species. Our analysis may provide a deeper insight into how these medically important pathogens emerged and evolved from a common ancestor.     

Osmoregulation in the model organismEscherichia coli: genes governing the synthesis of glycine betaine and trehalose and their use in metabolic engineering of stress tolerance

Glycine betaine is known to be the preferred osmoprotectant in many bacteria, and glycine betaine accumulation has also been correlated with increased cold tolerance. Trehalose is often a minor osmoprotectant in bacteria and it is a major determinant for desiccation tolerance in many so-called anhydrobiotic organisms such as baker's yeast(Saccharomyces cerevisiae). Escherichia coli has two pathways for synthesis of these protective molecules; i.e., a two-step conversion of UDP-glucose and glucose-6-phosphate to trehalose and a two-step oxidation of externally-supplied choline to glycine betaine. The genes governing the choline-to-glycine betaine pathway have been studied inE. coli and several other bacteria and higher plants. The genes governing UDP-glucose-dependent trehalose synthesis have been studied inE. coli andS. cerevisiae. Because of their well-documented function in stress protection, glycine betaine and trehalose have been identified as targets for metabolic engineering of stress tolerance. Examples of this experimental approach include the expression of theE. coli betA andArthrobacter globiformis codA genes for glycine betaine synthesis in plants and distantly related bacteria, and the expression of theE. coli otsA and yeastTPS1 genes for trehalose synthesis in plants. The published data show that glycine betaine synthesis protects transgenic plants and phototrophic bacteria against stress caused by salt and cold. Trehalose synthesis has been reported to confer increased drought tolerance in transgenic plants, but it causes negative side effects which is of concern. Thus, the much-used model organismE. coli has now become a gene resource for metabolic engineering of stress tolerance.     

APPEARANCE AND SWEEP OF A GENE DUPLICATION: ADAPTIVE RESPONSE AND POTENTIAL FOR NEW FUNCTIONS IN THE MOSQUITO CULEX PIPIENS

Evolution of a new gene function is a fundamental process of adaptation. Gene duplication followed by divergence due to relaxed selection on redundant copies has been viewed as the predominant mechanism involved in this process. At a macroevolutionary scale, evidence for this scenario came from the analysis of sequences of genes families. However, even if several genetic models have described the different potential microevolutionary scenario for a new function to evolve, little is really known about the initial evolutionary dynamics of such processes. We analyze such early dynamics in natural populations of the mosquito Culex pipiens polymorphic for a duplication at Ace.1, a locus involved in insecticide resistance. The date of occurrence and the selective advantages of the duplication were estimated using frequency data. We propose a scenario where the spread of a duplication is driven, from the very beginning, by selection due to insecticide treatment.