Chromosome architecture in an archaeal species naturally lacking structural maintenance of chromosomes proteins

Proteins in the structural maintenance of chromosomes (SMC) superfamily play key roles in chromosome organization and are ubiquitous across all domains of life. However, SMC proteins are notably absent in the Desulfurococcales of phylum Crenarchaeota. Intrigued by this observation, we performed chromosome conformation capture experiments in the model Desulfurococcales species Aeropyrum pernix. As in other archaea, we observe chromosomal interaction domains across the chromosome. The boundaries between chromosomal interaction domains show a dependence on transcription and translation for their definition. Importantly, however, we reveal an additional higher-order, bipartite organization of the chromosome-with a small high-gene-expression and self-interacting domain that is defined by transcriptional activity and loop structures. Viewing these data in the context of the distribution of SMC superfamily proteins in the Crenarchaeota, we suggest that the organization of the Aeropyrum genome represents an evolutionary antecedent of the compartmentalized architecture observed in the Sulfolobus lineage.

Identification of Up47 in three thermophilic archaea, one mesophilic archaeon, and one hyperthermophilic bacterium

Analysis of the profile of the tRNA modifications in several Archaea allowed us to observe a novel modified uridine in the V-loop of several tRNAs from two species: Pyrococcus furiosus and Sulfolobus acidocaldarius Recently, Ohira and colleagues characterized 2'-phosphouridine (Up) at position 47 in tRNAs of thermophilic Sulfurisphaera tokodaii, as well as in several other archaea and thermophilic bacteria. From the presence of the gene arkI corresponding to the RNA kinase responsible for Up47 formation, they also concluded that Up47 should be present in tRNAs of other thermophilic Archaea Reanalysis of our earlier data confirms that the unidentified residue in tRNAs of both P. furiosus and S. acidocaldarius is indeed 2'-phosphouridine followed by m5C48. Moreover, we find this modification in several tRNAs of other Archaea and of the hyperthermophilic bacterium Aquifex aeolicus.

Inactivation of Target RNA Cleavage of a III-B CRISPR-Cas System Induces Robust Autoimmunity in Saccharolobus islandicus

Type III CRISPR-Cas systems show the target (tg)RNA-activated indiscriminate DNA cleavage and synthesis of oligoadenylates (cOA) and a secondary signal that activates downstream nuclease effectors to exert indiscriminate RNA/DNA cleavage, and both activities are regulated in a spatiotemporal fashion. In III-B Cmr systems, cognate tgRNAs activate the two Cmr2-based activities, which are then inactivated via tgRNA cleavage by Cmr4, but how Cmr4 nuclease regulates the Cmr immunization remains to be experimentally characterized. Here, we conducted mutagenesis of Cmr4 conserved amino acids in Saccharolobus islandicus, and this revealed that Cmr4α RNase-dead (dCmr4α) mutation yields cell dormancy/death. We also found that plasmid-borne expression of dCmr4α in the wild-type strain strongly reduced plasmid transformation efficiency, and deletion of CRISPR arrays in the host genome reversed the dCmr4α inhibition. Expression of dCmr4α also strongly inhibited plasmid transformation with Cmr2α^HD and Cmr2α^Palm mutants, but the inhibition was diminished in Cmr2α^HD,Palm. Since dCmr4α-containing effectors lack spatiotemporal regulation, this allows an everlasting interaction between crRNA and cellular RNAs to occur. As a result, some cellular RNAs, which are not effective in mediating immunity due to the presence of spatiotemporal regulation, trigger autoimmunity of the Cmr-α system in the S. islandicus cells expressing dCmr4α. Together, these results pinpoint the crucial importance of tgRNA cleavage in autoimmunity avoidance and in the regulation of immunization of type III systems.

A Novel Family of Winged-Helix Single-Stranded DNA-Binding Proteins from Archaea

The winged helix superfamily comprises a large number of structurally related nucleic acid-binding proteins. While these proteins are often shown to bind dsDNA, few are known to bind ssDNA. Here, we report the identification and characterization of Sul7s, a novel winged-helix single-stranded DNA binding protein family highly conserved in Sulfolobaceae. Sul7s from Sulfolobus islandicus binds ssDNA with an affinity approximately 15-fold higher than that for dsDNA in vitro. It prefers binding oligo(dT)₃₀ over oligo(dC)₃₀ or a dG-rich 30-nt oligonucleotide, and barely binds oligo(dA)₃₀. Further, binding by Sul7s inhibits DNA strand annealing, but shows little effect on the melting temperature of DNA duplexes. The solution structure of Sul7s determined by NMR shows a winged helix-turn-helix fold, consisting of three α-helices, three β-strands, and two short wings. It interacts with ssDNA via a large positively charged binding surface, presumably resulting in ssDNA deformation. Our results shed significant light on not only non-OB fold single-stranded DNA binding proteins in Archaea, but also the divergence of the winged-helix proteins in both function and structure during evolution.

Intraspecific antagonism through viral toxin encoded by chronic Sulfolobus spindle-shaped virus

Virus-host interactions evolve along a symbiosis continuum from antagonism to mutualism. Long-term associations between virus and host, such as those in chronic infection, will select for traits that drive the interaction towards mutualism, especially when susceptible hosts are rare in the population. Virus-host mutualism has been demonstrated in thermophilic archaeal populations where Sulfolobus spindle-shaped viruses (SSVs) provide a competitive advantage to their host Sulfolobus islandicus by producing a toxin that kills uninfected strains. Here, we determine the genetic basis of this killing phenotype by identifying highly transcribed genes in cells that are chronically infected with a diversity of SSVs. We demonstrate that these genes alone confer growth inhibition by being expressed in uninfected cells via a Sulfolobus expression plasmid. Challenge of chronically infected strains with vector-expressed toxins revealed a nested network of cross-toxicity among divergent SSVs, with both broad and specific toxin efficacies. This suggests that competition between viruses and/or their hosts could maintain toxin diversity. We propose that competitive interactions among chronic viruses to promote their host fitness form the basis of virus-host mutualism. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

A Rapid Targeted Gene Inactivation Approach in Sulfolobus islandicus

Homologous recombination-based gene targeting is a powerful and classic reverse genetics approach to precisely elucidate in vivo gene functions in the organisms across all three domains of life. Gene function studies in Archaea, particularly for those flourishing in inhospitable natural environments that are anaerobic, usually hot, and acidic, have been a great challenge; however, this situation was recently overturned with the increasing availability of genetic manipulation systems in several cultivable archaeal species. In the present chapter, we describe a detailed procedure to rapidly generate gene disruption mutants in the hyperthermophilic crenarchaeon Sulfolobus islandicus via a recently developed Microhomology-Mediated Gene Inactivation (MMGI) approach. We highlight crucial experimental details required to be carefully considered when using the MMGI approach for genetic manipulations. We hope this highly reproducible procedure can supplement existing genetic tools for studying the biology of archaeal order Sulfolobales.

Transposon Insertion Mutagenesis in Hyperthermophilic Crenarchaeon Sulfolobus islandicus

Transposon insertion mutagenesis is a forward genetic approach that has been widely utilized for genetic characterization of bacteria and single-celled eukaryotes, and its applications are being rapidly expanded into a few archaeal model organisms for gene function analysis. Previously, we developed a Tn5-based in vivo transposon insertion mutagenesis system in the hyperthermophilic crenarchaeon S. islandicucs M.16.4 and defined the essential gene set under laboratory growth conditions. In this chapter, we will mainly focus on presenting details regarding the generation of a near-saturating transposon insertion mutant library in this crenarchaeal model. We envision that the traditional transposon-based forward mutagenesis screening paired with next generation sequencing will greatly speed up the exploration of archaeal genomic features.

Characterization of an archaeal inorganic pyrophosphatase from Sulfolobus islandicus using a [31P]-NMR-based assay

Inorganic pyrophosphatase catalyzes the conversion of pyrophosphate to phosphate and is often critical for driving reactions forward in cellular processes such as nucleic acid and protein synthesis. Commonly used methods for quantifying pyrophosphatase enzyme activity employ reacting liberated phosphate with a second molecule to produce absorbance changes or employing a second enzyme in coupled reactions to produce a product with a detectable absorbance. In this investigation, a novel [³¹P]-NMR spectroscopy-based assay was used to quantitatively measure the formation of phosphate and evaluate the activity of inorganic pyrophosphatase from the thermoacidophilic Crenarchaeota Sulfolobus islandicus. The enzymatic activity was directly measured via integration of the [³¹P] resonance associated with the phosphate product (δ = 2.1 ppm). Sulfolobus islandicus inorganic pyrophosphatase preferentially utilized Mg²⁺ as divalent cation and had pH and temperature optimums of 6.0 of 50 °C, respectively. The V_max value was 850 μmol/min/mg and the K_m for pyrophosphate was 1.02 mM. Sequence analysis indicates the enzyme is a Family I pyrophosphatase. Sulfolobus islandicus inorganic pyrophosphatase was shown to be inhibited by sodium fluoride with a IC₅₀ of 2.26 mM, compared to a IC₅₀ of 0.066 mM for yeast inorganic pyrophosphatase. These studies reveal that a [³¹P]-NMR spectroscopy-based assay is an effective method for analyzing catalysis by phosphate-producing enzymes.

High-resolution analysis of chromosome conformation in hyperthermophilic archaea

Chromosome conformation capture (3C) techniques are emerging as promising approaches to study genome organization in Archaea, the least understood domain of life in terms of chromosome biology. Here, we describe a 3C technique combined with deep sequencing for the hyperthermophilic archaeal genus Sulfolobus. Instead of using restriction enzymes compatible with fill-in labeling, this protocol uses the 4-bp blunt cutter AluI to generate high-resolution (up to 2 kb) contact maps from Sulfolobus species. For complete details on the use and execution of this protocol, please refer to Takemata and Bell (2021).

Structural basis for inhibition of an archaeal CRISPR-Cas type I-D large subunit by an anti-CRISPR protein

A hallmark of type I CRISPR-Cas systems is the presence of Cas3, which contains both the nuclease and helicase activities required for DNA cleavage during interference. In subtype I-D systems, however, the histidine-aspartate (HD) nuclease domain is encoded as part of a Cas10-like large effector complex subunit and the helicase activity in a separate Cas3' subunit, but the functional and mechanistic consequences of this organisation are not currently understood. Here we show that the Sulfolobus islandicus type I-D Cas10d large subunit exhibits an unusual domain architecture consisting of a Cas3-like HD nuclease domain fused to a degenerate polymerase fold and a C-terminal domain structurally similar to Cas11. Crystal structures of Cas10d both in isolation and bound to S. islandicus rod-shaped virus 3 AcrID1 reveal that the anti-CRISPR protein sequesters the large subunit in a non-functional state unable to form a cleavage-competent effector complex. The architecture of Cas10d suggests that the type I-D effector complex is similar to those found in type III CRISPR-Cas systems and that this feature is specifically exploited by phages for anti-CRISPR defence.

Archaeal Chromatin Proteins Cren7 and Sul7d Compact DNA by Bending and Bridging

Archaeal chromatin proteins Cren7 and Sul7d from Sulfolobus are DNA benders. To better understand their architectural roles in chromosomal DNA organization, we analyzed DNA compaction by Cren7 and Sis7d, a Sul7d family member, from Sulfolobus islandicus at the single-molecule (SM) level by total single-molecule internal reflection fluorescence microscopy (SM-TIRFM) and atomic force microscopy (AFM). We show that both Cren7 and Sis7d were able to compact singly tethered λ DNA into a highly condensed structure in a three-step process and that Cren7 was over an order of magnitude more efficient than Sis7d in DNA compaction. The two proteins were similar in DNA bending kinetics but different in DNA condensation patterns. At saturating concentrations, Sis7d formed randomly distributed clusters whereas Cren7 generated a single and highly condensed core on plasmid DNA. This observation is consistent with the greater ability of Cren7 than of Sis7d to bridge DNA. Our results offer significant insights into the mechanism and kinetics of chromosomal DNA organization in Crenarchaea.IMPORTANCE A long-standing question is how chromosomal DNA is packaged in Crenarchaeota, a major group of archaea, which synthesize large amounts of unique small DNA-binding proteins but in general contain no archaeal histones. In the present work, we tested our hypothesis that the two well-studied crenarchaeal chromatin proteins Cren7 and Sul7d compact DNA by both DNA bending and bridging. We show that the two proteins are capable of compacting DNA, albeit with different efficiencies and in different manners, at the single molecule level. We demonstrate for the first time that the two proteins, which have long been regarded as DNA binders and benders, are able to mediate DNA bridging, and this previously unknown property of the proteins allows DNA to be packaged into highly condensed structures. Therefore, our results provide significant insights into the mechanism and kinetics of chromosomal DNA organization in Crenarchaeota.

Killer Archaea: Virus-Mediated Antagonism to CRISPR-Immune Populations Results in Emergent Virus-Host Mutualism

Theory, simulation, and experimental evolution demonstrate that diversified CRISPR-Cas immunity to lytic viruses can lead to stochastic virus extinction due to a limited number of susceptible hosts available to each potential new protospacer escape mutation. Under such conditions, theory predicts that to evade extinction, viruses evolve toward decreased virulence and promote vertical transmission and persistence in infected hosts. To better understand the evolution of host-virus interactions in microbial populations with active CRISPR-Cas immunity, we studied the interaction between CRISPR-immune Sulfolobus islandicus cells and immune-deficient strains that are infected by the chronic virus SSV9. We demonstrate that Sulfolobus islandicus cells infected with SSV9, and with other related SSVs, kill uninfected, immune strains through an antagonistic mechanism that is a protein and is independent of infectious virus. Cells that are infected with SSV9 are protected from killing and persist in the population. We hypothesize that this infection acts as a form of mutualism between the host and the virus by removing competitors in the population and ensuring continued vertical transmission of the virus within populations with diversified CRISPR-Cas immunity.IMPORTANCE Multiple studies, especially those focusing on the role of lytic viruses in key model systems, have shown the importance of viruses in shaping microbial populations. However, it has become increasingly clear that viruses with a long host-virus interaction, such as those with a chronic lifestyle, can be important drivers of evolution and have large impacts on host ecology. In this work, we describe one such interaction with the acidic crenarchaeon Sulfolobus islandicus and its chronic virus Sulfolobus spindle-shaped virus 9. Our work expands the view in which this symbiosis between host and virus evolved, describing a killing phenotype which we hypothesize has evolved in part due to the high prevalence and diversity of CRISPR-Cas immunity seen in natural populations. We explore the implications of this phenotype in population dynamics and host ecology, as well as the implications of mutualism between this virus-host pair.

Functional Insights Into Protein Acetylation in the Hyperthermophilic Archaeon Sulfolobus islandicus

Proteins undergo acetylation at the Nε-amino group of lysine residues and the Nα-amino group of the N terminus in Archaea as in Bacteria and Eukarya. However, the extent, pattern and roles of the modifications in Archaea remain poorly understood. Here we report the proteomic analyses of a wild-type Sulfolobus islandicus strain and its mutant derivative strains lacking either a homolog of the protein acetyltransferase Pat (ΔSisPat) or a homolog of the Nt-acetyltransferase Ard1 (ΔSisArd1). A total of 1708 Nε-acetylated lysine residues in 684 proteins (26% of the total proteins), and 158 Nt-acetylated proteins (44% of the identified proteins) were found in S. islandicus ΔSisArd1 grew more slowly than the parental strain, whereas ΔSisPat showed no significant growth defects. Only 24 out of the 1503 quantifiable Nε-acetylated lysine residues were differentially acetylated, and all but one of the 24 residues were less acetylated by >1.3 fold in ΔSisPat than in the parental strain, indicating the narrow substrate specificity of the enzyme. Six acyl-CoA synthetases were the preferred substrates of SisPat in vivo, suggesting that Nε-acetylation by the acetyltransferase is involved in maintaining metabolic balance in the cell. Acetylation of acyl-CoA synthetases by SisPat occurred at a sequence motif conserved among all three domains of life. On the other hand, 92% of the acetylated N termini identified were acetylated by SisArd1 in the cell. The enzyme exhibited broad substrate specificity and could modify nearly all types of the target N termini of human NatA-NatF. The deletion of the SisArd1 gene altered the cellular levels of 18% of the quantifiable proteins (1518) by >1.5 fold. Consistent with the growth phenotype of ΔSisArd1, the cellular levels of proteins involved in cell division and cell cycle control, DNA replication, and purine synthesis were significantly lowered in the mutant than those in the parental strain.

Improvement of ST0452 N-Acetylglucosamine-1-Phosphate Uridyltransferase Activity by the Cooperative Effect of Two Single Mutations Identified through Structure-Based Protein Engineering

We showed previously that the Y97N mutant of the ST0452 protein, isolated from Sulfolobus tokodaii, exhibited over 4 times higher N-acetylglucosamine-1-phosphate (GlcNAc-1-P) uridyltransferase (UTase) activity, compared with that of the wild-type ST0452 protein. We determined the three-dimensional structure of the Y97N protein to explore the detailed mechanism underlying this increased activity. The overall structure was almost identical to that of the wild-type ST0452 protein (PDB ID 2GGO), with residue 97 (Asn) interacting with the O-5 atom of N-acetylglucosamine (GlcNAc) in the complex without metal ions. The same interaction was observed for Escherichia coli GlmU in the absence of metal ions. These observations indicated that the three-dimensional structure of the Y97N protein was not changed by this substitution but the interactions with the substrate were slightly modified, which might cause the activity to increase. The crystal structure of the Y97N protein also showed that positions 146 (Glu) and 80 (Thr) formed interactions with GlcNAc, and an engineering strategy was applied to these residues to increase activity. All proteins substituted at position 146 had drastically decreased activities, whereas several proteins substituted at position 80 showed higher GlcNAc-1-P UTase activity, compared to that of the wild-type protein. The substituted amino acids at positions 80 and 97 might result in optimized interactions with the substrate; therefore, we predicted that the combination of these two substitutions might cooperatively increase GlcNAc-1-P UTase activity. Of the four double mutant ST0452 proteins generated, T80S/Y97N showed 6.5-times-higher activity, compared to that of the wild-type ST0452 protein, revealing that these two substituted residues functioned cooperatively to increase GlcNAc-1-P UTase activity.IMPORTANCE We demonstrated that the enzymatic activity of a thermostable protein was over 4 times higher than that of the wild-type protein following substitution of a single amino acid, without affecting its thermostability. The three-dimensional structure of the improved mutant protein complexed with substrate was determined. The same overall structure and interaction between the substituted residue and the GlcNAc substrate as observed in the well-characterized bacterial enzyme suggested that the substitution of Tyr at position 97 by Asn might slightly change the interaction. This subtle change in the interaction might potentially increase the GlcNAc-1-P UTase activity of the mutant protein. These observations indicated that a drastic change in the structure of a natural thermostable enzyme is not necessary to increase its activity; a subtle change in the interaction with the substrate might be sufficient. Cooperative effects were observed in the appropriate double mutant protein. This work provides useful information for the future engineering of natural enzymes.

Exploration and isolation of novel thermophiles in frozen enrichment cultures derived from a terrestrial acidic hot spring

An isolation strategy, exploring novel microorganisms in frozen enrichment cultures (ENFE), which uses a combination of enrichment culture and 16S rRNA gene clone analysis, was evaluated for isolating uncultured thermophiles from a terrestrial acidic hot spring. The procedure comprised (a) multiple enrichment cultures under various conditions, (b) cryostorage of all enrichments, (c) microbial community analyses of the enrichments using 16S rRNA gene sequences, and (d) purification of microorganisms from enrichments containing previously uncultured microorganisms. The enrichments were performed under a total of 36 conditions, and 16 of these enrichments yielded positive microbial growth with the detection of three previously uncultured archaea. Two of the three previously uncultured archaea, strains HS-1 and HS-3, were successfully isolated. Strain HS-1 and HS-3 represented a novel lineage of the order Sulfolobales and novel species of the genus Sulfolobus, respectively. Although innovative isolation methods play strategic roles in isolating previously uncultured microorganisms, the ENFE strategy showed potential for characterizing and isolating such microorganisms using conventional media and techniques.

The CRISPR-associated Csx1 protein of Pyrococcus furiosus is an adenosine-specific endoribonuclease

Prokaryotes are frequently exposed to potentially harmful invasive nucleic acids from phages, plasmids, and transposons. One method of defense is the CRISPR-Cas adaptive immune system. Diverse CRISPR-Cas systems form distinct ribonucleoprotein effector complexes that target and cleave invasive nucleic acids to provide immunity. The Type III-B Cmr effector complex has been found to target the RNA and DNA of the invader in the various bacterial and archaeal organisms where it has been characterized. Interestingly, the gene encoding the Csx1 protein is frequently located in close proximity to the Cmr1-6 genes in many genomes, implicating a role for Csx1 in Cmr function. However, evidence suggests that Csx1 is not a stably associated component of the Cmr effector complex, but is necessary for DNA silencing by the Cmr system in Sulfolobus islandicus. To investigate the function of the Csx1 protein, we characterized the activity of recombinant Pyrococcus furiosus Csx1 against various nucleic acid substrates. We show that Csx1 is a metal-independent, endoribonuclease that acts selectively on single-stranded RNA and cleaves specifically after adenosines. The RNA cleavage activity of Csx1 is dependent upon a conserved HEPN motif located within the C-terminal domain of the protein. This motif is also key for activity in other known ribonucleases. Collectively, the findings indicate that invader silencing by Type III-B CRISPR-Cas systems relies both on RNA and DNA nuclease activities from the Cmr effector complex as well as on the affiliated, trans-acting Csx1 endoribonuclease.

Genetic analysis of the Holliday junction resolvases Hje and Hjc in Sulfolobus islandicus

The in vivo functions of Hje and Hjc, two Holliday junction resolvases in Sulfolobus islandicus were investigated. We found that deletion of either hje or hjc had no effect on normal cell growth, while deletion of both hje and hjc is lethal. Although Hjc is the conserved resolvase in all archaea, the hje deletion rather than hjc deletion rendered cells more sensitive to DNA-damaging agents such as hydroxyurea, cisplatin, and methyl methanesulfonate than the wild type (WT). Intriguingly, the sensitivity of Δhje could not be rescued by ectopic expression of Hje from a plasmid and Hje overexpression slowed growth and large cells appeared with more than two genome equivalents. We showed that Hje was maintained at a low level in WT cells. Furthermore, transcriptomic microarray analysis revealed that the abundance of transcripts of many genes including those involved in DNA replication, repair, transcription regulation, and cell division changed drastically in the Hje-overexpressed strain. However, only limited genes were up- or downregulated in the hje deletion strain. Our findings collectively suggest that Hje is the primary resolvase involved in DNA repair and its expression must be tightly controlled in cells.

An archaeal CRISPR type III-B system exhibiting distinctive RNA targeting features and mediating dual RNA and DNA interference

CRISPR-Cas systems provide a small RNA-based mechanism to defend against invasive genetic elements in archaea and bacteria. To investigate the in vivo mechanism of RNA interference by two type III-B systems (Cmr-α and Cmr-β) in Sulfolobus islandicus, a genetic assay was developed using plasmids carrying an artificial mini-CRISPR (AC) locus with a single spacer. After pAC plasmids were introduced into different strains, Northern analyses confirmed that mature crRNAs were produced from the plasmid-borne CRISPR loci, which then guided gene silencing to target gene expression. Spacer mutagenesis identified a trinucleotide sequence in the 3'-region of crRNA that was crucial for RNA interference. Studying mutants lacking Cmr-α or Cmr-β system showed that each Cmr complex exhibited RNA interference. Strikingly, these analyses further revealed that the two Cmr systems displayed distinctive interference features. Whereas Cmr-β complexes targeted transcripts and could be recycled in RNA cleavage, Cmr-α complexes probably targeted nascent RNA transcripts and remained associated with the substrate. Moreover, Cmr-β exhibited much stronger RNA cleavage activity than Cmr-α. Since we previously showed that S. islandicus Cmr-α mediated transcription-dependent DNA interference, the Cmr-α constitutes the first CRISPR system exhibiting dual targeting of RNA and DNA.

Comparative genome analysis reveals metabolic versatility and environmental adaptations of Sulfobacillus thermosulfidooxidans strain ST

The genus Sulfobacillus is a cohort of mildly thermophilic or thermotolerant acidophiles within the phylum Firmicutes and requires extremely acidic environments and hypersalinity for optimal growth. However, our understanding of them is still preliminary partly because few genome sequences are available. Here, the draft genome of Sulfobacillus thermosulfidooxidans strain ST was deciphered to obtain a comprehensive insight into the genetic content and to understand the cellular mechanisms necessary for its survival. Furthermore, the expressions of key genes related with iron and sulfur oxidation were verified by semi-quantitative RT-PCR analysis. The draft genome sequence of Sulfobacillus thermosulfidooxidans strain ST, which encodes 3225 predicted coding genes on a total length of 3,333,554 bp and a 48.35% G+C, revealed the high degree of heterogeneity with other Sulfobacillus species. The presence of numerous transposases, genomic islands and complete CRISPR/Cas defence systems testifies to its dynamic evolution consistent with the genome heterogeneity. As expected, S. thermosulfidooxidans encodes a suit of conserved enzymes required for the oxidation of inorganic sulfur compounds (ISCs). The model of sulfur oxidation in S. thermosulfidooxidans was proposed, which showed some different characteristics from the sulfur oxidation of Gram-negative A. ferrooxidans. Sulfur oxygenase reductase and heterodisulfide reductase were suggested to play important roles in the sulfur oxidation. Although the iron oxidation ability was observed, some key proteins cannot be identified in S. thermosulfidooxidans. Unexpectedly, a predicted sulfocyanin is proposed to transfer electrons in the iron oxidation. Furthermore, its carbon metabolism is rather flexible, can perform the transformation of pentose through the oxidative and non-oxidative pentose phosphate pathways and has the ability to take up small organic compounds. It encodes a multitude of heavy metal resistance systems to adapt the heavy metal-containing environments.

Identification and codon reading properties of 5-cyanomethyl uridine, a new modified nucleoside found in the anticodon wobble position of mutant haloarchaeal isoleucine tRNAs

Most archaea and bacteria use a modified C in the anticodon wobble position of isoleucine tRNA to base pair with A but not with G of the mRNA. This allows the tRNA to read the isoleucine codon AUA without also reading the methionine codon AUG. To understand why a modified C, and not U or modified U, is used to base pair with A, we mutated the C34 in the anticodon of Haloarcula marismortui isoleucine tRNA (tRNA2(Ile)) to U, expressed the mutant tRNA in Haloferax volcanii, and purified and analyzed the tRNA. Ribosome binding experiments show that although the wild-type tRNA2(Ile) binds exclusively to the isoleucine codon AUA, the mutant tRNA binds not only to AUA but also to AUU, another isoleucine codon, and to AUG, a methionine codon. The G34 to U mutant in the anticodon of another H. marismortui isoleucine tRNA species showed similar codon binding properties. Binding of the mutant tRNA to AUG could lead to misreading of the AUG codon and insertion of isoleucine in place of methionine. This result would explain why most archaea and bacteria do not normally use U or a modified U in the anticodon wobble position of isoleucine tRNA for reading the codon AUA. Biochemical and mass spectrometric analyses of the mutant tRNAs have led to the discovery of a new modified nucleoside, 5-cyanomethyl U in the anticodon wobble position of the mutant tRNAs. 5-Cyanomethyl U is present in total tRNAs from euryarchaea but not in crenarchaea, eubacteria, or eukaryotes.

Shining a light on dark sequencing: characterising errors in Ion Torrent PGM data

The Ion Torrent Personal Genome Machine (PGM) is a new sequencing platform that substantially differs from other sequencing technologies by measuring pH rather than light to detect polymerisation events. Using re-sequencing datasets, we comprehensively characterise the biases and errors introduced by the PGM at both the base and flow level, across a combination of factors, including chip density, sequencing kit, template species and machine. We found two distinct insertion/deletion (indel) error types that accounted for the majority of errors introduced by the PGM. The main error source was inaccurate flow-calls, which introduced indels at a raw rate of 2.84% (1.38% after quality clipping) using the OneTouch 200 bp kit. Inaccurate flow-calls typically resulted in over-called short-homopolymers and under-called long-homopolymers. Flow-call accuracy decreased with consecutive flow cycles, but we also found significant periodic fluctuations in the flow error-rate, corresponding to specific positions within the flow-cycle pattern. Another less common PGM error, high frequency indel (HFI) errors, are indels that occur at very high frequency in the reads relative to a given base position in the reference genome, but in the majority of instances were not replicated consistently across separate runs. HFI errors occur approximately once every thousand bases in the reference, and correspond to 0.06% of bases in reads. Currently, the PGM does not achieve the accuracy of competing light-based technologies. However, flow-call inaccuracy is systematic and the statistical models of flow-values developed here will enable PGM-specific bioinformatics approaches to be developed, which will account for these errors. HFI errors may prove more challenging to address, especially for polymorphism and amplicon applications, but may be overcome by sequencing the same DNA template across multiple chips.

DNA sequencing is used routinely within biology to reveal the genetic information of living organisms. In recent years, technological advances have led to the availability of high-throughput, low-cost DNA sequencing machines (‘sequencers’). In 2011, Life Sciences released a new sequencer, the Ion Torrent Personal Genome Machine (PGM). This is the first sequencer to measure changes in pH rather that emitted light to register sequencing reactions. Consequently, this unique technology is both cost-effective and advertised to have high accuracy, making it attractive for many laboratories. However, every sequencing technology introduces unique errors and biases into the resulting DNA sequences, and understanding PGM-specific characteristics is crucial to determining suitable applications for this new technology. We comprehensively examine the types of errors and biases in PGM-sequenced data across several experimental variables, including chip density, template kit, template DNA and across two machines. Using statistical approaches, we quantify the influence of experimental variables, as well as DNA sequence-specific effects, and find that the PGM has two types of technology-specific errors. We also find that the accuracy of the PGM is poorer than that of light-based technologies, and we make recommendations for this technology as well as provide statistical models for overcoming PGM sequencing errors.

PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes

New microbial genomes are constantly being sequenced, and it is crucial to accurately determine their taxonomic identities and evolutionary relationships. Here we report PhyloPhlAn, a new method to assign microbial phylogeny and putative taxonomy using >400 proteins optimized from among 3,737 genomes. This method measures the sequence diversity of all clades, classifies genomes from deep-branching candidate divisions through closely related subspecies and improves consistency between phylogenetic and taxonomic groupings. PhyloPhlAn improved taxonomic accuracy for existing and newly sequenced genomes, detecting 157 erroneous labels, correcting 46 and placing or refining 130 new genomes. We provide examples of accurate classifications from subspecies (Sulfolobus spp.) to phyla, and of preliminary rooting of deep-branching candidate divisions, including consistent statistical support for Caldiserica (formerly candidate division OP5). PhyloPhlAn will thus be useful for both phylogenetic assessment and taxonomic quality control of newly sequenced genomes. The final phylogenies, conserved protein sequences and open-source implementation are available online.

Molecular characterization of a thermostable aldehyde dehydrogenase (ALDH) from the hyperthermophilic archaeon Sulfolobus tokodaii strain 7

Aldehyde dehydrogenase (ALDH) is a widely distributed enzyme in nature. Although many ALDHs have been reported until now, the detailed enzymatic properties of ALDH from Archaea remain elusive. Herein, we describe the characterization of an ALDH from the hyperthermophilic archaeon Sulfolobus tokodaii. The enzyme (stALDH) could utilize various aldehydes as substrates, and maximal activity was found with acetaldehyde and the coenzyme NAD. The optimal temperature and pH were 80 °C and 8, respectively, and high thermostability was found with the half-life at 90 °C to be 4 h. The enzyme was considerably resistant to nitroglycerin (GTN) inhibition, which could be restored by reducing agent DTT or (±)-α-lipoic acid. Coenzyme NAD or NADP could regulate the enzymatic thermostability, as well as the esterase activity. Molecular modeling suggested that the enzyme harbored similar structural arrangement with its eukaryotic and bacterial counterparts. Sequence alignment showed the conserved catalytic residues E240 and C274 and cofactor interactive sites N142, K165, I168 and E370, the function of which were verified by site-directed mutagenesis analysis. This is the most thermostable ALDH reported until now and the unique property of this enzyme is potentially beneficial in the fields of biotechnology and biomedicine.

Structural and enzymatic characterization of the lactonase SisLac from Sulfolobus islandicus

Background

A new member of the Phosphotriesterase-Like Lactonases (PLL) family from the hyperthermophilic archeon Sulfolobus islandicus (SisLac) has been characterized. SisLac is a native lactonase that exhibits a high promiscuous phosphotriesterase activity. SisLac thus represents a promising target for engineering studies, exhibiting both detoxification and bacterial quorum quenching abilities, including human pathogens such as Pseudomonas aeruginosa.

Methodology/Principal Findings

Here, we describe the substrate specificity of SisLac, providing extensive kinetic studies performed with various phosphotriesters, esters, N-acyl-homoserine lactones (AHLs) and other lactones as substrates. Moreover, we solved the X-ray structure of SisLac and structural comparisons with the closely related SsoPox structure highlighted differences in the surface salt bridge network and the dimerization interface. SisLac and SsoPox being close homologues (91% sequence identity), we undertook a mutational study to decipher these structural differences and their putative consequences on the stability and the catalytic properties of these proteins.

Conclusions/Significance

We show that SisLac is a very proficient lactonase against aroma lactones and AHLs as substrates. Hence, data herein emphasize the potential role of SisLac as quorum quenching agent in Sulfolobus. Moreover, despite the very high sequence homology with SsoPox, we highlight key epistatic substitutions that influence the enzyme stability and activity.

In silico classification of proteins from acidic and neutral cytoplasms

Protein acidostability is a common problem in biopharmaceutical and other industries. However, it remains a great challenge to engineer proteins for enhanced acidostability because our knowledge of protein acidostabilization is still very limited. In this paper, we present a comparative study of proteins from bacteria with acidic (AP) and neutral cytoplasms (NP) using an integrated statistical and machine learning approach. We construct a set of 393 non-redundant AP-NP ortholog pairs and calculate a total of 889 sequence based features for these proteins. The pairwise alignments of these ortholog pairs are used to build a residue substitution propensity matrix between APs and NPs. We use Gini importance provided by the Random Forest algorithm to rank the relative importance of these features. A scoring function using the 10 most significant features is developed and optimized using a hill climbing algorithm. The accuracy of the score function is 86.01% in predicting AP-NP ortholog pairs and is 76.65% in predicting non-ortholog AP-NP pairs, suggesting that there are significant differences between APs and NPs which can be used to predict relative acidostability of proteins. The overall trends uncovered in the study can be used as general guidelines for designing acidostable proteins. To best of our knowledge, this work represents the first systematic comparative study of the acidostable proteins and their non-acidostable orthologs.

Modulation of CRISPR locus transcription by the repeat-binding protein Cbp1 in Sulfolobus

CRISPR loci are essential components of the adaptive immune system of archaea and bacteria. They consist of long arrays of repeats separated by DNA spacers encoding guide RNAs (crRNA), which target foreign genetic elements. Cbp1 (CRISPR DNA repeat binding protein) binds specifically to the multiple direct repeats of CRISPR loci of members of the acidothermophilic, crenarchaeal order Sulfolobales. cbp1 gene deletion from Sulfolobus islandicus REY15A produced a strong reduction in pre-crRNA yields from CRISPR loci but did not inhibit the foreign DNA targeting capacity of the CRISPR/Cas system. Conversely, overexpression of Cbp1 in S. islandicus generated an increase in pre-crRNA yields while the level of reverse strand transcripts from CRISPR loci remained unchanged. It is proposed that Cbp1 modulates production of longer pre-crRNA transcripts from CRISPR loci. A possible mechanism is that it minimizes interference from potential transcriptional signals carried on spacers deriving from A-T-rich genetic elements and, occasionally, on DNA repeats. Supporting evidence is provided by microarray and northern blotting analyses, and publicly available whole-transcriptome data for S. solfataricus P2.

Patterns of gene flow define species of thermophilic Archaea

Despite a growing appreciation of their vast diversity in nature, mechanisms of speciation are poorly understood in Bacteria and Archaea. Here we use high-throughput genome sequencing to identify ongoing speciation in the thermoacidophilic Archaeon Sulfolobus islandicus. Patterns of homologous gene flow among genomes of 12 strains from a single hot spring in Kamchatka, Russia, demonstrate higher levels of gene flow within than between two persistent, coexisting groups, demonstrating that these microorganisms fit the biological species concept. Furthermore, rates of gene flow between two species are decreasing over time in a manner consistent with incipient speciation. Unlike other microorganisms investigated, we do not observe a relationship between genetic divergence and frequency of recombination along a chromosome, or other physical mechanisms that would reduce gene flow between lineages. Each species has its own genetic island encoding unique physiological functions and a unique growth phenotype that may be indicative of ecological specialization. Genetic differentiation between these coexisting groups occurs in large genomic "continents," indicating the topology of genomic divergence during speciation is not uniform and is not associated with a single locus under strong diversifying selection. These data support a model where species do not require physical barriers to gene flow but are maintained by ecological differentiation.

Functional promiscuity of the COG0720 family

The biosynthesis of GTP derived metabolites such as tetrahydrofolate (THF), biopterin (BH(4)), and the modified tRNA nucleosides queuosine (Q) and archaeosine (G(+)) relies on several enzymes of the Tunnel-fold superfamily. A subset of these proteins includes the 6-pyruvoyltetrahydropterin (PTPS-II), PTPS-III, and PTPS-I homologues, all members of the COG0720 family that have been previously shown to transform 7,8-dihydroneopterin triphosphate (H(2)NTP) into different products. PTPS-II catalyzes the formation of 6-pyruvoyltetrahydropterin in the BH(4) pathway, PTPS-III catalyzes the formation of 6-hydroxylmethyl-7,8-dihydropterin in the THF pathway, and PTPS-I catalyzes the formation of 6-carboxy-5,6,7,8-tetrahydropterin in the Q pathway. Genes of these three enzyme families are often misannotated as they are difficult to differentiate by sequence similarity alone. Using a combination of physical clustering, signature motif, phylogenetic codistribution analyses, in vivo complementation studies, and in vitro enzymatic assays, a complete reannotation of the COG0720 family was performed in prokaryotes. Notably, this work identified and experimentally validated dual function PTPS-I/III enzymes involved in both THF and Q biosynthesis. Both in vivo and in vitro analyses showed that the PTPS-I family could tolerate a translation of the active site cysteine and was inherently promiscuous, catalyzing different reactions on the same substrate or the same reaction on different substrates. Finally, the analysis and experimental validation of several archaeal COG0720 members confirmed the role of PTPS-I in archaeosine biosynthesis and resulted in the identification of PTPS-III enzymes with variant signature sequences in Sulfolobus species. This study reveals an expanded versatility of the COG0720 family members and illustrates that for certain protein families extensive comparative genomic analysis beyond homology is required to correctly predict function.

Deletion of the topoisomerase III gene in the hyperthermophilic archaeon Sulfolobus islandicus results in slow growth and defects in cell cycle control

Topoisomerase III (topo III), a type IA topoisomerase, is widespread in hyperthermophilic archaea. In order to interrogate the in vivo role of archaeal topo III, we constructed and characterized a topo III gene deletion mutant of Sulfolobus islandicus. The mutant was viable but grew more slowly than the wild-type strain, especially in a nutrient-poor medium. Flow cytometry analysis revealed changes of the mutant in growth cycle characteristics including an increase in proportion of cells containing either more than two genome equivalents or less than one genome equivalent in exponentially-growing cultures. As shown by fluorescence microscopy, a fraction of mutant cells in the cultures were drastically enlarged, and at least some of the enlarged cells were apparently capable of resuming cell division. The mutant also shows a different transcriptional profile from that of the wild-type strain. Our results suggest that the enzyme may serve roles in chromosomal segregation and control of the level of supercoiling in the cell.

Conserved residues in membrane-bound acid pyrophosphatase from Sulfolobus tokodaii, a thermoacidophilic archaeon

A membrane-intrinsic acid pyrophosphatase (ST2226) from Sulfolobus tokodaii, a thermoacidophilic archaeon, is possibly involved in glycoprotein biosynthesis and belongs to the phosphatidic acid phosphatase class 2 superfamily, including both membrane-intrinsic and soluble enzymes with divergent functions ranging from dephosphorylation of undecaprenylpyrophosphate and phospho-monoesters such as glucose-6-phosphate to vanadium-containing chloroperoxidation. ST2226 is an archaeal ortholog of these enzymes sharing a common phosphatase motif. Through site-directed mutagenesis as to each of the conserved residues, the catalytic roles of the latter were deduced, as well as the transmembrane topology with all the conserved residues in close proximity to the outside of the membrane.

Horizontal transfer, not duplication, drives the expansion of protein families in prokaryotes

Gene duplication followed by neo- or sub-functionalization deeply impacts the evolution of protein families and is regarded as the main source of adaptive functional novelty in eukaryotes. While there is ample evidence of adaptive gene duplication in prokaryotes, it is not clear whether duplication outweighs the contribution of horizontal gene transfer in the expansion of protein families. We analyzed closely related prokaryote strains or species with small genomes (Helicobacter, Neisseria, Streptococcus, Sulfolobus), average-sized genomes (Bacillus, Enterobacteriaceae), and large genomes (Pseudomonas, Bradyrhizobiaceae) to untangle the effects of duplication and horizontal transfer. After removing the effects of transposable elements and phages, we show that the vast majority of expansions of protein families are due to transfer, even among large genomes. Transferred genes--xenologs--persist longer in prokaryotic lineages possibly due to a higher/longer adaptive role. On the other hand, duplicated genes--paralogs--are expressed more, and, when persistent, they evolve slower. This suggests that gene transfer and gene duplication have very different roles in shaping the evolution of biological systems: transfer allows the acquisition of new functions and duplication leads to higher gene dosage. Accordingly, we show that paralogs share most protein-protein interactions and genetic regulators, whereas xenologs share very few of them. Prokaryotes invented most of life's biochemical diversity. Therefore, the study of the evolution of biology systems should explicitly account for the predominant role of horizontal gene transfer in the diversification of protein families.

CRISPR associated diversity within a population of Sulfolobus islandicus

BACKGROUND

Predator-prey models for virus-host interactions predict that viruses will cause oscillations of microbial host densities due to an arms race between resistance and virulence. A new form of microbial resistance, CRISPRs (clustered regularly interspaced short palindromic repeats) are a rapidly evolving, sequence-specific immunity mechanism in which a short piece of invading viral DNA is inserted into the host's chromosome, thereby rendering the host resistant to further infection. Few studies have linked this form of resistance to population dynamics in natural microbial populations.

METHODOLOGY/PRINCIPAL FINDINGS

We examined sequence diversity in 39 strains of the archeaon Sulfolobus islandicus from a single, isolated hot spring from Kamchatka, Russia to determine the effects of CRISPR immunity on microbial population dynamics. First, multiple housekeeping genetic markers identify a large clonal group of identical genotypes coexisting with a diverse set of rare genotypes. Second, the sequence-specific CRISPR spacer arrays split the large group of isolates into two very different groups and reveal extensive diversity and no evidence for dominance of a single clone within the population.

CONCLUSIONS/SIGNIFICANCE

The evenness of resistance genotypes found within this population of S. islandicus is indicative of a lack of strain dominance, in contrast to the prediction for a resistant strain in a simple predator-prey interaction. Based on evidence for the independent acquisition of resistant sequences, we hypothesize that CRISPR mediated clonal interference between resistant strains promotes and maintains diversity in this natural population.

Hyperthermophilic phosphotriesterases/lactonases for the environment and human health

In the last decades the idea to use enzymes for environmental bioremediation has been more and more proposed and, in the light of this, new solutions have been suggested and detailed studies on some classes of enzymes have been performed. In particular, our attention in the last few years has been focused on the enzymes belonging to the amidohydrolase superfamily. Several members of this superfamily are endowed with promiscuous activities. The term 'catalytic promiscuity' describes the capability of an enzyme to catalyse different chemical reactions, called secondary activities, at the active site responsible for the main activity. Recently, a new family of microbial lactonases with promiscuous phosphotriesterase activity, dubbed PTE-Like Lactonase (PLL), has been ascribed to the amidohydrolase superfamily. Among members of this family are enzymes found in the archaea Sulfolobus solfataricus and Sulfolobus acidocaldarius, which show high thermophilicity and thermal resistance. Enzymes showing phosphotriesterase activity are attractive from a biotechnological point of view because they are capable of hydrolysing the organophosphate phosphotriesters (OPs), a class of synthetic compounds employed worldwide both as insecticides and chemical warfare agents. Furthermore, from a basic point of view, studies of catalytic promiscuity offer clues to understand natural evolution of enzymes and to translate this into in vitro adaptation of enzymes to specific human needs. Thermostable enzymes able to hydrolyse OPs are considered good candidates for the set-up of efficient detoxification tools.

Unmarked gene deletion and host-vector system for the hyperthermophilic crenarchaeon Sulfolobus islandicus

Sulfolobus islandicus is being used as a model for studying archaeal biology, geo-biology and evolution. However, no genetic system is available for this organism. To produce an S. islandicus mutant suitable for genetic analyses, we screened for colonies with a spontaneous pyrEF mutation. One mutant was obtained containing only 233 bp of the original pyrE sequence in the mutant allele and it was used as a host to delete the beta-glycosidase (lacS) gene. Two unmarked gene deletion methods were employed, namely plasmid integration and segregation, and marker replacement and looping out, and unmarked lacS mutants were obtained by each method. A new alternative recombination mechanism, i.e., marker circularization and integration, was shown to operate in the latter method, which did not yield the designed deletion mutation. Subsequently, Sulfolobus-E. coli plasmid shuttle vectors were constructed, which genetically complemented DeltapyrEFDeltalacS mutation after transformation. Thus, a complete set of genetic tools was established for S. islandicus with pyrEF and lacS as genetic markers.

Evolution of complex RNA polymerases: the complete archaeal RNA polymerase structure

The archaeal RNA polymerase (RNAP) shares structural similarities with eukaryotic RNAP II but requires a reduced subset of general transcription factors for promoter-dependent initiation. To deepen our knowledge of cellular transcription, we have determined the structure of the 13-subunit DNA-directed RNAP from Sulfolobus shibatae at 3.35 Å resolution. The structure contains the full complement of subunits, including RpoG/Rpb8 and the equivalent of the clamp-head and jaw domains of the eukaryotic Rpb1. Furthermore, we have identified subunit Rpo13, an RNAP component in the order Sulfolobales, which contains a helix-turn-helix motif that interacts with the RpoH/Rpb5 and RpoA'/Rpb1 subunits. Its location and topology suggest a role in the formation of the transcription bubble.

[Current status on proteomics of extremophilic microorganisms--a review]

We summarized the key handicap and troubleshooting when proteomic techniques were used to investigate extremophilic microorganisms, and the actual state of their proteomes research in recent years. Up to now, proteomics techniques keep developing and improving rapidly, but they has not been widely used to explore proteome of extremophilic microorganisms including halophiles, thermophiles, psychophiles, acidophiles and alkaliphiles due to specific problems including incomplete dissociation of protein-protein complexes of extremophiles, and a lot of proteins synthesized by extremophiles are resistant to the conditions which dissociated and denatured proteins synthesized by mesophilic organisms. However, the foreground of potential application of the techniques draws people on attempting zealously multifarious methods. At the present time, several technical problems for separating halophilic proteins, integral membrane proteins and predicting the function of new proteins have been solved availably. Proteomics data have validated some conclusions of genome predictions, and revealed many novel proteins and a few properties of extremophiles can not be resolved fully by genome data. The investigation of extremophiles proteomes indicated that a comprehensive view of protein expression profiles should rely on more than one proteomic method. In addition, the mutual verification of conclusions on the basis of genome and proteome and combination of these two techniques must accelerate the study of extremophilic microorganisms, and redound to uncover deeply and wholly the unique mechanisms of microorganisms adaptation to extreme environments. Moreover, it would clarify the mechanisms of their survival, and point out new direction of survey for improving damage result from stresses, finally contribute to human survival and healthy.

Structure of a putative molybdenum-cofactor biosynthesis protein C (MoaC) from Sulfolobus tokodaii (ST0472)

The crystal structure of a putative molybdenum-cofactor (Moco) biosynthesis protein C (MoaC) from Sulfolobus tokodaii (ST0472) was determined at 2.2 A resolution. The crystal belongs to the monoclinic space group C2, with unit-cell parameters a = 123.31, b = 78.58, c = 112.67 A, beta = 118.1 degrees . The structure was solved by molecular replacement using the structure of Escherichia coli MoaC as the probe model. The asymmetric unit is composed of a hexamer arranged as a trimer of dimers with noncrystallographic 32 symmetry. The structure of ST0472 is very similar to that of E. coli MoaC; however, in the ST0472 protein an additional loop formed by the insertion of seven residues participates in intermonomer interactions and the new structure also reveals the formation of an interdimer beta-sheet. These features may contribute to the stability of the oligomeric state.

Genetic tools for Sulfolobus spp.: vectors and first applications

Sulfolobus species belong to the best-studied archaeal organisms but have lacked powerful genetic methods. Recently, there has been considerable progress in the field of Sulfolobus genetics. Urgently needed basic genetic tools, such as targeted gene knockout techniques and shuttle vectors are being developed at an increasing pace. For S. solfataricus knockout systems as well as different shuttle vectors are available. For the genetically more stable S. acidocaldarius shuttle vectors have been recently developed. In this review we summarize the currently available genetic tools and methods for the genus Sulfolobus. Different transformation protocols are discussed, as well as all so far developed knockout systems and Sulfolobus-Escherichia coli shuttle vectors are summarized. Special emphasis is put on the important vector components, i.e., selectable markers and Sulfolobus replicons. Additionally, the information gathered on different Sulfolobus strains with respect to their use as recipient strains is reviewed. The advantages and disadvantages of the different systems are discussed and aims for further improvement of genetic systems are identified.

[Cloning, expression and radiation inducibility of RadA from the hyperthermophilic archaeon Sulfolobus tokodaii]

RecA/Rad51/RadA recombinases are important recombination proteins with conserved functions. Studies on the enzymes have shown that members of RecA/Rad51/RadA family from bacteria, eukaryota, methanogens and halophilic archaea have UV inducibility. However, the UV inducibility of RadA homologues from hyperthermophilic archaea is controversial. We analyzed the UV inducibility of Sulfolobus tokodaii RadA by RT-PCR and immune assays. Comparing with the mock, the transcription and expression of the radA increased 2 and 1.5 folds respectively after UV irradiation at 100 J/m2, or 3 and 1 fold at 200 J/m2. These results demonstrated that S. tokodaii RadA could be induced after UV treatment. In addition, proteome induction analysis proved that there existed a DNA damage induction response in S. tokodaii, which further supported RadA inductility in this hyperthermophilic archaeon.

Virus movement maintains local virus population diversity

Viruses are the largest reservoir of genetic material on the planet, yet little is known about the population dynamics of any virus within its natural environment. Over a 2-year period, we monitored the diversity of two archaeal viruses found in hot springs within Yellowstone National Park (YNP). Both temporal phylogeny and neutral biodiversity models reveal that virus diversity in these local environments is not being maintained by mutation but rather by high rates of immigration from a globally distributed metacommunity. These results indicate that geographically isolated hot springs are readily able to exchange viruses. The importance of virus movement is supported by the detection of virus particles in air samples collected over YNP hot springs and by their detection in metacommunity sequencing projects conducted in the Sargasso Sea. Rapid rates of virus movement are not expected to be unique to these archaeal viruses but rather a common feature among virus metacommunities. The finding that virus immigration rather than mutation can dominate community structure has significant implications for understanding virus circulation and the role that viruses play in ecology and evolution by providing a reservoir of mobile genetic material.

Carboxyl pK(a) values, ion pairs, hydrogen bonding, and the pH-dependence of folding the hyperthermophile proteins Sac7d and Sso7d

Sac7d and Sso7d are homologous, hyperthermophile proteins with a high density of charged surface residues and potential ion pairs. To determine the relative importance of specific amino acid side-chains in defining the stability and function of these Archaeal chromatin proteins, pK(a) values were measured for the acidic residues in both proteins using (13)C NMR chemical shifts. The stability of Sso7d enabled titrations to pH 1 under low-salt conditions. Two aspartate residues in Sso7d (D16 and D35) and a single glutamate residue (G54) showed significantly perturbed pK(a) values in low salt, indicating that the observed pH-dependence of stability was primarily due to these three residues. The pH-dependence of backbone amide NMR resonances demonstrated that perturbation of all three pK(a) values was primarily the result of side-chain to backbone amide hydrogen bonds. Few of the significantly perturbed acidic pK(a) values in Sac7d and Sso7d could be attributed to primarily ion pair or electrostatic interactions. A smaller perturbation of E48 (E47 in Sac7d) was ascribed to an ion pair interaction that may be important in defining the DNA binding surface. The small number (three) of significantly altered pK(a) values was in good agreement with a linkage analysis of the temperature, pH, and salt-dependence of folding. The linkage of the ionization of two or more side-chains to protein folding led to apparent cooperativity in the pH-dependence of folding, although each group titrated independently with a Hill coefficient near unity. These results demonstrate that the acid pH-dependence of protein stability in these hyperthermophile proteins is due to independent titration of acidic residues with pK(a) values perturbed primarily by hydrogen bonding of the side-chain to the backbone. This work demonstrates the need for caution in using structural data alone to argue the importance of ion pairs in stabilizing hyperthermophile proteins.

The cell cycle of Sulfolobus

Much of the current information about the archaeal cell cycle has been generated through studies of the genus Sulfolobus. The overall organization of the cell cycle in these species is well understood, and information about the regulatory principles that govern cell cycle progression is rapidly accumulating. Exciting progress regarding the control and molecular details of the chromosome replication process is evident, and the first insights into the elusive crenarchaeal mitosis and cytokinesis machineries are within reach.

Variation in gene content among geographically diverse Sulfolobus isolates

The ability of competitive (i.e., comparative) genomic hybridization (CGH) to assess similarity across entire microbial genomes suggests that it should reveal diversification within and between natural populations of free-living prokaryotes. We used CGH to measure relatedness of genomes drawn from Sulfolobus populations that had been shown in a previous study to be diversified along geographical lines. Eight isolates representing a wide range of spatial separation were compared with respect to gene-specific tags based on a closely related reference strain (Sulfolobus solfataricus P2). For the purpose of assessing genetic divergence, 232 loci identified as polymorphic were assigned one of two alleles based on the corresponding fluorescence intensities from the arrays. Clustering of these binary genotypes was stable with respect to changes in the threshold and similarity criteria, and most of the groupings were consistent with an isolation-by-distance model of diversification. These results indicate that increasing spatial separation of geothermal sites correlates not only with minor sequence polymorphisms in conserved genes of Sulfolobus (demonstrated in the previous study), but also with the regions of difference (RDs) that occur between genomes of conspecifics. In view of the abundance of RDs in prokaryotic genomes and the relevance that some RDs may have for ecological adaptation, the results further suggest that CGH on microarrays may have advantages for investigating patterns of diversification in other free-living archaea and bacteria.

An order–disorder twin crystal ofL-2-haloacid dehalogenase fromSulfolobus tokodaii

The L-2-haloacid dehalogenase enzymes catalyse the hydrolytic cleavage of a halogen from the C2 position of short-chain haloacids. The recombinant dehalogenase from the thermophilic archaeon Sulfolobus tokodaii has been cloned, overexpressed and purified to homogeneity. The 24 kDa enzyme was crystallized using the microbatch method in the monoclinic space group C2, with unit-cell parameters a = 127.6, b = 58.1, c = 51.2 A, beta = 97.2 degrees . Data were collected to 1.9 A resolution using synchrotron radiation and the structure was solved by molecular replacement. Analysis of the data and the preliminary refined model showed that the crystal was an order-disorder twin by reticular merohedry with a twin index of 10. It was possible to de-twin the experimental data utilizing the symmetry of the molecular layers from which the crystal is built.

Transcriptional analysis of the genetic element pSSVx: differential and temporal regulation of gene expression reveals correlation between transcription and replication

pSSVx from Sulfolobus islandicus strain REY15/4 is a hybrid between a plasmid and a fusellovirus. A systematic study performed by a combination of Northern blot analysis, primer extension, and reverse transcriptase PCR revealed the presence of nine major transcripts whose expression was differentially and temporally regulated over the growth cycle of S. islandicus. The map positions of the RNAs as well as the clockwise and the anticlockwise directions of their transcription were determined. Some genes were clustered and appeared to be transcribed as polycistronic messengers, among which one long transcriptional unit comprised the genes for the plasmid copy number control protein ORF60 (CopG), ORF91, and the replication protein ORF892 (RepA). We propose that a termination readthrough mechanism might be responsible for the formation of more than one RNA species from a single 5' end and therefore that the nine different RNAs corresponded to only seven different transcriptional starts. Three transcripts, ORF76 and two antisense RNAs, countertranscribed RNA1 (ctRNA1) and ctRNA2, were found to be specifically expressed during (and hence correlated to) the phase in which the pSSVx copy number is kept under stringent control, as they were completely switched off upon the onset of the induction of replication.

Cloning and characterization of glycogen-debranching enzyme from hyperthermophilic archaeon Sulfolobus shibatae

A gene encoding a putative glycogen-debranching enzyme in Sulfolobus shibatae (abbreviated as SSGDE) was cloned and expressed in Escherichia coli. The recombinant enzyme was purified to homogeneity by heat treatment and Ni-NTA affinity chromatography. The recombinant SSGDE was extremely thermostable, with an optimal temperature at 85 degrees C. The enzyme had an optimum pH of 5.5 and was highly stable from pH 4.5 to 6.5. The substrate specificity of SSGDE suggested that it possesses characteristics of both amylo-1,6-glucosidase and alpha-1,4-glucanotransferase. SSGDE clearly hydrolyzed pullulan to maltotrlose, and 6-O-alpha-maltosyl-beta-cyclodextrin (G2-beta-CD) to maltose and beta-cyclodextrin. At the same time, SSGDE transferred maltooligosyl residues to the maltooligosaccharides employed, and maltosyl residues to G2-beta-CD. The enzyme preferentially hydrolyzed amylopectin, followed in a decreasing order by glycogen, pullulan, and amylose. Therefore, the present results suggest that the glycogen-debranching enzyme from S. shibatae may have industrial application for the efficient debranching and modification of starch to dextrins at a high temperature.

Elucidating the transcription cycle of the UV-inducible hyperthermophilic archaeal virus SSV1 by DNA microarrays

The spindle-shaped Sulfolobus virus SSV1 was the first of a series of unusual and uniquely shaped viruses isolated from hyperthermophilic Archaea. Using whole-genome microarrays we show here that the circular 15.5 kb DNA genome of SSV1 exhibits a chronological regulation of its transcription upon UV irradiation, reminiscent to the life cycles of bacteriophages and eukaryotic viruses. The transcriptional cycle starts with a small UV-specific transcript and continues with early transcripts on both its flanks. The late transcripts appear after the onset of viral replication and are extended to their full lengths towards the end of the approximately 8.5 h cycle. While we detected only small differences in genome-wide analysis of the host Sulfolobus solfataricus comparing infected versus uninfected strains, we found a marked difference with respect to the strength and speed of the general UV response of the host. Models for the regulation of the virus cycle, and putative functions of genes in SSV1 are presented.

The N-terminal region is important for the nuclease activity and thermostability of the flap endonuclease-1 from Sulfolobus tokodaii

This paper reports the biochemical properties of two types of recombinant flap endonuclease-1 (FEN-1) proteins obtained from the thermophilic crenarchaeon, Sulfolobus tokodaii strain 7. One of the two FEN-1 proteins is a product of the gene with AUG as the translational start codon (StoS-FEN-1), which is originally assigned in the database. The other is a product of the gene with a new AUG start codon (StoL-FEN-1), which is inserted at 153 bases upstream of the original AUG codon. Although StoL-FEN-1 showed activity and thermostability, StoS-FEN-1 showed neither activity nor thermostability. The N-terminal region in StoL-FEN-1 was also conserved in all of the FEN-1 homologs deduced from genes from newly isolated Sulfolobus spp. These results strongly suggest that the actual start codon of the fen-1 gene from S. tokodaii is not the originally assigned AUG, but rather is located at about 100 bases upstream of this codon.

Role of the Sulfolobus shibatae viral T6 initiator in conferring promoter strength and in influencing transcription start site selection

Archaeal promoters contain a TATA-box, an adjacent upstream TFB-recognition element (BRE), and a downstream initiator (INR) region from which transcription originates. While the contribution of A-box and BRE to promoter strength is well established, the role of DNA sequences within the INR region and its vicinity on transcription efficiency and start site selection remains unclear. Here, I demonstrate using the strong Sulfolobus shibatae viral T6 promoter that either substitution of its natural sequence from -17 and beyond with plasmid DNA or introduction of point transversion mutations at +3, -2, -4, and -5 positions reduce promoter strength dramatically, whereas +1, -1, and -2 mutations influence the transcription start site. These data therefore reveal that the INR region plays a role as important as the BRE and the A-box in T6 gene transcription.

Enzymatic characterization and substrate specificity of thermostable beta-glycosidase from hyperthermophilic archaea, Sulfolobus shibatae, expressed in E. coli

Enzymatic properties and substrate specificity of recombinant beta-glycosidases from a hyperthermophilic archaeon, Sulfolobus shibatae (rSSG), were analyzed. rSSG showed its optimum temperature and pH at 95 degrees C and pH 5.0, respectively. Thermal inactivation of rSSG showed that its half-life of enzymatic activity at 75 degrees C was 15 h whereas it drastically decreased to 3.9 min at 95 degrees C. The addition of 10 mM of MnCl2 enhanced the hydrolysis activity of rSSG up to 23% whereas most metal ions did not show any considerable effect. Dithiothreitol (DTT) and 2-mercaptoethanol exhibited significant influence on the increase of the hydrolysis activity of rSSG. rSSG apparently preferred laminaribiose (beta1-->3Glc), followed by sophorose (beta1-->2Glc), gentiobiose (beta1-->6Glc), and cellobiose (beta1--4Glc). Various intermolecular transfer products were formed by rSSG in the lactose reaction, indicating that rSSG prefers lactose as a good acceptor as well as a donor. The strong intermolecular transglycosylation activity of rSSG can be applied in making functional oligosaccharides.