Participation of UV-regulated Genes in the Response to Helix-distorting DNA Damage in the Thermoacidophilic Crenarchaeon Sulfolobus acidocaldarius.

Several species of Sulfolobales have been used as model organisms in the study of response mechanisms to ultraviolet (UV) irradiation in hyperthermophilic crenarchaea. To date, the transcriptional responses of genes involved in the initiation of DNA replication, transcriptional regulation, protein phosphorylation, and hypothetical function have been observed in Sulfolobales species after UV irradiation. However, due to the absence of knockout experiments, the functions of these genes under in situ UV irradiation have not yet been demonstrated. In the present study, we constructed five gene knockout strains (cdc6-2, tfb3, rio1, and two genes encoding the hypothetical proteins, Saci_0951 and Saci_1302) of Sulfolobus acidocaldarius and examined their sensitivities to UV irradiation. The knockout strains exhibited significant sensitivities to UV-B irradiation, indicating that the five UV-regulated genes play an important role in responses to UV irradiation in vivo. Furthermore, Δcdc6-2, Δrio1, ΔSaci_0951, and Δtfb3 were sensitive to a wide variety of helix-distorting DNA lesions, including UV-induced DNA damage, an intra-strand crosslink, and bulky adducts. These results reveal that cdc6-2, tfb3, rio1, and Saci_0951 are play more important roles in broad responses to helix-distorting DNA damage than in specific responses to UV irradiation.

Loop-loop interactions govern multiple steps in indole-3-glycerol phosphate synthase catalysis.

Substrate binding, product release, and likely chemical catalysis in the tryptophan biosynthetic enzyme indole-3-glycerol phosphate synthase (IGPS) are dependent on the structural dynamics of the ß1α1 active-site loop. Statistical coupling analysis and molecular dynamic simulations had previously indicated that covarying residues in the ß1α1 and ß2α2 loops, corresponding to Arg54 and Asn90, respectively, in the Sulfolobus sulfataricus enzyme (ssIGPS), are likely important for coordinating functional motions of these loops. To test this hypothesis, we characterized site mutants at these positions for changes in catalytic function, protein stability and structural dynamics for the thermophilic ssIGPS enzyme. Although there were only modest changes in the overall steady-state kinetic parameters, solvent viscosity and solvent deuterium kinetic isotope effects indicated that these amino acid substitutions change the identity of the rate-determining step across multiple temperatures. Surprisingly, the N90A substitution had a dramatic effect on the general acid/base catalysis of the dehydration step, as indicated by the loss of the descending limb in the pH rate profile, which we had previously assigned to Lys53 on the ß1α1 loop. These changes in enzyme function are accompanied with a quenching of ps-ns and µs-ms timescale motions in the ß1α1 loop as measured by nuclear magnetic resonance studies. Altogether, our studies provide structural, dynamic and functional rationales for the coevolution of residues on the ß1α1 and ß2α2 loops, and highlight the multiple roles that the ß1α1 loop plays in IGPS catalysis. Thus, substitution of covarying residues in the active-site ß1α1 and ß2α2 loops of indole-3-glycerol phosphate synthase results in functional, structural, and dynamic changes, highlighting the multiple roles that the ß1α1 loop plays in enzyme catalysis and the importance of regulating the structural dynamics of this loop through noncovalent interactions with nearby structural elements.

A high-throughput screening method to reengineer DNA polymerases for random mutagenesis.

A screening system for directed evolution of DNA polymerases employing a fluorescent Scorpion probe as a reporter has been developed. The screening system has been validated in a directed evolution experiment of a distributive polymerase from the Y-polymerase family (Dpo4 from Sulfolobus solfataricus) which was improved in elongation efficiency of consecutive mismatches. The engineering campaign yielded improved Dpo4 polymerase variants one of which was successfully benchmarked in a sequence saturation mutagenesis experiment especially with regard to the desirable consecutive transversion mutations ([2.5-fold increase in frequency relative to a reference library prepared with Dpo4 WT). The Scorpion probe screening system enables to reengineer polymerases with low processivity and fidelity, and no secondary activities (i.e. exonuclease activity or strand displacement activity) to match demands in diversity generation for directed protein evolution.

Genome-scale analysis of gene function in the hydrogenotrophic methanogenic archaeon Methanococcus maripaludis.

A comprehensive whole-genome analysis of gene function by transposon mutagenesis and deep sequencing methodology has been implemented successfully in a representative of the Archaea domain. Libraries of transposon mutants were generated for the hydrogenotrophic, methanogenic archaeon Methanococcus maripaludis S2 using a derivative of the Tn5 transposon. About 89,000 unique insertions were mapped to the genome, which allowed for the classification of 526 genes or about 30% of the genome as possibly essential or strongly advantageous for growth in rich medium. Many of these genes were homologous to eukaryotic genes that encode fundamental processes in replication, transcription, and translation, providing direct evidence for their importance in Archaea. Some genes classified as possibly essential were unique to the archaeal or methanococcal lineages, such as that encoding DNA polymerase PolD. In contrast, the archaeal homolog to the gene encoding DNA polymerase B was not essential for growth, a conclusion confirmed by construction of an independent deletion mutation. Thus PolD, and not PolB, likely plays a fundamental role in DNA replication in methanococci. Similarly, 121 hypothetical ORFs were classified as possibly essential and likely play fundamental roles in methanococcal information processing or metabolism that are not established outside this group of prokaryotes.

Genetic studies on the virus-like regions in the genome of hyperthermophilic archaeon, Thermococcus kodakarensis.

Four virus-like integrated elements (TKV1, TKV2, TKV3, and TKV4) have been found in the genome of hyperthermophilic archaeon, Thermococcus kodakarensis, but virus particle formation has not been observed in the culture of T. kodakarensis. As the result of growth property analyses, mutants lacking each of the four virus-like regions exhibited decrease in the cell concentration and/or less growth rates compared to growth of parental strain (KU216), when the T. kodakarensis strains were grown at 85 °C in nutrient-rich medium. These results indicated that the genes in virus-like regions stimulated the cell growth under the observed growth condition. As the result of transcriptome analyses, genes involved in amino acid, energy or nucleotide metabolisms, and transport systems were up- or down-regulated in the cells of mutant strains. Interestingly, a decrease in transcriptional levels of glutamine synthetase (TK1796) gene (Tk-glnA) was observed in the cells of four mutant strains. Growths of TKV1 disrupted strain and TKV4 disrupted strain have shown no difference compared with that of KU216 by the addition of glutamate or glutamine, and the result suggested that TKV1 and TKV4 contributed to supply of amino acids to the cell.

Expression and characterization of an extremely thermostable ß-glycosidase (mannosidase) from the hyperthermophilic archaeon Pyrococcus furiosus DSM3638.

Genomic analysis of the hyperthermophilic archaeon Pyrococcus furiosus revealed the presence of an open reading frame (ORF PF0356) similar to the enzymes in glycoside hydrolase family 1. This ß-glycosidase, designated PFTG (P. furiosus thermostable glycosidase), was cloned and expressed in Escherichia coli. The expressed enzyme was purified by heat treatment and Ni-NTA affinity chromatography. The gene was composed of 1,452 bp encoding 483 amino acids for a protein with a predicted molecular mass of 56,326 Da. The temperature and pH optima were 100°C and 5.0 in sodium citrate buffer, respectively. The substrate specificity of PFTG suggests that it possesses characteristics of both ß-galactosidase and ß-mannosidase activities. However, through kinetic studies by ITC (Isothermal Titration Colorimetry) which is very sensitive method for enzyme kinetics, PF0356 enzyme revealed the highest catalytic efficiency toward p-nitrophenyl-ß-d-mannopyranoside (3.02 k(cat)/K(m)) and mannobiose (4.32 k(cat)/K(m)). The enzyme showed transglycosylation and transgalactosylation activities toward cellobiose, lactose and mannooligosaccharides that could produce GOS (galactooligosaccharides) and MOS (maltooligosaccharides). This novel hyperthermostable ß-glycosidase may be useful for food and pharmaceutical applications.

The 2-oxoacid dehydrogenase multienzyme complex of Haloferax volcanii.

Those aerobic archaea whose genomes have been sequenced possess four adjacent genes that, by sequence comparisons with bacteria and eukarya, appear to encode the component enzymes of a 2-oxoacid dehydrogenase multienzyme complex. However, no catalytic activity of any such complex has ever been detected in the archaea. In Thermoplasma acidophilum, evidence has been presented that the heterologously expressed recombinant enzyme possesses activity with the branched chain 2-oxoacids and, to a lesser extent, with pyruvate. In the current paper, we demonstrate that in Haloferax volcanii the four genes are transcribed as an operon in vivo. However, no functional complex or individual enzyme, except for the dihydrolipoamide dehydrogenase component, could be detected in this halophile grown on a variety of carbon sources. Dihydrolipoamide dehydrogenase is present at low catalytic activities, the level of which is increased three to fourfold when Haloferax volcanii is grown on the branched-chain amino acids valine, leucine and isoleucine.

Phosphoenolpyruvate synthetase and pyruvate, phosphate dikinase of Thermoproteus tenax: key pieces in the puzzle of archaeal carbohydrate metabolism.

The interconversion of phosphoenolpyruvate and pyruvate represents an important control point of the Embden-Meyerhof-Parnas (EMP) pathway in Bacteria and Eucarya, but little is known about this site of regulation in Archaea. Here we report on the coexistence of phosphoenolpyruvate synthetase (PEPS) and the first described archaeal pyruvate, phosphate dikinase (PPDK), which, besides pyruvate kinase (PK), are involved in the catalysis of this reaction in the hyperthermophilic crenarchaeote Thermoproteus tenax. The genes encoding T. tenax PEPS and PPDK were cloned and expressed in Escherichia coli, and the enzymic and regulatory properties of the recombinant gene products were analysed. Whereas PEPS catalyses the unidirectional conversion of pyruvate to phosphoenolpyruvate, PPDK shows a bidirectional activity with a preference for the catabolic reaction. In contrast to PK of T. tenax, which is regulated on transcript level but exhibits only limited regulatory potential on protein level, PEPS and PPDK activities are modulated by adenosine phosphates and intermediates of the carbohydrate metabolism. Additionally, expression of PEPS is regulated on transcript level in response to the offered carbon source as revealed by Northern blot analyses. The combined action of the differently regulated enzymes PEPS, PPDK and PK represents a novel way of controlling the interconversion of phosphoenolpyruvate and pyruvate in the reversible EMP pathway, allowing short-term and long-term adaptation to different trophic conditions. Comparative genomic analyses indicate the coexistence of PEPS, PPDK and PK in other Archaea as well, suggesting a similar regulation of the carbohydrate metabolism in these organisms.

Duplicated genes evolve slower than singletons despite the initial rate increase.

BACKGROUND: Gene duplication is an important mechanism that can lead to the emergence of new functions during evolution. The impact of duplication on the mode of gene evolution has been the subject of several theoretical and empirical comparative-genomic studies. It has been shown that, shortly after the duplication, genes seem to experience a considerable relaxation of purifying selection. RESULTS: Here we demonstrate two opposite effects of gene duplication on evolutionary rates. Sequence comparisons between paralogs show that, in accord with previous observations, a substantial acceleration in the evolution of paralogs occurs after duplication, presumably due to relaxation of purifying selection. The effect of gene duplication on evolutionary rate was also assessed by sequence comparison between orthologs that have paralogs (duplicates) and those that do not (singletons). It is shown that, in eukaryotes, duplicates, on average, evolve significantly slower than singletons. Eukaryotic ortholog evolutionary rates for duplicates are also negatively correlated with the number of paralogs per gene and the strength of selection between paralogs. A tally of annotated gene functions shows that duplicates tend to be enriched for proteins with known functions, particularly those involved in signaling and related cellular processes; by contrast, singletons include an over-abundance of poorly characterized proteins. CONCLUSIONS: These results suggest that whether or not a gene duplicate is retained by selection depends critically on the pre-existing functional utility of the protein encoded by the ancestral singleton. Duplicates of genes of a higher biological import, which are subject to strong functional constraints on the sequence, are retained relatively more often. Thus, the evolutionary trajectory of duplicated genes appears to be determined by two opposing trends, namely, the post-duplication rate acceleration and the generally slow evolutionary rate owing to the high level of functional constraints.

Arsenic resistance in Halobacterium sp. strain NRC-1 examined by using an improved gene knockout system.

The genome sequence of Halobacterium sp. strain NRC-1 encodes genes homologous to those responsible for conferring resistance to arsenic. These genes occur on both the large extrachromosomal replicon pNRC100 (arsADRC and arsR2M) and on the chromosome (arsB). We studied the role of these ars genes in arsenic resistance genetically by construction of gene knockouts. Deletion of the arsADRC gene cluster in a Halobacterium NRC-1 Deltaura3 strain resulted in increased sensitivity to arsenite and antimonite but not arsenate. In contrast, knockout of the chromosomal arsB gene did not show significantly increased sensitivity to arsenite or arsenate. We also found that knockout of the arsM gene produced sensitivity to arsenite, suggesting a second novel mechanism of arsenic resistance involving a putative arsenite(III)-methyltransferase. These results indicate that Halobacterium sp. strain NRC-1 contains an arsenite and antimonite extrusion system with significant differences from bacterial counterparts. Deletion analysis was facilitated by an improved method for gene knockouts/replacements in Halobacterium that relies on both selection and counterselection of ura3 using a uracil dropout medium and 5-fluoroorotic acid. The arsenite and antimonite resistance elements were shown to be regulated, with resistance to arsenic in the wild type inducible by exposure to a sublethal concentration of the metal. Northern hybridization and reverse transcription-PCR analyses showed that arsA, arsD, arsR, arsM, arsC, and arsB, but not arsR2, are inducible by arsenite and antimonite. We discuss novel aspects of arsenic resistance in this halophilic archaeon and technical improvements in our capability for gene knockouts in the genome.

Quantifying modularity in the evolution of biomolecular systems.

Functional modules are considered the primary building blocks of biomolecular systems. Here we study to what extent functional modules behave cohesively across genomes:That is, are functional modules also evolutionary modules? We probe this question by analyzing for a large collection of functional modules the phyletic patterns of their genes across 110 genomes. The majority of functional modules display limited evolutionary modularity. This result confirms certain comparative genome analyses, but is in contrast to implicit assumptions in the systems analysis of functional genomics data. We show that this apparent interspecies flexibility in the organization of functional modules depends more on functional differentiation within orthologous groups of genes, than on noise in the functional module definitions. When filtering out these sources of nonmodularity, even though very few functional modules behave perfectly modular in evolution, about half behave at least significantly more modular than a random set of genes. There are substantial differences in the evolutionary modularity between individual functional modules as well as between collections of functional modules, partly corresponding to conceptual differences in the functional module definition, which make comparisons between functional module collections biologically difficult to interpret. Analysis within one collection does not suffer from such differences, and we show that within the EcoCyc metabolic pathway database, biosynthetic pathways are evolutionarily more modular than catabolic pathways.

Analysis of two large functionally uncharacterized regions in the Methanopyrus kandleri AV19 genome.

BACKGROUND: For most sequenced prokaryotic genomes, about a third of the protein coding genes annotated are "orphan proteins", that is, they lack homology to known proteins. These hypothetical genes are typically short and randomly scattered throughout the genome. This trend is seen for most of the bacterial and archaeal genomes published to date. RESULTS: In contrast we have found that a large fraction of the genes coding for such orphan proteins in the Methanopyrus kandleri AV19 genome occur within two large regions. These genes have no known homologs except from other M. kandleri genes. However, analysis of their lengths, codon usage, and Ribosomal Binding Site (RBS) sequences shows that they are most likely true protein coding genes and not random open reading frames. CONCLUSIONS: Although these regions can be considered as candidates for massive lateral gene transfer, our bioinformatics analysis suggests that this is not the case. We predict many of the organism specific proteins to be transmembrane and belong to protein families that are non-randomly distributed between the regions. Consistent with this, we suggest that the two regions are most likely unrelated, and that they may be integrated plasmids.

Re-annotation of genome microbial coding-sequences: finding new genes and inaccurately annotated genes.

BACKGROUND: Analysis of any newly sequenced bacterial genome starts with the identification of protein-coding genes. Despite the accumulation of multiple complete genome sequences, which provide useful comparisons with close relatives among other organisms during the annotation process, accurate gene prediction remains quite difficult. A major reason for this situation is that genes are tightly packed in prokaryotes, resulting in frequent overlap. Thus, detection of translation initiation sites and/or selection of the correct coding regions remain difficult unless appropriate biological knowledge (about the structure of a gene) is imbedded in the approach. RESULTS: We have developed a new program that automatically identifies biologically significant candidate genes in a bacterial genome. Twenty-six complete prokaryotic genomes were analyzed using this tool, and the accuracy of gene finding was assessed by comparison with existing annotations. This analysis revealed that, despite the enormous effort of genome program annotators, a small but not negligible number of genes annotated within the framework of sequencing projects are likely to be partially inaccurate or plainly wrong. Moreover, the analysis of several putative new genes shows that, as expected, many short genes have escaped annotation. In most cases, these new genes revealed frameshifts that could be either artifacts or genuine frameshifts. Some entirely unexpected new genes have also been identified. This allowed us to get a more complete picture of prokaryotic genomes. The results of this procedure are progressively integrated into the SWISS-PROT reference databank. CONCLUSIONS: The results described in the present study show that our procedure is very satisfactory in terms of gene finding accuracy. Except in few cases, discrepancies between our results and annotations provided by individual authors can be accounted for by the nature of each annotation process or by specific characteristics of some genomes. This stresses that close cooperation between scientists, regular update and curation of the findings in databases are clearly required to reduce the level of errors in genome annotation (and hence in reducing the unfortunate spreading of errors through centralized data libraries).

Genome alignment, evolution of prokaryotic genome organization, and prediction of gene function using genomic context.

Gene order in prokaryotes is conserved to a much lesser extent than protein sequences. Only several operons, primarily those that code for physically interacting proteins, are conserved in all or most of the bacterial and archaeal genomes. Nevertheless, even the limited conservation of operon organization that is observed can provide valuable evolutionary and functional clues through multiple genome comparisons. A program for constructing gapped local alignments of conserved gene strings in two genomes was developed. The statistical significance of the local alignments was assessed using Monte Carlo simulations. Sets of local alignments were generated for all pairs of completely sequenced bacterial and archaeal genomes, and for each genome a template-anchored multiple alignment was constructed. In most pairwise genome comparisons, <10% of the genes in each genome belonged to conserved gene strings. When closely related pairs of species (i.e., two mycoplasmas) are excluded, the total coverage of genomes by conserved gene strings ranged from <5% for the cyanobacterium Synechocystis sp to 24% for the minimal genome of Mycoplasma genitalium, and 23% in Thermotoga maritima. The coverage of the archaeal genomes was only slightly lower than that of bacterial genomes. The majority of the conserved gene strings are known operons, with the ribosomal superoperon being the top-scoring string in most genome comparisons. However, in some of the bacterial-archaeal pairs, the superoperon is rearranged to the extent that other operons, primarily those subject to horizontal transfer, show the greatest level of conservation, such as the archaeal-type H+-ATPase operon or ABC-type transport cassettes. The level of gene order conservation among prokaryotic genomes was compared to the cooccurrence of genomes in clusters of orthologous genes (COGs) and to the conservation of protein sequences themselves. Only limited correlation was observed between these evolutionary variables. Gene order conservation shows a much lower variance than the cooccurrence of genomes in COGs, which indicates that intragenome homogenization via recombination occurs in evolution much faster than intergenome homogenization via horizontal gene transfer and lineage-specific gene loss. The potential of using template-anchored multiple-genome alignments for predicting functions of uncharacterized genes was quantitatively assessed. Functions were predicted or significantly clarified for approximately 90 COGs (approximately 4% of the total of 2414 analyzed COGs). The most significant predictions were obtained for the poorly characterized archaeal genomes; these include a previously uncharacterized restriction-modification system, a nuclease-helicase combination implicated in DNA repair, and the probable archaeal counterpart of the eukaryotic exosome. Multiple genome alignments are a resource for studies on operon rearrangement and disruption, which is central to our understanding of the evolution of prokaryotic genomes. Because of the rapid evolution of the gene order, the potential of genome alignment for prediction of gene functions is limited, but nevertheless, such predictions information significantly complements the results obtained through protein sequence and structure analysis.

Eight of fourteen gvp genes are sufficient for formation of gas vesicles in halophilic archaea.

The minimal number of genes required for the formation of gas vesicles in halophilic archaea has been determined. Single genes of the 14 gvp genes present in the p-vac region on plasmid pHH1 of Halobacterium salinarum (p-gvpACNO and p-gvpDEFGHIJKLM) were deleted, and the remaining genes were tested for the formation of gas vesicles in Haloferax volcanii transformants. The deletion of six gvp genes (p-gvpCN, p-gvpDE, and p-gvpHI) still enabled the production of gas vesicles in H. volcanii. The gas vesicles formed in some of these gvp gene deletion transformants were altered in shape (Delta I, Delta C) or strength (Delta H) but still functioned as flotation devices. A minimal p-vac region (minvac) containing the eight remaining genes (gvpFGJKLM-gvpAO) was constructed and tested for gas vesicle formation in H. volcanii. The minvac transformants did not form gas vesicles; however, minvac/gvpJKLM double transformants contained gas vesicles seen as light refractile bodies by phase-contrast microscopy. Transcript analyses demonstrated that minvac transformants synthesized regular amounts of gvpA mRNA, but the transcripts derived from gvpFGJKLM were mainly short and encompassed only gvpFG(J), suggesting that the gvpJKLM genes were not sufficiently expressed. Since gvpAO and gvpFGJKLM are the only gvp genes present in minvac/JKLM transformants containing gas vesicles, these gvp genes represent the minimal set required for gas vesicle formation in halophilic archaea. Homologs of six of these gvp genes are found in Anabaena flos-aquae, and homologs of all eight minimal halobacterial gvp genes are present in Bacillus megaterium and in the genome of Streptomyces coelicolor.