Glutamine synthetase gene evolution: a good molecular clock

glnE is an essential gene in Mycobacterium tuberculosis

Mycobacterium tuberculosis possesses a homologue of glnE, potentially encoding a regulator of glutamine synthetase activity. We attempted to construct glnE-disrupted mutants using a two-step strategy, whereby a single-crossover strain was first isolated, followed by sacB counterselection to isolate the double-crossover strain. Of 192 sucrose-resistant colonies tested, none were mutants, although the wild-type double crossover could be easily isolated. When a second copy of the wild-type glnE was integrated into the chromosome, we could isolate both wild-type and mutant double-crossover strains. Thus, the chromosomal gene could only be replaced with a disrupted copy when another functional copy of the gene was provided, demonstrating that this gene is essential under the conditions tested.

Two genes encoding distinct cytosolic glutamine synthetases are closely linked in the pine genome

The major isoenzyme of glutamine synthetase found in leaves of angiosperms is the chloroplastic form. However, pine seedlings contain two cytosolic glutamine synthetases in green cotyledons: GS1a, the predominant isoform, and GS1b, a minor enzyme whose relative amount is increased following phosphinotricin treatment. We have cloned a GS1b cDNA, and comparison with the previously reported GS1a cDNA sequence indicated that they correspond to separate cytosolic GS genes encoding distinct protein products. Phylogenetic analysis showed that the newly reported sequence is closer to cytosolic angiosperm GS than to GS1a, suggesting therefore that GS1a could be a divergent gymnospermous GS1 gene. Gene mapping using a F2 family of maritime pine showed co-localization of both GS genes on group 2 of the genetic linkage map. This result supports the proposed origin of different members of the GS1 family by adjacent gene duplication. The implications for gymnosperm genome organization are discussed.

Site-directed mutagenesis of Cys-92 from the alpha-polypeptide of Phaseolus vulgaris glutamine synthetase reveals that this highly conserved residue is not essential for enzyme activity but it is involved in thermal stability

The residue Cys-92 from the alpha-polypeptide of Phaseolus vulgaris glutamine synthetase is a highly conserved residue in prokaryotic and eukaryotic glutamine synthetase genes. This cysteine residue was previously proposed as a good candidate for being essential for enzyme activity. We have examined through heterologous expression in Escherichia coli and site-directed mutagenesis the functional importance of this residue. We have found that the thiol group of Cys-92 is not essential either for glutamine synthetase biosynthetic or transferase enzyme activities. The characteristic inhibition by p-hydroxymercuribenzoate (a specific sulphydryl reagent) was not substantially altered as a consequence of replacement of Cys-92 by Ala. Immunoreactivity of the glutamine synthetase mutant protein, examined both under native and denaturing conditions, was similar to the wild-type, indicating that no significant conformational changes were produced as a consequence of the introduced mutation. However, the mutant enzyme C92A was considerably less stable than the wild-type. These results indicate that Cys-92 is not an essential residue for enzyme activity but it is important for stability of the glutamine synthetase protein.

Cloning of two glutamate dehydrogenase cDNAs from Asparagus officinalis: sequence analysis and evolutionary implications

Two different amplification products, termed c1 and c2, showing a high similarity to glutamate dehydrogenase sequences from plants, were obtained from Asparagus officinalis using two degenerated primers and RT-PCR (reverse transcriptase polymerase chain reaction). The genes corresponding to these cDNA clones were designated aspGDHA and aspGDHB. Screening of a cDNA library resulted in the isolation of cDNA clones for aspGDHB only. Analysis of the deduced amino acid (aa) sequence from the full-length cDNA suggests that the gene product contains all regions associated with metabolic function of NAD glutamate dehydrogenase (NAD-GDH). A first phylogenetic analysis including only GDHs from plants suggested that the two GDH genes of A. officinalis arose by an ancient duplication event, pre-dating the divergence of monocots and dicots. Codon usage analysis showed a bias towards A/T ending codons. This tendency is likely due to the biased nucleotide composition of the asparagus genome, rather than to the translational selection for specific codons. Using principal coordinate analysis, the evolutionary relatedness of plant GDHs with homologous sequences from a large spectrum of organisms was investigated. The results showed a closer affinity of plant GDHs to GDHs of thermophilic archaebacterial and eubacterial species, when compared to those of unicellular eukaryotic fungi. Sequence analysis at specific amino acid signatures, known to affect the thermal stability of GDH, and assays of enzyme activity at non-physiological temperatures, showed a greater adaptation to heat-stress conditions for the asparagus and tobacco enzymes compared with the Saccharomyces cerevisiae enzyme.

The glutamine synthetases of rhizobia: phylogenetics and evolutionary implications

Glutamine synthetase exists in at least two related forms, GSI and GSII, the sequences of which have been used in evolutionary molecular clock studies. GSI has so far been found exclusively in bacteria, and GSII has been found predominantly in eukaryotes. To date, only a minority of bacteria, including rhizobia, have been shown to express both forms of GS. The sequences of equivalent internal fragments of the GSI and GSII genes for the type strains of 16 species of rhizobia have been determined and analyzed. The GSI and GSII data sets do not produce congruent phylogenies with either neighbor-joining or maximum-likelihood analyses. The GSI phylogeny is broadly congruent with the 16S rDNA phylogeny for the same bacteria; the GSII phylogeny is not. There are three striking rearrangements in the GSII phylograms, all of which might be explained by horizontal gene transfer to Bradyrhizobium (probably from Mesorhizobium), to Rhizobium galegae (from Rhizobium), and to Mesorhizobium huakuii (perhaps from Rhizobium). There is also evidence suggesting intrageneric DNA transfer within Mesorhizobium. Meta-analysis of both GS genes from the different genera of rhizobia and other reference organisms suggests that the divergence times of the different rhizobium genera predate the existence of legumes, their host plants.

Functional importance of Asp56 from the alpha-polypeptide of Phaseolus vulgaris glutamine synthetase. An essential residue for transferase but not for biosynthetic enzyme activity

Replacement of Asp56 by site-directed mutagenesis of the alpha-gene from Phaseolus vulgaris glutamine synthetase heterologously expressed in Escherichia coli produces a complete loss of transferase enzyme activity, thus revealing essentiality of the residue for this particular enzyme activity. This happens independent of Asp56 being replaced by Ala or Glu, suggesting that the essentiality of this residue cannot be attributed to its negative electrical charge. However, a high level of glutamine synthetase biosynthetic specific activity (referred to glutamine synthetase protein, as determined immunologically), is present in D56A and D56E mutants, suggesting that Asp56 is an example of a residue that has a different role in the catalytic mechanism of both enzyme activities of this protein. Km for ATP, glutamate and Mg2+, as well as energy of activation, can be altered as a consequence of the performed mutations. However, the Km and catalytic efficiency for ammonium remains unaffected. Therefore, the catalytic role of Asp56 in the alpha-polypeptide of higher plant glutamine synthetase is quite different from the role proposed for its highly conserved homologue in bacteria (Asp50 in E. coli), which has been associated with binding and deprotonation of ammonium. On the other hand, we also show other results indicating that Asp56 is important in the spatial conformation of the active site and/or the protein, Asp56 being a crucial residue in the salting-out aggregation properties of the enzyme.

Chloroplast-expressed glutamine synthetase (ncpGS): potential utility for phylogenetic studies with an example from Oxalis (Oxalidaceae)

Chloroplast-expressed glutamine synthetase (ncpGS), a nuclear-encoded gene containing several introns, is introduced as a tool for phylogenetic studies at lower taxonomic levels. This gene is a member of a multigene family, but it diverged long ago from the cytosolic-expressed members of the family and appears to be single copy in the majority of taxa examined to date. The conservation of both coding sequence and position of introns has allowed the design of primers for use in a broad range of dicot taxa to amplify and sequence a region of ncpGS that contains four introns. The utility of this region in phylogenetic studies of congeneric species is illustrated by an example using eight Oxalis species. The four introns in these taxa are typical in size (76 to 136 bp), base composition (high T content), and structure (e.g., sequence of splice sites and putative branch points) for plant internal introns. Levels of variation among these ncpGS sequences compare favorably with those of the internal transcribed spacer of nuclear ribosomal DNA (ITS) from the same taxa, and results of phylogenetic analysis of ncpGS data are generally congruent with previous results using ITS.

Site-directed mutagenesis of Glu-297 from the alpha-polypeptide of Phaseolus vulgaris glutamine synthetase alters kinetic and structural properties and confers resistance to L-methionine sulfoximine

In this paper we examine the functionality of Glu-297 from the alpha-polypeptide of Phaseolus vulgaris glutamine synthetase (EC 6.3.1.2). For this purpose, the gln alpha cDNA was recombinantly expressed in Escherichia coli, and site-directed mutants constructed, in which this residue was replaced by alanine. The level of glutamine synthetase transferase catalytic activity in the mutant strain was 70-fold lower while biosynthetic activity remained practically unaffected. Kinetic parameters for both enzyme activities were not greatly altered except for the Km for ammonium in biosynthetic activity, which increased 100-fold. A similar result was reported when mutagenizing Glu-327 from E. coli glutamine synthetase, a residue shown to be present at the active site. This suggests that the Glu residue mutated in the higher-plant enzyme could develop a similar catalytic role to that of bacteria. Another characteristic feature of the mutant protein was its higher resistance to inhibition of the biosynthetic activity by L-methionine sulfoximine, a typical inhibitor of glutamine synthetase. In addition, we show that immunoreactivity of the glutamine synthetase mutant protein, both under native and denaturing conditions, is similar to the wild type, indicating that no deep conformational changes were produced as a consequence of the introduced mutation. However, structural changes in the active site can be predicted from alterations detected in the behaviour of the mutant protein towards affinity chromatography on 2',5'-ADP-Sepharose, as compared to the wild type. Nevertheless, complementation of an E. coli glnA mutation indicated that the E297A mutant enzyme was physiologically functional.

Oxidative turnover of soybean root glutamine synthetase. In vitro and in vivo studies

Glutamine synthetase (GS) is the key enzyme in ammonia assimilation and catalyzes the ATP-dependent condensation of NH3 with glutamate to produce glutamine. GS in plants is an octameric enzyme. Recent work from our laboratory suggests that GS activity in plants may be regulated at the level of protein turnover (S.J. Temple, T.J. Knight, P.J. Unkefer, C. Sengupta-Gopalan [1993] Mol Gen Genet 236: 315-325; S.J. Temple, S. Kunjibettu, D. Roche, C. Sengupta-Gopalan [1996] Plant Physiol 112: 1723-1733; S.J. Temple, C. Sengupta-Gopalan [1997] In C.H. Foyer, W.P. Quick, eds, A Molecular Approach to Primary Metabolism in Higher Plants. Taylor & Francis, London, pp 155-177). Oxidative modification of GS has been implicated as the first step in the turnover of GS in bacteria. By incubating soybean (Glycine max) root extract enriched in GS in a metal-catalyzed oxidation system to produce the.OH radical, we have shown that GS is oxidized and that oxidized GS is inactive and more susceptible to degradation than nonoxidized GS. Histidine and cysteine protect GS from metal-catalyzed inactivation, indicating that oxidation modifies the GS active site and that cysteine and histidine residues are the site of modification. Similarly, ATP and particularly ATP/glutamate give the enzyme the greatest protection against oxidative inactivation. The roots of plants fed ammonium nitrate showed a 3-fold increase in the level of GS polypeptides and activity compared with plants not fed ammonium nitrate but without a corresponding increase in the GS transcript level. This would suggest either translational or posttranslational control of GS levels.

Nitrogen control of the glnN gene that codes for GS type III, the only glutamine synthetase in the cyanobacterium Pseudanabaena sp. PCC 6903

Pseudanabaena sp. strain PCC 6903 is the first cyanobacteria lacking the typical prokaryotic glutamine synthetase type I encoded by the glnA gene. The glnN gene product, glutamine synthetase type III, is the only glutamine synthetase activity present in this cyanobacterium. Analysis of glnN expression clearly indicated a nitrogen-dependent regulation. Pseudanabaena glnN gene expression and GSIII activity were upregulated under nitrogen starvation or using nitrate as a nitrogen source, while low levels of transcript and activity were found in ammonium-containing medium. Primer extension analysis showed that the glnN gene promoter structure resembled that of the NtcA-related promoters. Mobility shift assays demonstrated that Synechocystis sp. PCC 6803 NtcA protein, expressed and purified from Escherichia coli, bound to the promoter of the Pseudanabaena 6903 glnN gene. The NtcA control of the glnN gene in this cyanobacterium suggested that, in the absence of a glnA gene, NtcA took control of the only glutamine synthetase gene in a fashion similar to the way the glnA gene is governed in those cyanobacteria harbouring a glnA gene.

Regulation of the spatiotemporal pattern of expression of the glutamine synthetase gene

Glutamine synthetase, the enzyme that catalyzes the ATP-dependent conversion of glutamate and ammonia into glutamine, is expressed in a tissue-specific and developmentally controlled manner. The first part of this review focuses on its spatiotemporal pattern of expression, the factors that regulate its levels under (patho)physiological conditions, and its role in glutamine, glutamate, and ammonia metabolism in mammals. Glutamine synthetase protein stability is more than 10-fold reduced by its product glutamine and by covalent modifications. During late fetal development, translational efficiency increases more than 10-fold. Glutamine synthetase mRNA stability is negatively affected by cAMP, whereas glucocorticoids, growth hormone, insulin (all positive), and cAMP (negative) regulate its rate of transcription. The signal transduction pathways by which these factors may regulate the expression of glutamine synthetase are briefly discussed. The second part of the review focuses on the evolution, structure, and transcriptional regulation of the glutamine synthetase gene in rat and chicken. Two enhancers (at -6.5 and -2.5 kb) were identified in the upstream region and two enhancers (between +156 and +857 bp) in the first intron of the rat glutamine synthetase gene. In addition, sequence analysis suggests a regulatory role for regions in the 3' untranslated region of the gene. The immediate-upstream region of the chicken glutamine synthetase gene is responsible for its cell-specific expression, whereas the glucocorticoid-induced developmental appearance in the neural retina is governed by its far-upstream region.

Archaea and the prokaryote-to-eukaryote transition

Since the late 1970s, determining the phylogenetic relationships among the contemporary domains of life, the Archaea (archaebacteria), Bacteria (eubacteria), and Eucarya (eukaryotes), has been central to the study of early cellular evolution. The two salient issues surrounding the universal tree of life are whether all three domains are monophyletic (i.e., all equivalent in taxanomic rank) and where the root of the universal tree lies. Evaluation of the status of the Archaea has become key to answering these questions. This review considers our cumulative knowledge about the Archaea in relationship to the Bacteria and Eucarya. Particular attention is paid to the recent use of molecular phylogenetic approaches to reconstructing the tree of life. In this regard, the phylogenetic analyses of more than 60 proteins are reviewed and presented in the context of their participation in major biochemical pathways. Although many gene trees are incongruent, the majority do suggest a sisterhood between Archaea and Eucarya. Altering this general pattern of gene evolution are two kinds of potential interdomain gene transferrals. One horizontal gene exchange might have involved the gram-positive Bacteria and the Archaea, while the other might have occurred between proteobacteria and eukaryotes and might have been mediated by endosymbiosis.

Cloning of the glutamine synthetase gene from group B streptococci

The glnA gene from the human pathogen Streptococcus agalactiae was cloned from a genomic library prepared with the lambda phage vector lambdaDASHII. A 4.6-kb DNA fragment of one of the recombinant phages was subcloned in pUC18. This Escherichia coli clone expressed a 52-kDa protein encoded by a 1,341-bp open reading frame. The nucleotide sequence of the open reading frame and the deduced amino acid sequence shared a significant degree of homology with the sequences of other glutamine synthetases (GS). The highest homology was between our deduced protein and GS of gram-positive bacteria such as Bacillus subtilis, Bacillus cereus, and Staphylococcus aureus. Plasmids with the cloned streptococcal glnA were able to complement E. coli glnA mutants grown on minimal media. Rabbit antisera to streptococcal GS recombinant protein recognized not only the recombinant protein but also a similar-sized band in mutanolysin extracts of all group B streptococcal strains tested, regardless of polysaccharide type or surface protein profile. The amino acid sequence of the deduced protein had similarities to other streptococcal cell-surface-bound proteins. The possible functional role of the immunological features of streptococcal GS is discussed.

Molecular evolution of nitrate reductase genes

To understand the evolutionary mechanisms and relationships of nitrate reductases (NRs), the nucleotide sequences encoding 19 nitrate reductase (NR) genes from 16 species of fungi, algae, and higher plants were analyzed. The NR genes examined show substantial sequence similarity, particularly within functional domains, and large variations in GC content at the third codon position and intron number. The intron positions were different between the fungi and plants, but conserved within these groups. The overall and nonsynonymous substitution rates among fungi, algae, and higher plants were estimated to be 4.33 x 10(-10) and 3.29 x 10(-10) substitutions per site per year. The three functional domains of NR genes evolved at about one-third of the rate of the N-terminal and the two hinge regions connecting the functional domains. Relative rate tests suggested that the nonsynonymous substitution rates were constant among different lineages, while the overall nucleotide substitution rates varied between some lineages. The phylogenetic trees based on NR genes correspond well with the phylogeny of the organisms determined from systematics and other molecular studies. Based on the nonsynonymous substitution rate, the divergence time of monocots and dicots was estimated to be about 340 Myr when the fungi-plant or algae-higher plant divergence times were used as reference points and 191 Myr when the rice-barley divergence time was used as a reference point. These two estimates are consistent with other estimates of divergence times based on these reference points. The lack of consistency between these two values appears to be due to the uncertainty of the reference times.

Molecular classification of living organisms

Recent studies in molecular evolution have generated strong conflicts in opinion as to how world living organisms should be classified. The traditional classification of life into five kingdom has been challenged by the molecular analysis carried out mostly on rRNA sequences, which supported the division of the extant living organisms into three major groups: Archaebacteria, Eubacteria, and Eukaryota. As to the problem of placing the root of the tree of life, the analysis carried out on a few genes has provided discrepant results. In order to measure the genetic distances between species, we have carried out an evolutionary analysis of the glutamine synthetase genes, which previously have been revealed to be good molecular clocks, and of the small and large rRNA genes. All data demonstrate that archaebacteria are more closely related to eubacteria than to eukaryota, thus supporting the classical division of living organisms into two main superkingdoms, Prokaryota and Eukaryota.

Probing the catalytic roles of n2-site glutamate residues in Escherichia coli glutamine synthetase by mutagenesis

The contribution of metal ion ligand type and charge to catalysis and regulation at the lower affinity metal ion site (n2 site) of Escherichia coli glutamine synthetase (GS) was tested by mutagenesis and kinetic analysis. The 2 glutamate residues at the n2 site, E129 and E357, were changed to E129D, E129H, E357H, E357Q, and E357D, representing conservative and nonconservative alterations. Unadenylylated and fully adenylylated enzyme forms were studied. The Mn(2+)-KD values, UV-cis and fluorescence emission properties were similar for all mutants versus WTGS, except E129H. For kinetic determinations with both Mn2+ and Mg2+, nonconservative mutants (E357H, E129H, E357Q) showed lower biosynthetic activities than conservative mutants (E129D, E357D). Relative to WTGS, all the unadenylylated Mn(2+)-activated enzymes showed reduced kcat/Km values for ATP (> 7-fold) and for glutamate (> 10-fold). Of the unadenylylated Mg(2+)-activated enzymes, only E129D showed kinetic parameters competitive with WTGS, and adenylylated E129D was a 20-fold better catalyst than WTGS. We propose the n2-site metal ion activates ADP for departure in the phosphorylation of glutamate by ATP to generate gamma-glutamyl phosphate. Alteration of the charge density at this metal ion alters the transition-state energy for phosphoryl group transfer and may affect ATP binding and/or ADP release. Thus, the steady-state kinetic data suggest that modifying the charge density increases the transition-state energies for chemical steps. Importantly, the data demonstrate that each ligand position has a specialized spatial environment and the charge of the ligand modulates the catalytic steps occurring at the metal ion. The data are discussed in the context of the known X-ray structures of GS.

Evolution of glutamine synthetase genes is in accordance with the neutral theory of molecular evolution

Evolution of glutamine synthetase gene is discussed on the results of DNA sequence analysis of the gene. Thirty DNA sequences of various organisms spanning from prokaryotes to eukaryotes were collected from the DNA data banks and translated first, they were aligned next, then evolutionary distances were computed, and molecular phylogeny was finally estimated. The results of the alignment reveal that functionally important regions of glutamine synthetase have been evolutionarily more conserved than the remaining regions. The evolutionary distances computed also show that the rate of synonymous substitution is higher than that of nonsynonymous substitution. These are well in accordance with the neutral theory of molecular evolution. Besides, the molecular phylogeny obtained shows that the origin of glutamine synthetase gene is much earlier than the divergence between eukaryotes and prokaryotes, suggesting that the gene is one of the oldest genes functioning now.

Genetic basis for tissue isozymes of glutamine synthetase in elasmobranchs

Tissue-specific isozymes of glutamine synthetase are present in elasmobranchs. A larger isozyme occurs in tissues in which the enzyme is localized in mitochondria (liver, kidney) whereas a smaller form occurs in tissues in which it is cytosolic (brain, spleen, etc.). The nucleotide sequence of spiny dogfish shark (Squalus acanthias) liver glutamine synthetase mRNA, derived from its cDNA, shows there are two in-frame initiation codons (AUG) at the N-terminus which will account for the size differences between the two isozymes. Initiation at the up-stream and down-stream sites would yield peptides of 45,406 and 41,869 mol. wts. representing the precursor of the mitochondrial isozyme and the cytosolic isozyme, respectively. The additional N-terminal 29 amino acids present in the mitochondrial isozyme precursor contains two putative cleavage sites based on the Arg-X-(Phe,Ile,Leu) motif. The predicted two-step processing would remove 14 of the 29 N-terminal amino acids. These 14 amino acids can be predicted to form a very strong amphipathic mitochondrial targeting signal. Their removal would yield a mature peptide of 43,680 mol. wt. The calculated mol. wts. based on the derived amino acid sequence are therefore in good agreement with previous estimates of an approximately 1.5-2-kDa difference between the M(r)s of the mitochondrial and cytosolic isozymes. A model for the evolution of the mitochondrial targeting of glutamine synthetase in vertebrates is proposed.

Evolutionary relationships of bacterial and archaeal glutamine synthetase genes

Glutamine synthetase (GS), an essential enzyme in ammonia assimilation and glutamine biosynthesis, has three distinctive types: GSI, GSII and GSIII. Genes for GSI have been found only in bacteria (eubacteria) and archaea (archaebacteria), while GSII genes only occur in eukaryotes and a few soil-dwelling bacteria. GSIII genes have been found in only a few bacterial species. Recently, it has been suggested that several lateral gene transfers of archaeal GSI genes to bacteria may have occurred. In order to study the evolution of GS, we cloned and sequenced GSI genes from two divergent archaeal species: the extreme thermophile Pyrococcus furiosus and the extreme halophile Haloferax volcanii. Our phylogenetic analysis, which included most available GS sequences, revealed two significant prokaryotic GSI subdivisions: GSI-alpha and GSI-beta. GSI-alpha-genes are found in the thermophilic bacterium, Thermotoga maritima, the low G+C Gram-positive bacteria, and the Euryarchaeota (includes methanogens, halophiles, and some thermophiles). GSI-beta-type genes occur in all other bacteria. GSI-alpha- and GSI-beta-type genes also differ with respect to a specific 25-amino-acid insertion and adenylylation control of GS enzyme activity, both absent in the former but present in the latter. Cyanobacterial genes lack adenylylation regulation of GS and may have secondarily lost it. The GSI gene of Sulfolobus solfataricus, a member of the Crenarchaeota (extreme thermophiles), is exceptional and could not be definitely placed in either subdivision.

Genetic, molecular and developmental analysis of the glutamine synthetase isozymes of Drosophila melanogaster

The glutamine synthetase isozymes of Drosophila melanogaster offer an attractive model for the study of the molecular genetics and evolution of a small gene family encoding enzymatic isoforms that evolved to assume a variety of specific and sometimes essential biological functions. In Drosophila melanogaster two GS isozymes have been described which exhibit different cellular localisation and are coded by a two-member gene family. The mitochondrial GS structural gene resides at the 21B region of the second chromosome, the structural gene for the cytosolic isoform at the 10B region of the X chromosome. cDNA clones corresponding to the two genes have been isolated and sequenced. Evolutionary analysis data are in accord with the hypothesis that the two Drosophila glutamine synthetase genes are derived from a duplication event that occurred near the time of divergence between Insecta and Vertebrata. Both isoforms catalyse all reactions catalysed by other glutamine synthetases, but the different kinetic parameters and the different cellular compartmentalisation suggest strong functional specialisation. In fact, mutations of the mitochondrial GS gene produce embryo-lethal female sterility, defining a function of the gene product essential for the early stages of embryonic development. Preliminary results show strikingly distinct spatial and temporal patterns of expression of the two isoforms at later stages of development.

A new type of glutamine synthetase in cyanobacteria: the protein encoded by the glnN gene supports nitrogen assimilation in Synechocystis sp. strain PCC 6803

A new glutamine synthetase gene, glnN, which encodes a polypeptide of 724 amino acid residues (M(r), 79,416), has been identified in the unicellular cyanobacterium Synechocystis sp. strain PCC 6803; this is the second gene that encodes a glutamine synthetase (GS) in this cyanobacterium. The functionality of this gene was evidenced by its ability to complement an Escherichia coli glnA mutant and to support Synechocystis growth in a strain whose glnA gene was inactivated by insertional mutagenesis. In this mutant (strain SJCR3), as well as in the wild-type strain, the second GS activity was subject to regulation by the nitrogen source, being strongly enhanced in nitrogen-free medium. Transcriptional fusion of a chloramphenicol acetyltransferase (cat) gene with the 5'-upstream region of glnN suggested that synthesis of the second Synechocystis GS is regulated at the transcriptional level. Furthermore, the level of glnN mRNA, a transcript of about 2,300 bases, was found to be strongly increased in nitrogen-free medium. The glnN product is similar to the GS subunits of Bacteroides fragilis and Butyrivibrio fibrisolvens, two obligate anaerobic bacteria whose GSs are markedly different from other prokaryotic and eukaryotic GSs. However, significant similarity is evident in the five regions which are homologous in all of the GSs so far described. The new GS gene was also found in other cyanobacteria but not in N2-fixing filamentous species.

Evolution of HSP70 gene and its implications regarding relationships between archaebacteria, eubacteria, and eukaryotes

The 70-kDa heat-shock protein (HSP70) constitutes the most conserved protein present in all organisms that is known to date. Based on global alignment of HSP70 sequences from organisms representing all three domains, numerous sequence signatures that are specific for prokaryotic and eukaryotic homologs have been identified. HSP70s from the two archaebacterial species examined (viz., Halobacterium marismortui and Methanosarcina mazei) have been found to contain all eubacterial but no eukaryotic signature sequences. Based on several novel features of the HSP70 family of proteins (viz., presence of tandem repeats of a 9-amino-acid [a.a.] polypeptide sequence and structural similarity between the first and second quadrants of HSP70, homology of the N-terminal half of HSP70 to the bacterial MreB protein, presence of a conserved insert of 23-27 a.a. in all HSP70s except those from archaebacteria and gram-positive eubacteria) a model for the evolution of HSP70 gene from an early stage is proposed. The HSP70 homologs from archaebacteria and gram-positive bacteria lacking the insert in the N-terminal quadrants are indicated to be the ancestral form of the protein. Detailed phylogenetic analyses of HSP70 sequence data (viz., by bootstrap analyses, maximum parsimony, and maximum likelihood methods) provide evidence that archaebacteria are not monophyletic and show a close evolutionary linkage with the gram-positive eubacteria. These results do not support the traditional archaebacterial tree, where a close relationship between archaebacterial and eukaryotic homologs is observed. To explain the phylogenies based on HSP70 and other gene sequences, a model for the origin of eukaryotic cells involving fusion between archaebacteria and gram-negative eubacteria is proposed.

Time and biosequences

In this paper we discuss and demonstrate the importance of several factors relative to the relationship between time and evolution of biosequences. In both quantitative and qualitative measurements of the genetic distances, the compositional constraints of the nucleotide sequences play a very important role. We demonstrate that when homologous sequences significantly differ in base composition we get erratic branching order and/or wrong evaluation of the evolutionary rates. We must consider that every gene may have a different evolutionary dynamic along its sequence, generally linked to its functional constraints; this too can seriously affect its clock-like behavior. We report some cases showing how these factors can affect the quantitative measurements of the genetic distances of biosequences.

Early evolution of photosynthesis: clues from nitrogenase and chlorophyll iron proteins

Chlorophyll (Chl) is often viewed as having preceded bacteriochlorophyll (BChl) as the primary photoreceptor pigment in early photosynthetic systems because synthesis of Chl requires one fewer enzymatic reduction than does synthesis of BChl. We have conducted statistical DNA sequence analyses of the two reductases involved in Chl and BChl synthesis, protochlorophyllide reductase and chlorin reductase. Both are three-subunit enzymes in which each subunit from one reductase shares significant amino acid identity with a subunit of the other, indicating that the two enzymes are derived from a common three-subunit ancestral reductase. The "chlorophyll iron protein" subunits, encoded by the bchL and bchX genes in the purple bacterium Rhodobacter capsulatus, also share amino acid sequence identity with the nitrogenase iron protein, encoded by nifH. When nitrogenase iron proteins are used as outgroups, the chlorophyll iron protein tree is rooted on the chlorine reductase lineage. This rooting suggests that the last common ancestor of all extant photosynthetic eubacteria contained BChl, not Chl, in its reaction center, and implies that Chl-containing reaction centers were a late invention unique to the cyanobacteria/chloroplast lineage.

Close linkage of genes encoding glutamine synthetases I and II in Frankia alni CpI1

Frankia alni CpI1 has two glutamine synthetases (GSs), GSI and GSII. The GSI gene (glnA) was isolated from a cosmid library of F. alni CpI1 DNA by heterologous probing with glnA from Streptomyces coelicolor. The glnA gene was shown to be located upstream of the GSII gene (glnII) by DNA-DNA hybridization. The nucleotide sequences of the 1,422-bp CpI1 glnA gene and of the 449-bp intervening region between glnA and glnII were determined, and the glnA amino acid sequence was deduced. In common with GSIs from other organisms, CpI1 GSI contains five conserved regions near the active site and a conserved tyrosine at the adenylylation site. F. alni CpI1 glnA complemented the glutamine growth requirement of the Escherichia coli glnA deletion strain YMC11 but only when expressed from an E. coli lac promoter. While the functional significance of maintaining two GSs adjacent to one another remains unclear, this arrangement in F. alni provides support for the recently proposed origin of GSI and GSII as resulting from a gene duplication early in the evolution of life.

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

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

Time-resolved fluorescence and computational studies of adenylylated glutamine synthetase: analysis of intersubunit interactions

Adenylylation of Tyr-397 of each subunit of Escherichia coli glutamine synthetase (GS) down-regulates enzymatic activity in vivo. The overall structure of the enzyme consists of 12 subunits arranged as two hexamers, face to face. Research reported in this paper addresses the question of whether the covalently attached adenylyl group interacts with neighboring amino acid residues to produce the regulatory phenomenon. Wild-type GS has two Trp residues (positions 57 and 158) and the adenylylation site lies within 7-8 A of the Trp-57 loop in the adjacent subunit of the same hexameric ring; Trp-158 is about 35 A from the site of adenylylation. Fluorescence lifetimes and quantum yields have been determined for two fluorophores with wild-type and mutant GS. One fluorophore is epsilon-AMP adenylylated GS (at Tyr-397), and the other fluorophore is the intrinsic protein residue Trp-57. These experiments were conducted in order to detect possible intersubunit interactions between adenylyl groups and the neighboring Trp-57 to search for a role for the Trp-57 loop in the regulation of GS. The fluorescence due to epsilon-AMP of two adenylylated enzymes, wild-type GS and the W158F mutant, exhibits heterogeneous decay kinetics; the data adequately fit to a double exponential decay model with recovered average lifetime values of 18.2 and 2.1 ns, respectively. The pre-exponential factors range from 0.66 to 0.73 for the long lifetime component, at five emission wavelengths. The W57L-epsilon-AMP enzyme yields longer average lifetime values of 19.5 and 2.4 ns, and the pre-exponential factors range from 0.82 to 0.85 for the long lifetime component. An additional residue in the Trp-57 loop, Lys-58, has been altered and the K58C mutant enzyme has been adenylylated with epsilon-AMP on Tyr-397. Lys-58 is near the ATP binding site and may represent a link by which the adenylyl group controls the activity of GS. The fluorescence of epsilon-AMP-adenylylated K58C mutant GS is best described by a triple exponential decay with average recovered lifetime values of 19.9, 4.6, and 0.58 ns, with the largest fraction being the median lifetime component. Relative quantum yields of epsilon-AMP-Tyr-397 were measured in order to determine if static quenching occurs from adenine-indole stacking in the wild-type GS. The relative quantum yield of the epsilon-AMP-adenylylated W57L mutant is larger than the wild-type protein by the amount predicted from the difference in lifetime values: thus, no static quenching is evident.(ABSTRACT TRUNCATED AT 400 WORDS)

Evolution of the glutamine synthetase gene, one of the oldest existing and functioning genes

We performed molecular phylogenetic analyses of glutamine synthetase (GS) genes in order to investigate their evolutionary history. The analyses were done on 30 DNA sequences of the GS gene which included both prokaryotes and eukaryotes. Two types of GS genes are known at present: the GSI gene found so far only in prokaryotes and the GSII gene found in both prokaryotes and eukaryotes. Our study has shown that the two types of GS gene were produced by a gene duplication which preceded, perhaps by > 1000 million years, the divergence of eukaryotes and prokaryotes. The results are consistent with the facts that (i) GS is a key enzyme of nitrogen metabolism found in all extant life forms and (ii) the oldest biological fossils date back 3800 million years. Thus, we suggest that GS genes are one of the oldest existing and functioning genes in the history of gene evolution and that GSI genes should also exist in eukaryotes. Furthermore, our study may stimulate investigation on the evolution of "preprokaryotes," by which we mean the organisms that existed during the era between the origin of life and the divergence of prokaryotes and eukaryotes.

Statistical tests of models of DNA substitution

Penny et al. have written that "The most fundamental criterion for a scientific method is that the data must, in principle, be able to reject the model. Hardly any [phylogenetic] tree-reconstruction methods meet this simple requirement." The ability to reject models is of such great importance because the results of all phylogenetic analyses depend on their underlying models--to have confidence in the inferences, it is necessary to have confidence in the models. In this paper, a test statistic suggested by Cox is employed to test the adequacy of some statistical models of DNA sequence evolution used in the phylogenetic inference method introduced by Felsenstein. Monte Carlo simulations are used to assess significance levels. The resulting statistical tests provide an objective and very general assessment of all the components of a DNA substitution model; more specific versions of the test are devised to test individual components of a model. In all cases, the new analyses have the additional advantage that values of phylogenetic parameters do not have to be assumed in order to perform the tests.

Altered regulation of the glnRA operon in a Bacillus subtilis mutant that produces methionine sulfoximine-tolerant glutamine synthetase

A Bacillus subtilis mutant that produced glutamine synthetase (GS) with altered sensitivity to DL-methionine sulfoximine was isolated. The mutation, designated glnA33, was due to a T.A-to-C.G transition, changing valine to alanine at codon 190 within the active-site C domain. Altered regulation was observed for GS activity and antigen and mRNA levels in a B. subtilis glnA33 strain. The mutant enzyme was 28-fold less sensitive to DL-methionine sulfoximine and had a 13.0-fold-higher Km for hydroxylamine and a 4.8-fold-higher Km for glutamate than wild-type GS did.

Modulation of glutamine synthetase gene expression in tobacco by the introduction of an alfalfa glutamine synthetase gene in sense and antisense orientation: molecular and biochemical analysis

A glutamine synthetase (GS) cDNA isolated from an alfalfa cell culture cDNA library was found to represent a cytoplasmic GS. The full-length alfalfa GS1 coding sequence, in both sense and antisense orientation and under the transcriptional control of the cauliflower mosaic virus 35S promoter, was introduced into tobacco. Leaves of tobacco plants transformed with the sense construct contained greatly elevated levels of GS transcript and GS polypeptide which assembled into active enzyme. Leaves of the plants transformed with the antisense GS1 construct showed a significant decrease in the level of both GS1 and GS2 polypeptides and GS activity, but did not show any significant decrease in the level of endogenous GS mRNA. We have proposed that antisense inhibition using a heterologous antisense GS RNA occurs at the level of translation. Our results also suggest that the post-translational assembly of GS subunits into a holoenzyme requires an additional factor(s) and is under regulatory control.

The molecular clock ticks regularly in muroid rodents and hamsters

Extensive DNA sequence data are used to compare the rates of nucleotide substitution in the mouse, rat, and hamster lineages. A relative rate test using hamster sequences as references shows that the rates of synonymous and nonsynonymous substitution in the mouse and rat lineages are nearly equal and a test using human sequences as references shows that the rates in the mouse, rat, and hamster lineages are also nearly equal. Under the assumptions that the guinea pig lineage and the myomorph (mouse, rat, and hamster) lineage diverged 70-100 million years (Myr) ago and that the rate of nucleotide substitution has been constant in all these lineages since their divergence, the date of the mouse-rat split is estimated to be between 20 and 29 Myr ago, which is considerably older than the date (approximately 12 Myr) suggested by available rodent fossils and considerably younger than the date (approximately 35 Myr) suggested by Wilson and colleagues. The murid-hamster split is estimated to be 1.6 times older than the mouse-rat split.

Diversification of the Wnt gene family on the ancestral lineage of vertebrates

Diversification of the Wnt genes, a family of powerful developmental regulator molecules, is inferred by molecular evolutionary analyses. Fifty-five recently determined partial sequences from a variety of vertebrates and invertebrates, together with 17 published sequences, mostly from the mouse and Drosophila melanogaster, are analyzed. Wnt-1 through -7 originated before the last common ancestor of arthropods and deuterostomes lived. Another round of gene duplication, involving Wnt-3, -5, -7, and -10, occurred after the echinoderm lineage arose, on the ancestral lineage of jawed vertebrates. Increased constraints were imposed on the Wnt genes when jawed vertebrates originated, as indicated by an overall 4-fold lower rate of amino acid replacements in jawed vertebrates compared with invertebrates and jawless vertebrates. The Wnt genes are thus inferred to have undergone a disproportionately high amount of structural and functional evolution in the relatively short time (approximately 100 million years) between the origin of the echinoderm lineage and the first diversification of jawed vertebrates. A model is presented for the relationship of functional diversification of developmental regulators and their rates of amino acid replacement.

Sequence of the GLN1 gene of Saccharomyces cerevisiae: role of the upstream region in regulation of glutamine synthetase expression

The GLN1 gene, encoding glutamine synthetase in Saccharomyces cerevisiae, was sequenced, and its encoded polypeptide was shown to have significant homology to other eukaryotic glutamine synthetases. S1 analysis has defined the transcriptional start site of the gene. Upstream analysis of the gene using lacZ fusions has verified transcriptional control of the gene and has identified a nitrogen upstream activation sequence which is required for the increased transcription of GLN1 seen when glutamine is replaced by glutamate as the nitrogen source. cis-acting sites required for the increased transcription in response to purine starvation also have been localized.

A simple method for global sequence comparison

A simple method of sequence comparison, based on a correlation analysis of oligonucleotide frequency distributions, is here shown to be a reliable test of overall sequence similarity. The method does not involve sequence alignment procedures and permits the rapid screening of large amounts of sequence data. It identifies those sequences which deserve more careful analysis of sequence similarity at the level of resolution of the single nucleotide. It uses observed quantities only and does not involve the adoption of any theoretical model.

Nucleotide polymorphism and evolution in the glyceraldehyde-3-phosphate dehydrogenase gene (gapA) in natural populations of Salmonella and Escherichia coli

Nucleotide sequences of the gapA gene, encoding the glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase, were determined for 16 strains of Salmonella and 13 strains of Escherichia coli recovered from natural populations. Pairs of sequences from strains representing the eight serovar groups of Salmonella differed, on average, at 3.8% of nucleotide sites and 1.1% of inferred amino acids, and comparable values for E. coli were an order of magnitude smaller (0.2% and 0.1%, respectively). The rate of substitution at synonymous sites was significantly higher for codons specifying the catalytic domain of the enzyme than for those encoding the NAD(+)-binding domain, but the nonsynonymous substitution rate showed the opposite relationship. For Salmonella, statistical tests for nonrandom clustering of polymorphic sites failed to provide evidence that intragenic recombination or gene conversion has contributed to the generation of allelic diversity. The topology of a tree constructed from the gapA sequences was generally similar to that of phylogenetic trees of the strains based on multilocus enzyme electrophoresis, but the level of divergence of gapA in Salmonella group V from other Salmonella and E. coli strains is much greater than that indicated by DNA hybridization for the genome as a whole.