Genome sequence of the hyperthermophilic crenarchaeon Pyrobaculum aerophilum

Antibacterial Mechanism of d-Cysteine/Polyethylene Glycol-Functionalized Gold Nanoparticles and Their Potential for the Treatment of Bacterial Infections

Bacterial infection has always posed a severe threat to public health. Gold nanoparticles (Au NPs) exhibit exceptional biocompatibility and hold immense potential in biomedical applications. However, their antibacterial effectiveness is currently unsatisfactory. Herein, a chiral antibacterial agent with high stability was prepared by the modification of Au NPs with d-cysteine with the assistance of polyethylene glycol (PEG). The as-synthesized d-cysteine/PEG-Au NPs (D/P-Au NPs) exhibited a stronger (99.5-99.9%) and more stable (at least 14 days) antibacterial performance against Gram-negative (Escherichia coli and Listeria monocytogenes) and Gram-positive (Salmonella enteritidis and Staphylococcus aureus) bacteria, compared with other groups. The analysis of the antibacterial mechanism revealed that the D/P-Au NPs mainly affected the assembly of ribosomes, the biosynthesis of amino acids and proteins, as well as the DNA replication and mismatch repair, ultimately leading to bacterial death, which is significantly different from the mechanism of reactive oxygen species-activated metallic antibacterial NPs. In particular, the D/P-Au NPs were shown to effectively accelerate the healing of S. aureus-infected wounds in mice to a rate comparable to or slightly higher than that of vancomycin. This work provides a novel approach to effectively design chiral antibacterial agents for bacterial infection treatment.

MiDSystem: A comprehensive online system for de novo assembly and analysis of microbial genomes

The substantial reduction in experimental cost of next-generation sequencing techniques makes it feasible to assemble a bacterial genome of unknown species de novo and acquire substantial genetic information from environmental samples. Many bioinformatics tools and algorithms have also been developed for prokaryotes, but complex parameter settings and command line-based user interfaces cause a significant entry barrier for novices. Efficient construction of pipelines that integrate all the available genomic data poses a major challenge to the understanding of unknown pathogens. MiDSystem is a comprehensive online system for analyzing genomic data from microbiomes. With a user-friendly interface, MiDSystem supports both de novo assembly and metagenomic analysis pipelines. It is designed to automatically analyze whole genome shotgun sequencing data of bacteria submitted by users. Multiple analytical steps can be performed directly on the system, and the results generated from the embedded tools are visualized in an online summary report to make it more interpretable. Constructing a genome de novo has gradually become the foundation of bacterial studies. Taking both single species and metagenomic samples into consideration, MiDSystem can greatly reduce the time and effort for analysis of bacterial genomic data. Use of MiDSystem will enable more focus to be placed on understanding the etiology of bacterial infections and microorganism ecologies.

Comprehensive predictions of secondary structures for comparative analysis in different species

Protein structures are directly linked to biological functions. However, there is a gap of knowledge between the decoded genome and the structure. To bridge the gap, we focused on the secondary structure (SS). From a comprehensive analysis of predicted SS of proteins in different types of organisms, we have arrived at the following findings: The proportions of SS in genomes were different among phylogenic domains. The distributions of strand lengths indicated structural limitations in all of the species. Different from bacteria and archaea, eukaryotes have an abundance of α-helical and random coil proteins. Interestingly, there was a relationship between SS and post-translational modifications. By calculating hydrophobicity moments of helices and strands, highly amphipathic fragments of SS were found, which might be related to the biological functions. In conclusion, comprehensive predictions of SS will provide valuable perspectives to understand the entire protein structures in genomes and will help one to discover or design functional proteins.

Functional characteristic of microbial communities in large-scale biotreatment systems of food waste

In order to evaluate microbial community structure dominated metabolic function profiles in large-scale food waste (FW) biotreatment systems, bacterial, archaeal and fungal community associated with metabolic function in high-temperature aerobic fermentation (AF) and anaerobic co-digestion (AcoD) processes were comprehensively investigated in this study. The qPCR results showed the higher gene copies of bacteria and fungi in initial and AF-treated FW compared with AcoD-treated FW, as well as bacteria and archaea in AcoD-treated FW were highly abundant among detected samples. Furthermore, the total abundances of archaea ((1.18-4.88) × 10⁶ copies/ng DNA) in AcoD system were 2-3 orders of magnitude higher than that in other samples (P < 0.01), indicating active archaeal activity in AcoD system. Correlation analysis of microbial community and metabolic function indicated that the higher abundances of Kazachstania, Pyrobaculum, Sulfophobococcus, Lactobacillus and Candida in initial FW had close linkages with lipid metabolism (P < 0.05). Abundant Aspergillus, Staphylococcus, Pelomonas, Corynebacterium, Faecalibacterium, Methanobacterium and Xeromyces in AF system were positively and significantly correlated with high metabolic activities of energy metabolism, carbohydrate metabolism, amino acid metabolism, fatty acid metabolism, glycosaminoglycan degradation, sulfur metabolism and nitrogen metabolism. As for AcoD system, dominant genera Methanosaeta, Methanoculleus, Methanobacterium, Fastidiosipila, Rikenellaceae RC9, Bifidobacterium and Xeromyces had close relationships with metabolism of cofactors and vitamins, energy metabolism, methane metabolism, carbohydrate metabolism and glycosaminoglycan degradation (P < 0.05). These results are expected to improve the metabolic efficiency by functional microorganism in different large-scale FW treatment systems.

Improving Arsenic Tolerance of Pyrococcus furiosus by Heterologous Expression of a Respiratory Arsenate Reductase

Arsenate is a notorious toxicant that is known to disrupt multiple biochemical pathways. Many microorganisms have developed mechanisms to detoxify arsenate using the ArsC-type arsenate reductase, and some even use arsenate as a terminal electron acceptor for respiration involving arsenate respiratory reductase (Arr). ArsC-type reductases have been studied extensively, but the phylogenetically unrelated Arr system is less investigated and has not been characterized from Archaea Here, we heterologously expressed the genes encoding Arr from the crenarchaeon Pyrobaculum aerophilum in the euryarchaeon Pyrococcus furiosus, both of which grow optimally near 100°C. Recombinant P. furiosus was grown on molybdenum (Mo)- or tungsten (W)-containing medium, and two types of recombinant Arr enzymes were purified, one containing Mo (Arr-Mo) and one containing W (Arr-W). Purified Arr-Mo had a 140-fold higher specific activity in arsenate [As(V)] reduction than Arr-W, and Arr-Mo also reduced arsenite [As(III)]. The P. furiosus strain expressing Arr-Mo (the Arr strain) was able to use arsenate as a terminal electron acceptor during growth on peptides. In addition, the Arr strain had increased tolerance compared to that of the parent strain to arsenate and also, surprisingly, to arsenite. Compared to the parent, the Arr strain accumulated intracellularly almost an order of magnitude more arsenic when cells were grown in the presence of arsenite. X-ray absorption spectroscopy (XAS) results suggest that the Arr strain of P. furiosus improves its tolerance to arsenite by increasing production of less-toxic arsenate and nontoxic methylated arsenicals compared to that by the parent.IMPORTANCE Arsenate respiratory reductases (Arr) are much less characterized than the detoxifying arsenate reductase system. The heterologous expression and characterization of an Arr from Pyrobaculum aerophilum in Pyrococcus furiosus provides new insights into the function of this enzyme. From in vivo studies, production of Arr not only enabled P. furiosus to use arsenate [As(V)] as a terminal electron acceptor, it also provided the organism with a higher resistance to arsenate and also, surprisingly, to arsenite [As(III)]. In contrast to the tungsten-containing oxidoreductase enzymes natively produced by P. furiosus, recombinant P. aerophilum Arr was much more active with molybdenum than with tungsten. It is also, to our knowledge, the only characterized Arr to be active with both molybdenum and tungsten in the active site.

Enzymatic Antioxidant Signatures in Hyperthermophilic Archaea

To fight reactive oxygen species (ROS) produced by both the metabolism and strongly oxidative habitats, hyperthermophilic archaea are equipped with an array of antioxidant enzymes whose role is to protect the biological macromolecules from oxidative damage. The most common ROS, such as superoxide radical (O₂^-.) and hydrogen peroxide (H₂O₂), are scavenged by superoxide dismutase, peroxiredoxins, and catalase. These enzymes, together with thioredoxin, protein disulfide oxidoreductase, and thioredoxin reductase, which are involved in redox homeostasis, represent the core of the antioxidant system. In this review, we offer a panorama of progression of knowledge on the antioxidative system in aerobic or microaerobic (hyper)thermophilic archaea and possible industrial applications of these enzymes.

Methylation guide RNA evolution in archaea: structure, function and genomic organization of 110 C/D box sRNA families across six Pyrobaculum species

Archaeal homologs of eukaryotic C/D box small nucleolar RNAs (C/D box sRNAs) guide precise 2′-O-methyl modification of ribosomal and transfer RNAs. Although C/D box sRNA genes constitute one of the largest RNA gene families in archaeal thermophiles, most genomes have incomplete sRNA gene annotation because reliable, fully automated detection methods are not available. We expanded and curated a comprehensive gene set across six species of the crenarchaeal genus Pyrobaculum, particularly rich in C/D box sRNA genes. Using high-throughput small RNA sequencing, specialized computational searches and comparative genomics, we analyzed 526 Pyrobaculum C/D box sRNAs, organizing them into 110 families based on synteny and conservation of guide sequences which determine methylation targets. We examined gene duplications and rearrangements, including one family that has expanded in a pattern similar to retrotransposed repetitive elements in eukaryotes. New training data and inclusion of kink-turn secondary structural features enabled creation of an improved search model. Our analyses provide the most comprehensive, dynamic view of C/D box sRNA evolutionary history within a genus, in terms of modification function, feature plasticity, and gene mobility.

The catalytic mechanism of electron-bifurcating electron transfer flavoproteins (ETFs) involves an intermediary complex with NAD<sup/>

Electron bifurcation plays a key role in anaerobic energy metabolism, but it is a relatively new discovery, and only limited mechanistic information is available on the diverse enzymes that employ it. Herein, we focused on the bifurcating electron transfer flavoprotein (ETF) from the hyperthermophilic archaeon Pyrobaculum aerophilum The EtfABCX enzyme complex couples NADH oxidation to the endergonic reduction of ferredoxin and exergonic reduction of menaquinone. We developed a model for the enzyme structure by using nondenaturing MS, cross-linking, and homology modeling in which EtfA, -B, and -C each contained FAD, whereas EtfX contained two [4Fe-4S] clusters. On the basis of analyses using transient absorption, EPR, and optical titrations with NADH or inorganic reductants with and without NAD+, we propose a catalytic cycle involving formation of an intermediary NAD+-bound complex. A charge transfer signal revealed an intriguing interplay of flavin semiquinones and a protein conformational change that gated electron transfer between the low- and high-potential pathways. We found that despite a common bifurcating flavin site, the proposed EtfABCX catalytic cycle is distinct from that of the genetically unrelated bifurcating NADH-dependent ferredoxin NADP+ oxidoreductase (NfnI). The two enzymes particularly differed in the role of NAD+, the resting and bifurcating-ready states of the enzymes, how electron flow is gated, and the two two-electron cycles constituting the overall four-electron reaction. We conclude that P. aerophilum EtfABCX provides a model catalytic mechanism that builds on and extends previous studies of related bifurcating ETFs and can be applied to the large bifurcating ETF family.

Structure and analysis of nucleoside diphosphate kinase from Borrelia burgdorferi prepared in a transition-state complex with ADP and vanadate moieties

Nucleoside diphosphate kinases (NDKs) are implicated in a wide variety of cellular functions owing to their enzymatic conversion of NDP to NTP. NDK from Borrelia burgdorferi (BbNDK) was selected for functional and structural analysis to determine whether its activity is required for infection and to assess its potential for therapeutic inhibition. The Seattle Structural Genomics Center for Infectious Diseases (SSGCID) expressed recombinant BbNDK protein. The protein was crystallized and structures were solved of both the apoenzyme and a liganded form with ADP and vanadate ligands. This provided two structures and allowed the elucidation of changes between the apo and ligand-bound enzymes. Infectivity studies with ndk transposon mutants demonstrated that NDK function was important for establishing a robust infection in mice, and provided a rationale for therapeutic targeting of BbNDK. The protein structure was compared with other NDK structures found in the Protein Data Bank and was found to have similar primary, secondary, tertiary and quaternary structures, with conserved residues acting as the catalytic pocket, primarily using His132 as the phosphohistidine-transfer residue. Vanadate and ADP complexes model the transition state of this phosphoryl-transfer reaction, demonstrating that the pocket closes when bound to ADP, while allowing the addition or removal of a γ-phosphate. This analysis provides a framework for the design of potential therapeutics targeting BbNDK inhibition.

Genomic Comparison of Two Family-Level Groups of the Uncultivated NAG1 Archaeal Lineage from Chemically and Geographically Disparate Hot Springs

Recent progress based on single-cell genomics and metagenomic investigations of archaea in a variety of extreme environments has led to significant advances in our understanding of the diversity, evolution, and metabolic potential of archaea, yet the vast majority of archaeal diversity remains undersampled. In this work, we coordinated single-cell genomics with metagenomics in order to construct a near-complete genome from a deeply branching uncultivated archaeal lineage sampled from Great Boiling Spring (GBS) in the U.S. Great Basin, Nevada. This taxon is distantly related (distinct families) to an archaeal genome, designated "Novel Archaeal Group 1" (NAG1), which was extracted from a metagenome recovered from an acidic iron spring in Yellowstone National Park (YNP). We compared the metabolic predictions of the NAG1 lineage to better understand how these archaea could inhabit such chemically distinct environments. Similar to the NAG1 population previously studied in YNP, the NAG1 population from GBS is predicted to utilize proteins as a primary carbon source, ferment simple carbon sources, and use oxygen as a terminal electron acceptor under oxic conditions. However, GBS NAG1 populations contained distinct genes involved in central carbon metabolism and electron transfer, including nitrite reductase, which could confer the ability to reduce nitrite under anaerobic conditions. Despite inhabiting chemically distinct environments with large variations in pH, GBS NAG1 populations shared many core genomic and metabolic features with the archaeon identified from YNP, yet were able to carve out a distinct niche at GBS.

The Structure, Function and Roles of the Archaeal ESCRT Apparatus

Although morphologically resembling bacteria, archaea constitute a distinct domain of life with a closer affiliation to eukaryotes than to bacteria. This similarity is seen in the machineries for a number of essential cellular processes, including DNA replication and gene transcription. Perhaps surprisingly, given their prokaryotic morphology, some archaea also possess a core cell division apparatus that is related to that involved in the final stages of membrane abscission in vertebrate cells, the ESCRT machinery.

Pyrobaculum yellowstonensis Strain WP30 Respires on Elemental Sulfur and/or Arsenate in Circumneutral Sulfidic Geothermal Sediments of Yellowstone National Park

Thermoproteales (phylum Crenarchaeota) populations are abundant in high-temperature (>70°C) environments of Yellowstone National Park (YNP) and are important in mediating the biogeochemical cycles of sulfur, arsenic, and carbon. The objectives of this study were to determine the specific physiological attributes of the isolate Pyrobaculum yellowstonensis strain WP30, which was obtained from an elemental sulfur sediment (Joseph's Coat Hot Spring [JCHS], 80°C, pH 6.1, 135 μM As) and relate this organism to geochemical processes occurring in situ. Strain WP30 is a chemoorganoheterotroph and requires elemental sulfur and/or arsenate as an electron acceptor. Growth in the presence of elemental sulfur and arsenate resulted in the formation of thioarsenates and polysulfides. The complete genome of this organism was sequenced (1.99 Mb, 58% G+C content), revealing numerous metabolic pathways for the degradation of carbohydrates, amino acids, and lipids. Multiple dimethyl sulfoxide-molybdopterin (DMSO-MPT) oxidoreductase genes, which are implicated in the reduction of sulfur and arsenic, were identified. Pathways for the de novo synthesis of nearly all required cofactors and metabolites were identified. The comparative genomics of P. yellowstonensis and the assembled metagenome sequence from JCHS showed that this organism is highly related (∼95% average nucleotide sequence identity) to in situ populations. The physiological attributes and metabolic capabilities of P. yellowstonensis provide an important foundation for developing an understanding of the distribution and function of these populations in YNP.

Understanding DNA Repair in Hyperthermophilic Archaea: Persistent Gaps and Other Reasons to Focus on the Fork

Although hyperthermophilic archaea arguably have a great need for efficient DNA repair, they lack members of several DNA repair protein families broadly conserved among bacteria and eukaryotes. Conversely, the putative DNA repair genes that do occur in these archaea often do not generate the expected phenotype when deleted. The prospect that hyperthermophilic archaea have some unique strategies for coping with DNA damage and replication errors has intellectual and technological appeal, but resolving this question will require alternative coping mechanisms to be proposed and tested experimentally. This review evaluates a combination of four enigmatic properties that distinguishes the hyperthermophilic archaea from all other organisms: DNA polymerase stalling at dU, apparent lack of conventional NER, lack of MutSL homologs, and apparent essentiality of homologous recombination proteins. Hypothetical damage-coping strategies that could explain this set of properties may provide new starting points for efforts to define how archaea differ from conventional models of DNA repair and replication fidelity.

Properties of bacterial laccases and their application in bioremediation of industrial wastes

The bioremediation process of industrial waste can be made more efficient using ligninolytic laccase enzymes, which are obtained from fungi, bacteria, higher plants, insects, and also in lichen. Laccase are catalyzed in the mono-electronic oxidation of a substrate from the expenditure of molecular oxygen. This enzyme belongs to the multicopper oxidases and participates in the cross linking of monomers, involved in the degradation of wide range industrial pollutants. In recent years, these enzymes have gained application in pulp and paper, textile and food industries. There are numerous reviews on laccases; however, a lot of information is still unknown due to their broad range of functions and applications. In this review, the bacterial laccases are focused for the bioremediation of various industrial pollutants. A brief description on structural molecular and physicochemical properties has been made. Moreover, the mechanism by which the reaction is catalyzed, the physical basis of thermostability and enantioselectivity, which requires more attention from researchers, and applications of laccase in various fields of biotechnology are pointed out.

Small nucleolar RNAs and RNA-guided post-transcriptional modification

snoRNAs (small nucleolar RNAs) constitute one of the largest and best-studied classes of non-coding RNAs that confer enzymatic specificity. With associated proteins, these snoRNAs form ribonucleoprotein complexes that can direct 2'-O-methylation or pseudouridylation of target non-coding RNAs. Aided by computational methods and high-throughput sequencing, new studies have expanded the diversity of known snoRNA functions. Complexes incorporating snoRNAs have dynamic specificity, and include diverse roles in RNA silencing, telomerase maintenance and regulation of alternative splicing. Evidence that dysregulation of snoRNAs can cause human disease, including cancer, indicates that the full scope of snoRNA roles remains an unfinished story. The diversity in structure, genomic origin and function between snoRNAs found in different complexes and among different phyla illustrates the surprising plasticity of snoRNAs in evolution. The ability of snoRNAs to direct highly specific interactions with other RNAs is a consistent thread in their newly discovered functions. Because they are ubiquitous throughout Eukarya and Archaea, it is likely they were a feature of the last common ancestor of these two domains, placing their origin over two billion years ago. In the present chapter, we focus on recent advances in our understanding of these ancient, but functionally dynamic RNA-processing machines.

Comparative genomic analysis reveals 2-oxoacid dehydrogenase complex lipoylation correlation with aerobiosis in archaea

Metagenomic analyses have advanced our understanding of ecological microbial diversity, but to what extent can metagenomic data be used to predict the metabolic capacity of difficult-to-study organisms and their abiotic environmental interactions? We tackle this question, using a comparative genomic approach, by considering the molecular basis of aerobiosis within archaea. Lipoylation, the covalent attachment of lipoic acid to 2-oxoacid dehydrogenase multienzyme complexes (OADHCs), is essential for metabolism in aerobic bacteria and eukarya. Lipoylation is catalysed either by lipoate protein ligase (LplA), which in archaea is typically encoded by two genes (LplA-N and LplA-C), or by a lipoyl(octanoyl) transferase (LipB or LipM) plus a lipoic acid synthetase (LipA). Does the genomic presence of lipoylation and OADHC genes across archaea from diverse habitats correlate with aerobiosis? First, analyses of 11,826 biotin protein ligase (BPL)-LplA-LipB transferase family members and 147 archaeal genomes identified 85 species with lipoylation capabilities and provided support for multiple ancestral acquisitions of lipoylation pathways during archaeal evolution. Second, with the exception of the Sulfolobales order, the majority of species possessing lipoylation systems exclusively retain LplA, or either LipB or LipM, consistent with archaeal genome streamlining. Third, obligate anaerobic archaea display widespread loss of lipoylation and OADHC genes. Conversely, a high level of correspondence is observed between aerobiosis and the presence of LplA/LipB/LipM, LipA and OADHC E2, consistent with the role of lipoylation in aerobic metabolism. This correspondence between OADHC lipoylation capacity and aerobiosis indicates that genomic pathway profiling in archaea is informative and that well characterized pathways may be predictive in relation to abiotic conditions in difficult-to-study extremophiles. Given the highly variable retention of gene repertoires across the archaea, the extension of comparative genomic pathway profiling to broader metabolic and homeostasis networks should be useful in revealing characteristics from metagenomic datasets related to adaptations to diverse environments.

Characterization of cleavage intermediate and star sites of RM.Tth111II

Tth111II is a thermostable Type IIGS restriction enzyme that recognizes DNA sites CAARCA (R = A or G) and cleaves downstream at N11/N9. Here, the tth111IIRM gene was cloned and expressed in E. coli, and Tth111II was purified. The purified enzyme contains internally-bound S-adenosylmethionine (SAM). When the internal SAM was removed, the endonuclease activity was stimulated by adding SAM or its analog sinefungin. The cleavage intermediate is mostly top-strand nicked DNA on a single-site plasmid. Addition of duplex oligos with a cognate site stimulates cleavage activity of the one-site substrate. Tth111II cleaves a two-site plasmid DNA with equal efficiency regardless of site orientation. We propose the top-strand nicking is carried out by a Tth111II monomer and bottom-strand cleavage is carried out by a transient dimer. Tth111II methylates cleavage product-like duplex oligos CAAACAN9, but the modification rate is estimated to be much slower than the top-strand nicking rate. We cloned and sequenced a number of Tth111II star sites which are 1-bp different from the cognate sites. A biochemical pathway is proposed for the restriction and methylation activities of Tth111II.

Evolution of the thermopsin peptidase family (A5)

Thermopsin is a peptidase from Sulfolobus acidocaldarius that is active at low pH and high temperature. From reversible inhibition with pepstatin, thermopsin is thought to be an aspartic peptidase. It is a member of the only family of peptidases to be restricted entirely to the archaea, namely peptidase family A5. Evolution within this family has been mapped, using a taxonomic tree based on the known classification of archaea. Homologues are found only in archaeans that are both hyperthermophiles and acidophiles, and this implies lateral transfer of genes between archaea, because species with homologues are not necessarily closely related. Despite the remarkable stability and activity in extreme conditions, no tertiary structure has been solved for any member of the family, and the catalytic mechanism is unknown. Putative catalytic residues have been predicted here by examination of aligned sequences.

The impact of genomics on research in diversity and evolution of archaea

Since the definition of archaea as a separate domain of life along with bacteria and eukaryotes, they have become one of the most interesting objects of modern microbiology, molecular biology, and biochemistry. Sequencing and analysis of archaeal genomes were especially important for studies on archaea because of a limited availability of genetic tools for the majority of these microorganisms and problems associated with their cultivation. Fifteen years since the publication of the first genome of an archaeon, more than one hundred complete genome sequences of representatives of different phylogenetic groups have been determined. Analysis of these genomes has expanded our knowledge of biology of archaea, their diversity and evolution, and allowed identification and characterization of new deep phylogenetic lineages of archaea. The development of genome technologies has allowed sequencing the genomes of uncultivated archaea directly from enrichment cultures, metagenomic samples, and even from single cells. Insights have been gained into the evolution of key biochemical processes in archaea, such as cell division and DNA replication, the role of horizontal gene transfer in the evolution of archaea, and new relationships between archaea and eukaryotes have been revealed.

Physicochemical characterization of a thermostable alcohol dehydrogenase from Pyrobaculum aerophilum

In this work we characterize an alcohol dehydrogenase (ADH) from the hyperthermophilic archaeon Pyrobaculum aerophilum (PyAeADHII). We have previously found that PyAeADHII has no activity when standard ADH substrates are used but is active when α-tetralone is used as substrate. Here, to gain insights into enzyme function, we screened several chemical libraries for enzymatic modulators using an assay employing α-tetralone. The results indicate that PyAeADHII activity in the presence of α-tetralone was inhibited by compounds such as flunarizine. We also examined metal coordination of the enzyme in solution by performing metal substitution of the enzyme-bound zinc (Zn²⁺) with cobalt. The solution-based absorption spectra for cobalt substituted PyAeADHII supports substitution at the structural Zn²⁺ site. To gain structural insight, we obtained the crystal structure of both wild-type and cobalt-substituted PyAeADHII at 1.75 Å and 2.20 Å resolution, respectively. The X-ray data confirmed one metal ion per monomer present only at the structural site with otherwise close conservation to other ADH enzymes. We next determined the co-crystal structure of the NADPH-bound form of the enzyme at 2.35 Å resolution to help define the active site region of the enzyme and this data shows close structural conservation with horse ADH, despite the lack of a catalytic Zn²⁺ ion in PyAeADHII. Modeling of α-tetralone into the NADPH bound structure suggests an arginine as a possible catalytic residue. The data presented here can yield a better understanding of alcohol dehydrogenases lacking the catalytic zinc as well as the structural features inherent to thermostable enzymes.

Phylogenetic and Functional Analysis of Metagenome Sequence from High-Temperature Archaeal Habitats Demonstrate Linkages between Metabolic Potential and Geochemistry

Geothermal habitats in Yellowstone National Park (YNP) provide an unparalleled opportunity to understand the environmental factors that control the distribution of archaea in thermal habitats. Here we describe, analyze, and synthesize metagenomic and geochemical data collected from seven high-temperature sites that contain microbial communities dominated by archaea relative to bacteria. The specific objectives of the study were to use metagenome sequencing to determine the structure and functional capacity of thermophilic archaeal-dominated microbial communities across a pH range from 2.5 to 6.4 and to discuss specific examples where the metabolic potential correlated with measured environmental parameters and geochemical processes occurring in situ. Random shotgun metagenome sequence (∼40–45 Mb Sanger sequencing per site) was obtained from environmental DNA extracted from high-temperature sediments and/or microbial mats and subjected to numerous phylogenetic and functional analyses. Analysis of individual sequences (e.g., MEGAN and G + C content) and assemblies from each habitat type revealed the presence of dominant archaeal populations in all environments, 10 of whose genomes were largely reconstructed from the sequence data. Analysis of protein family occurrence, particularly of those involved in energy conservation, electron transport, and autotrophic metabolism, revealed significant differences in metabolic strategies across sites consistent with differences in major geochemical attributes (e.g., sulfide, oxygen, pH). These observations provide an ecological basis for understanding the distribution of indigenous archaeal lineages across high-temperature systems of YNP.

Interactions of archaeal chromatin proteins Alba1 and Alba2 with nucleic acids

BACKGROUND

Architectural proteins have important roles in compacting and organising chromosomal DNA. There are two potential histone counterpart peptide sequences (Alba1 and Alba2) in the Aeropyrum pernix genome (APE1832.1 and APE1823).

METHODOLOGY/PRINCIPAL FINDINGS

THESE TWO PEPTIDES WERE EXPRESSED AND THEIR INTERACTIONS WITH VARIOUS DNAS WERE STUDIED USING A COMBINATION OF VARIOUS EXPERIMENTAL TECHNIQUES: surface plasmon resonance, UV spectrophotometry, circular dichroism-spectropolarimetry, gel-shift assays, and isothermal titration calorimetry.

CONCLUSIONS/SIGNIFICANCE

Our data indicate that there are significant differences in the properties of the Alba1 and Alba2 proteins. Both of these Alba proteins can thermally stabilise DNA polynucleotides, as seen from UV melting curves. Alba2 and equimolar mixtures of Alba1/Alba2 have greater effects on the thermal stability of poly(dA-dT).poly(dA-dT). Surface plasmon resonance sensorgrams for binding of Alba1, Alba2, and equimolar mixtures of Alba1/Alba2 to DNA oligonucleotides show different binding patterns. Circular dichroism indicates that Alba2 has a less-ordered secondary structure than Alba1. The secondary structures of the Alba proteins are not significantly influenced by DNA binding, even at high temperatures. Based on these data, we conclude that Alba1, Alba2, and equimolar mixtures of Alba1/Alba2 show different properties in their binding to various DNAs.

Reclassification of Thermoproteus neutrophilus Stetter and Zillig 1989 as Pyrobaculum neutrophilum comb. nov. based on phylogenetic analysis

The hyperthermophilic crenarchaeon Thermoproteus neutrophilus V24StaT was originally classified before sequence-based phylogenetic analysis became standard for bacterial taxonomy. Subsequent phylogenetic analyses by various groups have shown that strain V24StaT groups more closely with strains of the genus Pyrobaculum than with those in the genus Thermoproteus . Based on phylogenetic comparison of rRNA gene sequences and ribosomal proteins, we propose that strain V24StaT be reclassified as Pyrobaculum neutrophilum comb. nov., with the type strain V24StaT ( = DSM 2338T = JCM 9278T). An emended description of the genus Pyrobaculum is also presented.

The prokaryotic Mo/W-bisPGD enzymes family: a catalytic workhorse in bioenergetic

Over the past two decades, prominent importance of molybdenum-containing enzymes in prokaryotes has been put forward by studies originating from different fields. Proteomic or bioinformatic studies underpinned that the list of molybdenum-containing enzymes is far from being complete with to date, more than fifty different enzymes involved in the biogeochemical nitrogen, carbon and sulfur cycles. In particular, the vast majority of prokaryotic molybdenum-containing enzymes belong to the so-called dimethylsulfoxide reductase family. Despite its extraordinary diversity, this family is characterized by the presence of a Mo/W-bis(pyranopterin guanosine dinucleotide) cofactor at the active site. This review highlights what has been learned about the properties of the catalytic site, the modular variation of the structural organization of these enzymes, and their interplay with the isoprenoid quinones. In the last part, this review provides an integrated view of how these enzymes contribute to the bioenergetics of prokaryotes. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Complete genome sequence of Pyrobaculum oguniense

Pyrobaculum oguniense TE7 is an aerobic hyperthermophilic crenarchaeon isolated from a hot spring in Japan. Here we describe its main chromosome of 2,436,033 bp, with three large-scale inversions and an extra-chromosomal element of 16,887 bp. We have annotated 2,800 protein-coding genes and 145 RNA genes in this genome, including nine H/ACA-like small RNA, 83 predicted C/D box small RNA, and 47 transfer RNA genes. Comparative analyses with the closest known relative, the anaerobe Pyrobaculum arsenaticum from Italy, reveals unexpectedly high synteny and nucleotide identity between these two geographically distant species. Deep sequencing of a mixture of genomic DNA from multiple cells has illuminated some of the genome dynamics potentially shared with other species in this genus.

Diversity of Antisense and Other Non-Coding RNAs in Archaea Revealed by Comparative Small RNA Sequencing in Four Pyrobaculum Species

A great diversity of small, non-coding RNA (ncRNA) molecules with roles in gene regulation and RNA processing have been intensely studied in eukaryotic and bacterial model organisms, yet our knowledge of possible parallel roles for small RNAs (sRNA) in archaea is limited. We employed RNA-seq to identify novel sRNA across multiple species of the hyperthermophilic genus Pyrobaculum, known for unusual RNA gene characteristics. By comparing transcriptional data collected in parallel among four species, we were able to identify conserved RNA genes fitting into known and novel families. Among our findings, we highlight three novel cis-antisense sRNAs encoded opposite to key regulatory (ferric uptake regulator), metabolic (triose-phosphate isomerase), and core transcriptional apparatus genes (transcription factor B). We also found a large increase in the number of conserved C/D box sRNA genes over what had been previously recognized; many of these genes are encoded antisense to protein coding genes. The conserved opposition to orthologous genes across the Pyrobaculum genus suggests similarities to other cis-antisense regulatory systems. Furthermore, the genus-specific nature of these sRNAs indicates they are relatively recent, stable adaptations.

Expression, and molecular and enzymatic characterization of Cu-containing nitrite reductase from a marine ammonia-oxidizing gammaproteobacterium, Nitrosococcus oceani

Ammonia-oxidizing bacteria (AOB) remove intracellular nitrite to prevent its toxicity by a nitrifier denitrification pathway involving two denitrifying enzymes, nitrite reductase and nitric oxide reductase. Here, a Cu-containing nitrite reductase from Nitrosococcus oceani strain NS58, a gammaproteobacterial marine AOB, was expressed in Escherichia coli and purified to homogeneity. Sequence homology analysis indicated that the nitrite reductase from N. oceani was phylogenetically closer to its counterparts from denitrifying bacteria than that of the betaproteobacterium Nitrosomonas europaea. The recombinant enzyme was a homotrimer of a 32 kDa subunit molecule. The enzyme was green in the oxidized state with absorption peaks at 455 nm and 575 nm. EPR spectroscopy indicated the presence of type 2 Cu. Molecular activities and the affinity constant for the nitrite were determined to be 1.6×10(3) s(-1) and 52 μM, respectively.

The sub-cellular localization of Sulfolobus DNA replication

Analyses of the DNA replication-associated proteins of hyperthermophilic archaea have yielded considerable insight into the structure and biochemical function of these evolutionarily conserved factors. However, little is known about the regulation and progression of DNA replication in the context of archaeal cells. In the current work, we describe the generation of strains of Sulfolobus solfataricus and Sulfolobus acidocaldarius that allow the incorporation of nucleoside analogues during DNA replication. We employ this technology, in conjunction with immunolocalization analyses of replisomes, to investigate the sub-cellular localization of nascent DNA and replisomes. Our data reveal a peripheral localization of replisomes in the cell. Furthermore, while the two replication forks emerging from any one of the three replication origins in the Sulfolobus chromosome remain in close proximity, the three origin loci are separated.

Complete genome sequence of strain 1860, a crenarchaeon of the genus Pyrobaculum able to grow with various electron acceptors

Strain 1860, a novel member of the genus Pyrobaculum, is a hyperthermophilic organotrophic crenarchaeon growing anaerobically with various electron acceptors. The complete genome sequence reveals genes for several membrane-bound oxidoreductases, the Embden-Meyerhof and Entner-Doudoroff pathways for glucose metabolism, the tricarboxylic acid cycle, the glyoxylate cycle, and the dicarboxylate/4-hydroxybutyrate cycle.

Discovery of Pyrobaculum small RNA families with atypical pseudouridine guide RNA features

In the Eukarya and Archaea, small RNA-guided pseudouridine modification is believed to be an essential step in ribosomal RNA maturation. While readily modeled and identified by computational methods in eukaryotic species, these guide RNAs have not been found in most archaeal genomes. Using high-throughput transcriptome sequencing and comparative genomics, we have identified ten novel small RNA families that appear to function as H/ACA pseudouridylation guide sRNAs, yet surprisingly lack several expected canonical features. The new RNA genes are transcribed and highly conserved across at least six species in the archaeal hyperthermophilic genus Pyrobaculum. The sRNAs exhibit a single hairpin structure interrupted by a conserved kink-turn motif, yet only two of ten families contain the complete canonical structure found in all other H/ACA sRNAs. Half of the sRNAs lack the conserved 3'-terminal ACA sequence, and many contain only a single 3' guide region rather than the canonical 5' and 3' bipartite guides. The predicted sRNA structures contain guide sequences that exhibit strong complementarity to ribosomal RNA or transfer RNA. Most of the predicted targets of pseudouridine modification are structurally equivalent to those known in other species. One sRNA appears capable of guiding pseudouridine modification at positions U54 and U55 in most or all Pyrobaculum tRNAs. We experimentally tested seven predicted pseudouridine modifications in ribosomal RNA, and all but one was confirmed. The structural insights provided by this new set of Pyrobaculum sRNAs will augment existing models and may facilitate the identification and characterization of new guide sRNAs in other archaeal species.

Gene re-annotation in genome of the extremophile Pyrobaculum aerophilum by using bioinformatics methods

In this paper, we re-annotated the genome of Pyrobaculum aerophilum str. IM2, particularly for hypothetical ORFs. The annotation process includes three parts. Firstly and most importantly, 23 new genes, which were missed in the original annotation, are found by combining similarity search and the ab initio gene finding approaches. Among these new genes, five have significant similarities with function-known genes and the rest have significant similarities with hypothetical ORFs contained in other genomes. Secondly, the coding potentials of the 1645 hypothetical ORFs are re-predicted by using 33 Z curve variables combined with Fisher linear discrimination method. With the accuracy being 99.68%, 25 originally annotated hypothetical ORFs are recognized as non-coding by our method. Thirdly, 80 hypothetical ORFs are assigned with potential functions by using similarity search with BLAST program. Re-annotation of the genome will benefit related researches on this hyperthermophilic crenarchaeon. Also, the re-annotation procedure could be taken as a reference for other archaeal genomes. Details of the revised annotation are freely available at http://cobi.uestc.edu.cn/resource/paero/

The complete genome sequence of Thermoproteus tenax: a physiologically versatile member of the Crenarchaeota

Here, we report on the complete genome sequence of the hyperthermophilic Crenarchaeum Thermoproteus tenax (strain Kra1, DSM 2078^T) a type strain of the crenarchaeotal order Thermoproteales. Its circular 1.84-megabase genome harbors no extrachromosomal elements and 2,051 open reading frames are identified, covering 90.6% of the complete sequence, which represents a high coding density. Derived from the gene content, T. tenax is a representative member of the Crenarchaeota. The organism is strictly anaerobic and sulfur-dependent with optimal growth at 86°C and pH 5.6. One particular feature is the great metabolic versatility, which is not accompanied by a distinct increase of genome size or information density as compared to other Crenarchaeota. T. tenax is able to grow chemolithoautotrophically (CO₂/H₂) as well as chemoorganoheterotrophically in presence of various organic substrates. All pathways for synthesizing the 20 proteinogenic amino acids are present. In addition, two presumably complete gene sets for NADH:quinone oxidoreductase (complex I) were identified in the genome and there is evidence that either NADH or reduced ferredoxin might serve as electron donor. Beside the typical archaeal A₀A₁-ATP synthase, a membrane-bound pyrophosphatase is found, which might contribute to energy conservation. Surprisingly, all genes required for dissimilatory sulfate reduction are present, which is confirmed by growth experiments. Mentionable is furthermore, the presence of two proteins (ParA family ATPase, actin-like protein) that might be involved in cell division in Thermoproteales, where the ESCRT system is absent, and of genes involved in genetic competence (DprA, ComF) that is so far unique within Archaea.

Emerging vaccine informatics

Vaccine informatics is an emerging research area that focuses on development and applications of bioinformatics methods that can be used to facilitate every aspect of the preclinical, clinical, and postlicensure vaccine enterprises. Many immunoinformatics algorithms and resources have been developed to predict T- and B-cell immune epitopes for epitope vaccine development and protective immunity analysis. Vaccine protein candidates are predictable in silico from genome sequences using reverse vaccinology. Systematic transcriptomics and proteomics gene expression analyses facilitate rational vaccine design and identification of gene responses that are correlates of protection in vivo. Mathematical simulations have been used to model host-pathogen interactions and improve vaccine production and vaccination protocols. Computational methods have also been used for development of immunization registries or immunization information systems, assessment of vaccine safety and efficacy, and immunization modeling. Computational literature mining and databases effectively process, mine, and store large amounts of vaccine literature and data. Vaccine Ontology (VO) has been initiated to integrate various vaccine data and support automated reasoning.

Uracil-DNA glycosylase of Thermoplasma acidophilum directs long-patch base excision repair, which is promoted by deoxynucleoside triphosphates and ATP/ADP, into short-patch repair

Hydrolytic deamination of cytosine to uracil in DNA is increased in organisms adapted to high temperatures. Hitherto, the uracil base excision repair (BER) pathway has only been described in two archaeons, the crenarchaeon Pyrobaculum aerophilum and the euryarchaeon Archaeoglobus fulgidus, which are hyperthermophiles and use single-nucleotide replacement. In the former the apurinic/apyrimidinic (AP) site intermediate is removed by the sequential action of a 5'-acting AP endonuclease and a 5'-deoxyribose phosphate lyase, whereas in the latter the AP site is primarily removed by a 3'-acting AP lyase, followed by a 3'-phosphodiesterase. We describe here uracil BER by a cell extract of the thermoacidophilic euryarchaeon Thermoplasma acidophilum, which prefers a similar short-patch repair mode as A. fulgidus. Importantly, T. acidophilumcell extract also efficiently executes ATP/ADP-stimulated long-patch BER in the presence of deoxynucleoside triphosphates, with a repair track of ∼15 nucleotides. Supplementation of recombinant uracil-DNA glycosylase (rTaUDG; ORF Ta0477) increased the formation of short-patch at the expense of long-patch repair intermediates, and additional supplementation of recombinant DNA ligase (rTalig; Ta1148) greatly enhanced repair product formation. TaUDG seems to recruit AP-incising and -excising functions to prepare for rapid single-nucleotide insertion and ligation, thus excluding slower and energy-costly long-patch BER.

Molecular signatures for the Crenarchaeota and the Thaumarchaeota

Crenarchaeotes found in mesophilic marine environments were recently placed into a new phylum of Archaea called the Thaumarchaeota. However, very few molecular characteristics of this new phylum are currently known which can be used to distinguish them from the Crenarchaeota. In addition, their relationships to deep-branching archaeal lineages are unclear. We report here detailed analyses of protein sequences from Crenarchaeota and Thaumarchaeota that have identified many conserved signature indels (CSIs) and signature proteins (SPs) (i.e., proteins for which all significant blast hits are from these groups) that are specific for these archaeal groups. Of the identified signatures 6 CSIs and 13 SPs are specific for the Crenarchaeota phylum; 6 CSIs and >250 SPs are uniquely found in various Thaumarchaeota (viz. Cenarchaeum symbiosum, Nitrosopumilus maritimus and a number of uncultured marine crenarchaeotes) and 3 CSIs and ~10 SPs are found in both Thaumarchaeota and Crenarchaeota species. Some of the molecular signatures are also present in Korarchaeum cryptofilum, which forms the independent phylum Korarchaeota. Although some of these molecular signatures suggest a distant shared ancestry between Thaumarchaeota and Crenarchaeota, our identification of large numbers of Thaumarchaeota-specific proteins and their deep branching between the Crenarchaeota and Euryarchaeota phyla in phylogenetic trees shows that they are distinct from both Crenarchaeota and Euryarchaeota in both genetic and phylogenetic terms. These observations support the placement of marine mesophilic archaea into the separate phylum Thaumarchaeota. Additionally, many CSIs and SPs have been found that are specific for different orders within Crenarchaeota (viz. Sulfolobales-3 CSIs and 169 SPs, Thermoproteales-5 CSIs and 25 SPs, Desulfurococcales-4 SPs, and Sulfolobales and Desulfurococcales-2 CSIs and 18 SPs). The signatures described here provide novel means for distinguishing the Crenarchaeota and the Thaumarchaeota and for the classification of related and novel species in different environments. Functional studies on these signature proteins could lead to discovery of novel biochemical properties that are unique to these groups of archaea.

Properties and evolution of an alcohol dehydrogenase from the Crenarchaeota Pyrobaculum aerophilum

The gene encoding a novel alcohol dehydrogenase (ADH) that belongs to the medium chain dehydrogenase/reductase (MDR) superfamily was identified in the hyperthermophilic archaeon, Pyrobaculum aerophilum. The P. aerophilum ADH gene (Pae2687) was over-expressed in Escherichia coli, and the protein (PyAeADHII) was purified to homogeneity and characterized. The PyAeADHII belongs to a medium chain class because its monomer size is 330 residues and even if it is structurally similar to other enzymes belonging to MDR superfamily, it lacks key residues involved in the coordination of the catalytic Zn ion and in the binding of alcoholic substrates typical of other ADHs. Consistently, PyAeADHII does not show activity on a large number of alcohols, aldheydes or ketones. It is active only when alpha-tetralone is used as a substrate. The enzyme has a strict requirement for NADP(H) as the coenzyme and has remarkable thermophilicity, displaying activity at temperatures up to 95 degrees C. The study of the metabolic pathways of P. aerophilum can provide information on the evolution of genes and enzymes and may be crucial for understanding the evolution of eukaryotic cells.

Metagenomes from high-temperature chemotrophic systems reveal geochemical controls on microbial community structure and function

The Yellowstone caldera contains the most numerous and diverse geothermal systems on Earth, yielding an extensive array of unique high-temperature environments that host a variety of deeply-rooted and understudied Archaea, Bacteria and Eukarya. The combination of extreme temperature and chemical conditions encountered in geothermal environments often results in considerably less microbial diversity than other terrestrial habitats and offers a tremendous opportunity for studying the structure and function of indigenous microbial communities and for establishing linkages between putative metabolisms and element cycling. Metagenome sequence (14–15,000 Sanger reads per site) was obtained for five high-temperature (>65°C) chemotrophic microbial communities sampled from geothermal springs (or pools) in Yellowstone National Park (YNP) that exhibit a wide range in geochemistry including pH, dissolved sulfide, dissolved oxygen and ferrous iron. Metagenome data revealed significant differences in the predominant phyla associated with each of these geochemical environments. Novel members of the Sulfolobales are dominant in low pH environments, while other Crenarchaeota including distantly-related Thermoproteales and Desulfurococcales populations dominate in suboxic sulfidic sediments. Several novel archaeal groups are well represented in an acidic (pH 3) Fe-oxyhydroxide mat, where a higher O₂ influx is accompanied with an increase in archaeal diversity. The presence or absence of genes and pathways important in S oxidation-reduction, H₂-oxidation, and aerobic respiration (terminal oxidation) provide insight regarding the metabolic strategies of indigenous organisms present in geothermal systems. Multiple-pathway and protein-specific functional analysis of metagenome sequence data corroborated results from phylogenetic analyses and clearly demonstrate major differences in metabolic potential across sites. The distribution of functional genes involved in electron transport is consistent with the hypothesis that geochemical parameters (e.g., pH, sulfide, Fe, O₂) control microbial community structure and function in YNP geothermal springs.

Biochemical properties and base excision repair complex formation of apurinic/apyrimidinic endonuclease from Pyrococcus furiosus

Apurinic/apyrimidinic (AP) sites are the most frequently found mutagenic lesions in DNA, and they arise mainly from spontaneous base loss or modified base removal by damage-specific DNA glycosylases. AP sites are cleaved by AP endonucleases, and the resultant gaps in the DNA are repaired by DNA polymerase/DNA ligase reactions. We identified the gene product that is responsible for the AP endonuclease activity in the hyperthermophilic euryarchaeon, Pyrococcus furiosus. Furthermore, we detected the physical interaction between P. furiosus AP endonuclease (PfuAPE) and proliferating cell nuclear antigen (PCNA; PfuPCNA) by a pull-down assay and a surface plasmon resonance analysis. Interestingly, the associated 3′–5′ exonuclease activity, but not the AP endonuclease activity, of PfuAPE was stimulated by PfuPCNA. Immunoprecipitation experiments using the P. furiosus cell extracts supported the interaction between PfuAPE and PfuPCNA in the cells. This is the first report describing the physical and functional interactions between an archaeal AP endonuclease and PCNA. We also detected the ternary complex of PfuPCNA, PfuAPE and Pfu uracil-DNA glycosylase. This complex probably functions to enhance the repair of uracil-containing DNA in P. furiosus cells.

Gain and loss of an intron in a protein-coding gene in Archaea: the case of an archaeal RNA pseudouridine synthase gene

Background

We previously found the first examples of splicing of archaeal pre-mRNAs for homologs of the eukaryotic CBF5 protein (also known as dyskerin in humans) in Aeropyrum pernix, Sulfolobus solfataricus, S. tokodaii, and S. acidocaldarirus, and also showed that crenarchaeal species in orders Desulfurococcales and Sulfolobales, except for Hyperthermus butylicus, Pyrodictium occultum, Pyrolobus fumarii, and Ignicoccus islandicus, contain the (putative) cbf5 intron. However, the exact timing of the intron insertion was not determined and verification of the putative secondary loss of the intron in some lineages was not performed.

Results

In the present study, we determined approximately two-thirds of the entire coding region of crenarchaeal Cbf5 sequences from 43 species. A phylogenetic analysis of our data and information from the available genome sequences suggested that the (putative) cbf5 intron existed in the common ancestor of the orders Desulfurococcales and Sulfolobales and that probably at least two independent lineages in the order Desulfurococcales lost the (putative) intron.

Conclusion

This finding is the first observation of a lineage-specific loss of a pre-mRNA intron in Archaea. As the insertion or deletion of introns in protein-coding genes in Archaea has not yet been seriously considered, our finding suggests the possible difficulty of accurately and completely predicting protein-coding genes in Archaea.

Functional importance of crenarchaea-specific extra-loop revealed by an X-ray structure of a heterotetrameric crenarchaeal splicing endonuclease

Archaeal splicing endonucleases (EndAs) are currently classified into three groups. Two groups require a single subunit protein to form a homodimer or homotetramer. The third group requires two nonidentical protein components for the activity. To elucidate the molecular architecture of the two-subunit EndA system, we studied a crenarchaeal splicing endonuclease from Pyrobaculum aerophilum. In the present study, we solved a crystal structure of the enzyme at 1.7-Å resolution. The enzyme adopts a heterotetrameric form composed of two catalytic and two structural subunits. By connecting the structural and the catalytic subunits of the heterotetrameric EndA, we could convert the enzyme to a homodimer that maintains the broad substrate specificity that is one of the characteristics of heterotetrameric EndA. Meanwhile, a deletion of six amino acids in a Crenarchaea-specific loop abolished the endonuclease activity even on a substrate with canonical BHB motif. These results indicate that the subunit architecture is not a major factor responsible for the difference of substrate specificity between single- and two-subunit EndA systems. Rather, the structural basis for the broad substrate specificity is built into the crenarchaeal splicing endonuclease itself.

Comprehensive comparative-genomic analysis of type 2 toxin-antitoxin systems and related mobile stress response systems in prokaryotes

Background

The prokaryotic toxin-antitoxin systems (TAS, also referred to as TA loci) are widespread, mobile two-gene modules that can be viewed as selfish genetic elements because they evolved mechanisms to become addictive for replicons and cells in which they reside, but also possess "normal" cellular functions in various forms of stress response and management of prokaryotic population. Several distinct TAS of type 1, where the toxin is a protein and the antitoxin is an antisense RNA, and numerous, unrelated TAS of type 2, in which both the toxin and the antitoxin are proteins, have been experimentally characterized, and it is suspected that many more remain to be identified.

Results

We report a comprehensive comparative-genomic analysis of Type 2 toxin-antitoxin systems in prokaryotes. Using sensitive methods for distant sequence similarity search, genome context analysis and a new approach for the identification of mobile two-component systems, we identified numerous, previously unnoticed protein families that are homologous to toxins and antitoxins of known type 2 TAS. In addition, we predict 12 new families of toxins and 13 families of antitoxins, and also, predict a TAS or TAS-like activity for several gene modules that were not previously suspected to function in that capacity. In particular, we present indications that the two-gene module that encodes a minimal nucleotidyl transferase and the accompanying HEPN protein, and is extremely abundant in many archaea and bacteria, especially, thermophiles might comprise a novel TAS. We present a survey of previously known and newly predicted TAS in 750 complete genomes of archaea and bacteria, quantitatively demonstrate the exceptional mobility of the TAS, and explore the network of toxin-antitoxin pairings that combines plasticity with selectivity.

Conclusion

The defining properties of the TAS, namely, the typically small size of the toxin and antitoxin genes, fast evolution, and extensive horizontal mobility, make the task of comprehensive identification of these systems particularly challenging. However, these same properties can be exploited to develop context-based computational approaches which, combined with exhaustive analysis of subtle sequence similarities were employed in this work to substantially expand the current collection of TAS by predicting both previously unnoticed, derived versions of known toxins and antitoxins, and putative novel TAS-like systems. In a broader context, the TAS belong to the resistome domain of the prokaryotic mobilome which includes partially selfish, addictive gene cassettes involved in various aspects of stress response and organized under the same general principles as the TAS. The "selfish altruism", or "responsible selfishness", of TAS-like systems appears to be a defining feature of the resistome and an important characteristic of the entire prokaryotic pan-genome given that in the prokaryotic world the mobilome and the "stable" chromosomes form a dynamic continuum.

Reviewers

This paper was reviewed by Kenn Gerdes (nominated by Arcady Mushegian), Daniel Haft, Arcady Mushegian, and Andrei Osterman. For full reviews, go to the Reviewers' Reports section.