Orphan SelD proteins and selenium-dependent molybdenum hydroxylases

Clostridioides difficile exploits xanthine and uric acid as nutrients by utilizing a selenium-dependent catabolic pathway

Selenium is a trace element that plays critical roles in redox biology; it is typically incorporated into "selenoproteins" as the 21st amino acid selenocysteine. Additionally, selenium exists as a labile non-selenocysteine cofactor in a small subset of selenoproteins known as selenium-dependent molybdenum hydroxylases (SDMHs). In purinolytic clostridia, SDMHs are implicated in the degradation of hypoxanthine, xanthine, and uric acid for carbon and nitrogen. While SDMHs have been biochemically analyzed, the genes responsible for the insertion and maturation of the selenium cofactor lack characterization. In this study, we utilized the nosocomial pathogen Clostridioides difficile as a genetic model to begin characterizing this poorly understood selenium utilization pathway and its role in the catabolism of host-derived purines. We first observed that C. difficile could utilize hypoxanthine, xanthine, or uric acid to overcome a growth defect in a minimal medium devoid of glycine and threonine. However, strains lacking selenophosphate synthetase (selD mutants) still grew poorly in the presence of xanthine and uric acid, suggesting a selenium-dependent purinolytic process. Previous computational studies have identified yqeB and yqeC as potential candidates for cofactor maturation, so we subsequently deleted each gene using CRISPR-Cas9 technology. We surprisingly found that the growth of the ΔyqeB mutant in response to each purine was similar to the behavior of the selD mutants, while the ΔyqeC mutant exhibited no obvious phenotype. Our results suggest an important role for YqeB in selenium-dependent purine catabolism and also showcase C. difficile as an appropriate model organism to study the biological use of selenium.IMPORTANCEThe apparent modification of bacterial molybdenum hydroxylases with a catalytically essential selenium cofactor is the least understood mechanism of selenium incorporation. Selenium-dependent molybdenum hydroxylases play an important role in scavenging carbon and nitrogen from purines for purinolytic clostridia. Here, we used Clostridioides difficile as a genetic platform to begin dissecting the selenium cofactor trait and found genetic evidence for a selenium-dependent purinolytic pathway. The absence of selD or yqeB-a predicted genetic marker for the selenium cofactor trait-resulted in impaired growth on xanthine and uric acid, known substrates for selenium-dependent molybdenum hydroxylases. Our findings provide a genetic foundation for future research of this pathway and suggest a novel metabolic strategy for C. difficile to scavenge host-derived purines from the gut.

Enzymatic strategies for selenium incorporation into biological molecules

The trace element selenium (Se) is essential to the physiology of most organisms on the planet. The most well documented of Se's biological forms are selenoproteins, where selenocysteine often serves as the catalytic center for crucial redox processes. Se is also found in several other classes of biological molecules, including nucleic acids, sugars, and modified amino acids, although its role in the function of these metabolites is less understood. Despite its prevalence, only a small number of Se-specific biosynthetic pathways have been discovered. Around half of these were first characterized in the past three years, suggesting that the selenometabolome may be more diverse than previously appreciated. Here, we review the recent advances in our understanding of this intriguing biochemical space, and discuss prospects for future discovery efforts.

Proteomic and transcriptomic analysis of selenium utilization in Methanococcus maripaludis

Methanococcus maripaludis utilizes selenocysteine- (Sec-) containing proteins (selenoproteins), mostly active in the organism's primary energy metabolism, methanogenesis. During selenium depletion, M. maripaludis employs a set of enzymes containing cysteine (Cys) instead of Sec. The genes coding for these Sec-/Cys-containing isoforms were the only genes known of which expression is influenced by the selenium status of the cell. Using proteomics and transcriptomics, approx. 7% and 12%, respectively, of all genes/proteins were found differentially expressed/synthesized in response to the selenium supply. Some of the genes identified involve methanogenesis, nitrogenase functions, and putative transporters. An increase of transcript abundance for putative transporters under selenium depletion indicated the organism's effort to tap into alternative sources of selenium. M. maripaludis is known to utilize selenite and dimethylselenide as selenium sources. To expand this list, a selenium-responsive reporter strain was assessed with nine other, environmentally relevant selenium species. While the effect of some was very similar to that of selenite, others were effectively utilized at lower concentrations. Conversely, selenate and seleno-amino acids were only utilized at unphysiologically high concentrations and two compounds were not utilized at all. To address the role of the selenium-regulated putative transporters, M. maripaludis mutant strains lacking one or two of the putative transporters were tested for the capability to utilize the different selenium species. Of the five putative transporters analyzed by loss-of-function mutagenesis, none appeared to be absolutely required for utilizing any of the selenium species tested, indicating they have redundant and/or overlapping specificities or are not dedicated selenium transporters.

IMPORTANCE

While selenium metabolism in microorganisms has been studied intensively in the past, global gene expression approaches have not been employed so far. Furthermore, the use of different selenium sources, widely environmentally interconvertible via biotic and abiotic processes, was also not extensively studied before. Methanococcus maripaludis JJ is ideally suited for such analyses, thanks to its known selenium usage and available genetic tools. Thus, an overall view on the selenium regulon of M. maripaludis was obtained via transcriptomic and proteomic analyses, which inspired further experimentation. This led to demonstrating the use of selenium sources M. maripaludis was previously not known to employ. Also, an attempt-although so far unsuccessful-was made to pinpoint potential selenium transporter genes, in order to deepen our understanding of trace element utilization in this important model organism.

Gut bacterial metabolism contributes to host global purine homeostasis

The microbes and microbial pathways that influence host inflammatory disease progression remain largely undefined. Here, we show that variation in atherosclerosis burden is partially driven by gut microbiota and is associated with circulating levels of uric acid (UA) in mice and humans. We identify gut bacterial taxa spanning multiple phyla, including Bacillota, Fusobacteriota, and Pseudomonadota, that use multiple purines, including UA as carbon and energy sources anaerobically. We identify a gene cluster that encodes key steps of anaerobic purine degradation and that is widely distributed among gut-dwelling bacteria. Furthermore, we show that colonization of gnotobiotic mice with purine-degrading bacteria modulates levels of UA and other purines in the gut and systemically. Thus, gut microbes are important drivers of host global purine homeostasis and serum UA levels, and gut bacterial catabolism of purines may represent a mechanism by which gut bacteria influence health.

Purine catabolism by enterobacteria

Purines are abundant among organic nitrogen sources and have high nitrogen content. Accordingly, microorganisms have evolved different pathways to catabolize purines and their metabolic products such as allantoin. Enterobacteria from the genera Escherichia, Klebsiella and Salmonella have three such pathways. First, the HPX pathway, found in the genus Klebsiella and very close relatives, catabolizes purines during aerobic growth, extracting all four nitrogen atoms in the process. This pathway includes several known or predicted enzymes not previously observed in other purine catabolic pathways. Second, the ALL pathway, found in strains from all three species, catabolizes allantoin during anaerobic growth in a branched pathway that also includes glyoxylate assimilation. This allantoin fermentation pathway originally was characterized in a gram-positive bacterium, and therefore is widespread. Third, the XDH pathway, found in strains from Escherichia and Klebsiella spp., at present is ill-defined but likely includes enzymes to catabolize purines during anaerobic growth. Critically, this pathway may include an enzyme system for anaerobic urate catabolism, a phenomenon not previously described. Documenting such a pathway would overturn the long-held assumption that urate catabolism requires oxygen. Overall, this broad capability for purine catabolism during either aerobic or anaerobic growth suggests that purines and their metabolites contribute to enterobacterial fitness in a variety of environments.

Eight Unexpected Selenoprotein Families in Organometallic Biochemistry in Clostridium difficile, in ABC Transport, and in Methylmercury Biosynthesis

The bioinformatics of a nine-gene locus, designated selenocysteine-assisted organometallic (SAO), was investigated after identifying six new selenoprotein families and constructing hidden Markov models (HMMs) that find and annotate members of those families. Four are selenoproteins in most SAO loci, including Clostridium difficile. They include two ABC transporter subunits, namely, permease SaoP, with selenocysteine (U) at the channel-gating position, and substrate-binding subunit SaoB. Cytosolic selenoproteins include SaoL, homologous to MerB organomercurial lyases from mercury resistance loci, and SaoT, related to thioredoxins. SaoL, SaoB, and surface protein SaoC (an occasional selenoprotein) share an unusual CU dipeptide motif, which is something rare in selenoproteins but found in selenoprotein variants of mercury resistance transporter subunit MerT. A nonselenoprotein, SaoE, shares homology with Cu/Zn efflux and arsenical efflux pumps. The organization of the SAO system suggests substrate interaction with surface-exposed selenoproteins, followed by import, metabolism that may cleave a carbon-to-heavy metal bond, and finally metal efflux. A novel type of mercury resistance is possible, but SAO instead may support fermentative metabolism, with selenocysteine-mediated formation of organometallic intermediates, followed by import, degradation, and metal efflux. Phylogenetic profiling shows SOA loci consistently co-occur with Stickland fermentation markers but even more consistently with 8Fe-9S cofactor-type double-cubane proteins. Hypothesizing that the SAO system forms organometallic intermediates, we investigated the known methylmercury formation protein families HgcA and HgcB. Both families contained overlooked selenoproteins. Most HgcAs have a CU motif N terminal to their previously accepted start sites. Seeking additional rare and overlooked selenoproteins may help reveal more cryptic aspects of microbial biochemistry. IMPORTANCE This work adds 8 novel prokaryotic selenoproteins to the 80 or so families previously known. It describes the SAO (selenocysteine-assisted organometallic) locus, with the most selenoproteins of any known system. The rare CU motif recurs throughout, suggesting the formation and degradation of organometallic compounds. That suggestion triggered a reexamination of HgcA and HcgB, which are methylmercury formation proteins that can adversely impact food safety. Both are selenoproteins, once corrected, with HgcA again showing a CU motif. The SAO system is plausibly a mercury resistance locus for selenium-dependent anaerobes. But instead, it may exploit heavy metals as cofactors in organometallic intermediate-forming pathways that circumvent high activation energies and facilitate the breakdown of otherwise poorly accessible nutrients. SAO could provide an edge that helps Clostridium difficile, an important pathogen, establish disease.

The selenophosphate synthetase family: A review

Selenophosphate synthetases use selenium and ATP to synthesize selenophosphate. This is required for biological utilization of selenium, most notably for the synthesis of the non-canonical amino acid selenocysteine (Sec). Therefore, selenophosphate synthetases underlie all functions of selenoproteins, which include redox homeostasis, protein quality control, hormone regulation, metabolism, and many others. This protein family comprises two groups, SelD/SPS2 and SPS1. The SelD/SPS2 group represent true selenophosphate synthetases, enzymes central to selenium metabolism which are present in all Sec-utilizing organisms across the tree of life. Notably, many SelD/SPS2 proteins contain Sec as catalytic residue in their N-terminal flexible selenium-binding loop, while others replace it with cysteine (Cys). The SPS1 group comprises proteins originated through gene duplications of SelD/SPS2 in metazoa in which the Sec/Cys-dependent catalysis was disrupted. SPS1 proteins do not synthesize selenophosphate and are not required for Sec synthesis. They have essential regulatory functions related to redox homeostasis and pyridoxal phosphate, which affect signaling pathways for growth and differentiation. In this review, we summarize the knowledge about the selenophosphate synthetase family acquired through decades of research, encompassing their structure, mechanism, function, and evolution.

Selenium Metabolism and Selenoproteins in Prokaryotes: A Bioinformatics Perspective

Selenium (Se) is an important trace element that mainly occurs in the form of selenocysteine in selected proteins. In prokaryotes, Se is also required for the synthesis of selenouridine and Se-containing cofactor. A large number of selenoprotein families have been identified in diverse prokaryotic organisms, most of which are thought to be involved in various redox reactions. In the last decade or two, computational prediction of selenoprotein genes and comparative genomics of Se metabolic pathways and selenoproteomes have arisen, providing new insights into the metabolism and function of Se and their evolutionary trends in bacteria and archaea. This review aims to offer an overview of recent advances in bioinformatics analysis of Se utilization in prokaryotes. We describe current computational strategies for the identification of selenoprotein genes and generate the most comprehensive list of prokaryotic selenoproteins reported to date. Furthermore, we highlight the latest research progress in comparative genomics and metagenomics of Se utilization in prokaryotes, which demonstrates the divergent and dynamic evolutionary patterns of different Se metabolic pathways, selenoprotein families, and selenoproteomes in sequenced organisms and environmental samples. Overall, bioinformatics analyses of Se utilization, function, and evolution may contribute to a systematic understanding of how this micronutrient is used in nature.

Historical Roles of Selenium and Selenoproteins in Health and Development: The Good, the Bad and the Ugly

Selenium is a fascinating element that has a long history, most of which documents it as a deleterious element to health. In more recent years, selenium has been found to be an essential element in the diet of humans, all other mammals, and many other life forms. It has many health benefits that include, for example, roles in preventing heart disease and certain forms of cancer, slowing AIDS progression in HIV patients, supporting male reproduction, inhibiting viral expression, and boosting the immune system, and it also plays essential roles in mammalian development. Elucidating the molecular biology of selenium over the past 40 years generated an entirely new field of science which encompassed the many novel features of selenium. These features were (1) how this element makes its way into protein as the 21st amino acid in the genetic code, selenocysteine (Sec); (2) the vast amount of machinery dedicated to synthesizing Sec uniquely on its tRNA; (3) the incorporation of Sec into protein; and (4) the roles of the resulting Sec-containing proteins (selenoproteins) in health and development. One of the research areas receiving the most attention regarding selenium in health has been its role in cancer prevention, but further research has also exposed the role of this element as a facilitator of various maladies, including cancer.

Macrophage Selenoproteins Restrict Intracellular Replication of Francisella tularensis and Are Essential for Host Immunity

The essential micronutrient Selenium (Se) is co-translationally incorporated as selenocysteine into proteins. Selenoproteins contain one or more selenocysteines and are vital for optimum immunity. Interestingly, many pathogenic bacteria utilize Se for various biological processes suggesting that Se may play a role in bacterial pathogenesis. A previous study had speculated that Francisella tularensis, a facultative intracellular bacterium and the causative agent of tularemia, sequesters Se by upregulating Se-metabolism genes in type II alveolar epithelial cells. Therefore, we investigated the contribution of host vs. pathogen-associated selenoproteins in bacterial disease using F. tularensis as a model organism. We found that F. tularensis was devoid of any Se utilization traits, neither incorporated elemental Se, nor exhibited Se-dependent growth. However, 100% of Se-deficient mice (0.01 ppm Se), which express low levels of selenoproteins, succumbed to F. tularensis-live vaccine strain pulmonary challenge, whereas 50% of mice on Se-supplemented (0.4 ppm Se) and 25% of mice on Se-adequate (0.1 ppm Se) diet succumbed to infection. Median survival time for Se-deficient mice was 8 days post-infection while Se-supplemented and -adequate mice was 11.5 and >14 days post-infection, respectively. Se-deficient macrophages permitted significantly higher intracellular bacterial replication than Se-supplemented macrophages ex vivo, corroborating in vivo observations. Since Francisella replicates in alveolar macrophages during the acute phase of pneumonic infection, we hypothesized that macrophage-specific host selenoproteins may restrict replication and systemic spread of bacteria. F. tularensis infection led to an increased expression of several macrophage selenoproteins, suggesting their key role in limiting bacterial replication. Upon challenge with F. tularensis, mice lacking selenoproteins in macrophages (TrspM) displayed lower survival and increased bacterial burden in the lung and systemic tissues in comparison to WT littermate controls. Furthermore, macrophages from TrspM mice were unable to restrict bacterial replication ex vivo in comparison to macrophages from littermate controls. We herein describe a novel function of host macrophage-specific selenoproteins in restriction of intracellular bacterial replication. These data suggest that host selenoproteins may be considered as novel targets for modulating immune response to control a bacterial infection.

Functions of elements in soil microorganisms

The soil microbial community fulfils various functions, such as nutrient cycling and carbon (C) sequestration, therefore contributing to maintenance of soil fertility and mitigation of global warming. In this context, a major focus of research has been on C, nitrogen (N) and phosphorus (P) cycling. However, from aquatic and other environments, it is well known that other elements beyond C, N, and P are essential for microbial functioning. Nonetheless, for soil microorganisms this knowledge has not yet been synthesised. To gain a better mechanistic understanding of microbial processes in soil systems, we aimed at summarising the current knowledge on the function of a range of essential or beneficial elements, which may affect the efficiency of microbial processes in soil. This knowledge is discussed in the context of microbial driven nutrient and C cycling. Our findings may support future investigations and data evaluation, where other elements than C, N, and P affect microbial processes.

The physiology and evolution of microbial selenium metabolism

Selenium is an essential trace element whose compounds are widely metabolized by organisms from all three domains of life. Moreover, phylogenetic evidence indicates that selenium species, along with iron, molybdenum, tungsten, and nickel, were metabolized by the last universal common ancestor of all cellular lineages, primarily for the synthesis of the 21st amino acid selenocysteine. Thus, selenium metabolism is both environmentally ubiquitous and a physiological adaptation of primordial life. Selenium metabolic reactions comprise reductive transformations both for assimilation into macromolecules and dissimilatory reduction of selenium oxyanions and elemental selenium during anaerobic respiration. This review offers a comprehensive overview of the physiology and evolution of both assimilatory and dissimilatory selenium metabolism in bacteria and archaea, highlighting mechanisms of selenium respiration. This includes a thorough discussion of our current knowledge of the physiology of selenocysteine synthesis and incorporation into proteins in bacteria obtained from structural biology. Additionally, this is the first comprehensive discussion in a review of the incorporation of selenium into the tRNA nucleoside 5-methylaminomethyl-2-selenouridine and as an inorganic cofactor in certain molybdenum hydroxylase enzymes. Throughout, conserved mechanisms and derived features of selenium metabolism in both domains are emphasized and discussed within the context of the global selenium biogeochemical cycle.

Exploring the selenium-over-sulfur substrate specificity and kinetics of a bacterial selenocysteine lyase

Selenium is a vital micronutrient in many organisms. While traces are required for microbial utilization, excess amounts are toxic; thus, selenium can be regarded as a biological double-edged sword. Selenium is chemically similar to the essential element sulfur, but curiously, evolution has selected the former over the latter for a subset of oxidoreductases. Enzymes involved in sulfur metabolism are less discriminate in terms of preventing selenium incorporation; however, its specific incorporation into selenoproteins reveals a highly discriminate process that is not completely understood. We have identified SclA, a NifS-like protein in the nosocomial pathogen, Enterococcus faecalis, and characterized its enzymatic activity and specificity for l-selenocysteine over l-cysteine. It is known that Asp-146 is required for selenocysteine specificity in the human selenocysteine lyase. Thus, using computational biology, we compared the bacterial and mammalian enzymes and identified His-100, an Asp-146 ortholog in SclA, and generated site-directed mutants in order to study the residue's potential role in the l-selenocysteine discrimination mechanism. The proteins were overexpressed, purified, and characterized for their biochemical properties. All mutants exhibited varying Michaelis-Menten behavior towards l-selenocysteine, but His-100 was not found to be essential for this activity. Additionally, l-cysteine acted as a competitive inhibitor of all enzymes with higher affinity than l-selenocysteine. Finally, we discovered that SclA exhibited low activity with l-cysteine as a poor substrate regardless of mutations. We conclude that His-100 is not required for l-selenocysteine specificity, underscoring the inherent differences in discriminatory mechanisms between bacterial NifS-like proteins and mammalian selenocysteine lyases.

Genes Contributing to the Unique Biology and Intrinsic Antibiotic Resistance of Enterococcus faecalis

Enterococci are leading causes of antibiotic-resistant infection transmitted in hospitals. The intrinsic hardiness of these organisms allows them to survive disinfection practices and then proliferate in the gastrointestinal tracts of antibiotic-treated patients. The objective of this study was to identify the underlying genetic basis for its unusual hardiness. Using a functional genomic approach, we identified traits and pathways of general importance for enterococcal survival and growth that distinguish them from closely related pathogens as well as ancestrally related species. We further identified unique traits that enable them to survive antibiotic challenge, revealing a large set of genes that contribute to intrinsic antibiotic resistance and a smaller set of uniquely important genes that are rare outside enterococci.

The enterococci, which are among the leading causes of multidrug-resistant (MDR) hospital infection, are notable for their environmental ruggedness, which extends to intrinsic antibiotic resistance. To identify genes that confer this unique property, we used Tn-seq to comprehensively explore the genome of MDR Enterococcus faecalis strain MMH594 for genes important for growth in nutrient-containing medium and with low-level antibiotic challenge. As expected, a large core of genes for DNA replication, expression, and central metabolism, shared with other bacteria, are intolerant to transposon disruption. However, genes were identified that are important to E. faecalis that are either absent from or unimportant for Staphylococcus aureus and Streptococcus pneumoniae fitness when similarly tested. Further, 217 genes were identified that when challenged by sub-MIC antibiotic levels exhibited reduced tolerance to transposon disruption, including those previously shown to contribute to intrinsic resistance, and others not previously ascribed this role. E. faecalis is one of the few Gram-positive bacteria experimentally shown to possess a functional Entner-Doudoroff pathway for carbon metabolism, a pathway that contributes to stress tolerance in other microbes. Through functional genomics and network analysis we defined the unusual structure of this pathway in E. faecalis and assessed its importance. These approaches also identified toxin-antitoxin and related systems that are unique and active in E. faecalis. Finally, we identified genes that are absent in the closest nonenterococcal relatives, the vagococci, and that contribute importantly to fitness with and without antibiotic selection, advancing an understanding of the unique biology of enterococci.

Selenium nanovirus and its cytotoxicity in selenite-exposed higher living organisms

Selenium (Se) is an essential micronutrient in living organisms, having a narrow margin between essential and potentially toxic intake/exposure. Thus, the biochemistry of Se in living organisms must be studied in-depth to determine the underlying mechanism of Se cytotoxicity. In this study, we report the emergence of selenium nanovirus (SeNVs) in selenite-exposed fish (freshwater and saltwater) and plants (dryland) and its toxicity in them. SeNVs were found in both the abdomen and tail of Oryzias melastigma and saltwater Rhodeus ocellatus, which led to their death. The occurrence of the intracellular assembly of SeNVs was observed in the roots and leaves of corn Zea mays, but not in those of Limnobium laevigatum. SeNVs led to the death of Z. mays but caused chronic toxicity in L. laevigatum. SeNVs should be a system or structure that dissipates the intracellular redox gradients of the host cells, with simple information consisting Se-O, Se-N, or Se-S bond, that would ensure elemental Se ligand binding with nearly specific biomolecules in host cells, thereby maintaining their composition and stabilizing their structure. The multiple toxic effects of Se, therefore, could be the consequence of increase of entropy in the host cells caused by the intracellular assembly of SeNVs. This study may provide an insight into the underlying mechanism of Se in environmental toxicology and its applications in human health.

Comparative genomic analysis of selenium utilization traits in different marine environments

Selenium (Se) is an essential trace element for many organisms, which is required in the biosynthesis of proteins with selenocysteine, tRNAs with selenouridine, and certain enzymes with Se as a cofactor. Recent large-scale metagenomics projects provide a unique opportunity for studying the global trends of Se utilization in marine environments. Here, we analyzed samples from different marine microbial communities, revealed by the Tara Oceans project, to characterize the Se utilization traits. We found that the selenophosphate synthetase gene, which defines the overall Se utilization, and Se utilization traits are present in all samples. Regions with samples rich and poor in Se utilization traits were categorized. From the analysis of environmental factors, the mesopelagic zone and high temperature (> 15°C) of water are favorable, while geographical location has little influence on Se utilization. All Se utilization traits showed a relatively independent occurrence. The taxonomic classification of Se traits shows that most of the sequences corresponding to Se utilization traits belong to the phylum Proteobacteria. Overall, our study provides useful insights into the general features of Se utilization in ocean samples and may help to understand the evolutionary dynamics of Se utilization in different marine environments.

The complete genome sequence of Eubacterium limosum SA11, a metabolically versatile rumen acetogen

Acetogens are a specialized group of anaerobic bacteria able to produce acetate from CO2 and H2 via the Wood-Ljungdahl pathway. In some gut environments acetogens can compete with methanogens for H2, and as a result rumen acetogens are of interest in the development of microbial approaches for methane mitigation. The acetogen Eubacterium limosum SA11 was isolated from the rumen of a New Zealand sheep and its genome has been sequenced to examine its potential application in methane mitigation strategies, particularly in situations where hydrogenotrophic methanogens are inhibited resulting in increased H2 levels in the rumen. The 4.15 Mb chromosome of SA11 has an average G + C content of 47 %, and encodes 3805 protein-coding genes. There is a single prophage inserted in the chromosome, and several other gene clusters appear to have been acquired by horizontal transfer. These include genes for cell wall glycopolymers, a type VII secretion system, cell surface proteins and chemotaxis. SA11 is able to use a variety of organic substrates in addition to H2/CO2, with acetate and butyrate as the principal fermentation end-products, and genes involved in these metabolic pathways have been identified. An unusual feature is the presence of 39 genes encoding trimethylamine methyltransferase family proteins, more than any other bacterial genome. Overall, SA11 is a metabolically versatile organism, but its ability to grow on such a wide range of substrates suggests it may not be a suitable candidate to take the place of hydrogen-utilizing methanogens in the rumen.

Comparative genomics reveals new evolutionary and ecological patterns of selenium utilization in bacteria

Selenium (Se) is an important micronutrient for many organisms, which is required for the biosynthesis of selenocysteine, selenouridine and Se-containing cofactor. Several key genes involved in different Se utilization traits have been characterized; however, systematic studies on the evolution and ecological niches of Se utilization are very limited. Here, we analyzed more than 5200 sequenced organisms to examine the occurrence patterns of all Se traits in bacteria. A global species map of all Se utilization pathways has been generated, which demonstrates the most detailed understanding of Se utilization in bacteria so far. In addition, the selenophosphate synthetase gene, which is used to define the overall Se utilization, was also detected in some organisms that do not have any of the known Se traits, implying the presence of a novel Se form in this domain. Phylogenetic analyses of components of different Se utilization traits revealed new horizontal gene transfer events for each of them. Moreover, by characterizing the selenoproteomes of all organisms, we found a new selenoprotein-rich phylum and additional selenoprotein-rich species. Finally, the relationship between ecological environments and Se utilization was investigated and further verified by metagenomic analysis of environmental samples, which indicates new macroevolutionary trends of each Se utilization trait in bacteria. Our data provide insights into the general features of Se utilization in bacteria and should be useful for a further understanding of the evolutionary dynamics of Se utilization in nature.

Genome-Scale Mapping of Escherichia coli σ54 Reveals Widespread, Conserved Intragenic Binding

Bacterial RNA polymerases must associate with a σ factor to bind promoter DNA and initiate transcription. There are two families of σ factor: the σ⁷⁰ family and the σ⁵⁴ family. Members of the σ⁵⁴ family are distinct in their ability to bind promoter DNA sequences, in the context of RNA polymerase holoenzyme, in a transcriptionally inactive state. Here, we map the genome-wide association of Escherichia coli σ⁵⁴, the archetypal member of the σ⁵⁴ family. Thus, we vastly expand the list of known σ⁵⁴ binding sites to 135. Moreover, we estimate that there are more than 250 σ⁵⁴ sites in total. Strikingly, the majority of σ⁵⁴ binding sites are located inside genes. The location and orientation of intragenic σ⁵⁴ binding sites is non-random, and many intragenic σ⁵⁴ binding sites are conserved. We conclude that many intragenic σ⁵⁴ binding sites are likely to be functional. Consistent with this assertion, we identify three conserved, intragenic σ⁵⁴ promoters that drive transcription of mRNAs with unusually long 5ʹ UTRs.

Bacterial RNA polymerases must associate with a σ factor to bind to promoter DNA sequences upstream of genes and initiate transcription. There are two families of σ factor: σ⁷⁰ and σ⁵⁴. Members of the σ⁵⁴ family are distinct from members of the σ⁷⁰ family in their ability to bind promoter DNA sequences, in association with RNA polymerase, in a transcriptionally inactive state. We have determined positions in the Escherichia coli genome that are bound by σ⁵⁴, the archetypal member of the σ⁵⁴ family. Surprisingly, we identified 135 binding sites for σ⁵⁴, a huge increase over the number of previously described sites. Our data suggest that there are more than 250 σ⁵⁴ sites in total. Strikingly, most σ⁵⁴ binding sites are located inside genes, whereas only one intragenic σ⁵⁴ binding site has previously been described. The location and orientation of intragenic σ⁵⁴ binding sites is non-random, and many intragenic σ⁵⁴ binding sites are conserved in other bacterial species. We conclude that many intragenic σ⁵⁴ binding sites are likely to be functional. Consistent with this notion, we identify three σ⁵⁴ promoters in E. coli that are located inside genes but drive transcription of unusual mRNAs for the neighboring genes.

Evolution of selenophosphate synthetases: emergence and relocation of function through independent duplications and recurrent subfunctionalization

Selenoproteins are proteins that incorporate selenocysteine (Sec), a nonstandard amino acid encoded by UGA, normally a stop codon. Sec synthesis requires the enzyme Selenophosphate synthetase (SPS or SelD), conserved in all prokaryotic and eukaryotic genomes encoding selenoproteins. Here, we study the evolutionary history of SPS genes, providing a map of selenoprotein function spanning the whole tree of life. SPS is itself a selenoprotein in many species, although functionally equivalent homologs that replace the Sec site with cysteine (Cys) are common. Many metazoans, however, possess SPS genes with substitutions other than Sec or Cys (collectively referred to as SPS1). Using complementation assays in fly mutants, we show that these genes share a common function, which appears to be distinct from the synthesis of selenophosphate carried out by the Sec- and Cys- SPS genes (termed SPS2), and unrelated to Sec synthesis. We show here that SPS1 genes originated through a number of independent gene duplications from an ancestral metazoan selenoprotein SPS2 gene that most likely already carried the SPS1 function. Thus, in SPS genes, parallel duplications and subsequent convergent subfunctionalization have resulted in the segregation to different loci of functions initially carried by a single gene. This evolutionary history constitutes a remarkable example of emergence and evolution of gene function, which we have been able to trace thanks to the singular features of SPS genes, wherein the amino acid at a single site determines unequivocally protein function and is intertwined to the evolutionary fate of the entire selenoproteome.

Comparative genomics reveals new candidate genes involved in selenium metabolism in prokaryotes

Selenium (Se) is an important micronutrient that mainly occurs in proteins in the form of selenocysteine and in tRNAs in the form of selenouridine. In the past 20 years, several genes involved in Se utilization have been characterized in both prokaryotes and eukaryotes. However, Se homeostasis and the associated regulatory network are not fully understood. In this study, we conducted comparative genomics and phylogenetic analyses to examine the occurrence of all known Se utilization traits in prokaryotes. Our results revealed a highly mosaic pattern of species that use Se (in different forms) in spite that most organisms do not use this element. Further investigation of genomic context of known Se-related genes in different organisms suggested novel candidate genes that may participate in Se metabolism in bacteria and/or archaea. Among them, a membrane protein, YedE, which contains ten transmembrane domains and shows distant similarity to a sulfur transporter, is exclusively found in Se-utilizing organisms, suggesting that it may be involved in Se transport. A LysR-like transcription factor subfamily might be important for the regulation of Sec biosynthesis and/or other Se-related genes. In addition, a small protein family DUF3343 is widespread in Se-utilizing organisms, which probably serves as an important chaperone for Se trafficking within the cells. Finally, we proposed a simple model of Se homeostasis based on our findings. Our study reveals new candidate genes involved in Se metabolism in prokaryotes and should be useful for a further understanding of the complex metabolism and the roles of Se in biology.

Computational identification of a new SelD-like family that may participate in sulfur metabolism in hyperthermophilic sulfur-reducing archaea

BACKGROUND

Selenium (Se) and sulfur (S) are closely related elements that exhibit similar chemical properties. Some genes related to S metabolism are also involved in Se utilization in many organisms. However, the evolutionary relationship between the two utilization traits is unclear.

RESULTS

In this study, we conducted a comparative analysis of the selenophosphate synthetase (SelD) family, a key protein for all known Se utilization traits, in all sequenced archaea. Our search showed a very limited distribution of SelD and Se utilization in this kingdom. Interestingly, a SelD-like protein was detected in two orders of Crenarchaeota: Sulfolobales and Thermoproteales. Sequence and phylogenetic analyses revealed that SelD-like protein contains the same domain and conserved functional residues as those of SelD, and might be involved in S metabolism in these S-reducing organisms. Further genome-wide analysis of patterns of gene occurrence in different thermoproteales suggested that several genes, including SirA-like, Prx-like and adenylylsulfate reductase, were strongly related to SelD-like gene. Based on these findings, we proposed a simple model wherein SelD-like may play an important role in the biosynthesis of certain thiophosphate compound.

CONCLUSIONS

Our data suggest novel genes involved in S metabolism in hyperthermophilic S-reducing archaea, and may provide a new window for understanding the complex relationship between Se and S metabolism in archaea.

The genome of Spironucleus salmonicida highlights a fish pathogen adapted to fluctuating environments

Spironucleus salmonicida causes systemic infections in salmonid fish. It belongs to the group diplomonads, binucleated heterotrophic flagellates adapted to micro-aerobic environments. Recently we identified energy-producing hydrogenosomes in S. salmonicida. Here we present a genome analysis of the fish parasite with a focus on the comparison to the more studied diplomonad Giardia intestinalis. We annotated 8067 protein coding genes in the ∼12.9 Mbp S. salmonicida genome. Unlike G. intestinalis, promoter-like motifs were found upstream of genes which are correlated with gene expression, suggesting a more elaborate transcriptional regulation. S. salmonicida can utilise more carbohydrates as energy sources, has an extended amino acid and sulfur metabolism, and more enzymes involved in scavenging of reactive oxygen species compared to G. intestinalis. Both genomes have large families of cysteine-rich membrane proteins. A cluster analysis indicated large divergence of these families in the two diplomonads. Nevertheless, one of S. salmonicida cysteine-rich proteins was localised to the plasma membrane similar to G. intestinalis variant-surface proteins. We identified S. salmonicida homologs to cyst wall proteins and showed that one of these is functional when expressed in Giardia. This suggests that the fish parasite is transmitted as a cyst between hosts. The extended metabolic repertoire and more extensive gene regulation compared to G. intestinalis suggest that the fish parasite is more adapted to cope with environmental fluctuations. Our genome analyses indicate that S. salmonicida is a well-adapted pathogen that can colonize different sites in the host.

Studies of model organisms are very powerful. However, to appreciate the enormous diversity of genetic and cell biological processes we need to extend the number of available model organisms. For example, there are very few model organisms for diverse microbial eukaryotes, a group of organisms which indeed represents the vast majority of the eukaryotic diversity. To this end, we have developed a system to do genetic modification on the Atlantic salmon pathogen Spironucleus salmonicida. Using this system we could show that the organism is capable of producing hydrogen within specialised compartments. Here we present the genome sequence of S. salmonicida together with a thorough annotation. We compare the results with the closest available model organism, the human intestinal parasite Giardia intestinalis. The fish parasite has a more elaborate system for regulation of gene expression, as well as a larger metabolic capacity. This indicates that S. salmonicida is a well-adapted pathogen that can deal with fluctuating environments, an important trait to be able to establish systemic infections in the host. The development of S. salmonicida into a model system will benefit the studies of fish infections, as well as cell biological processes.

Selenium effects on the metabolism of a Se-metabolizingLactobacillus reuteri: analysis of envelope-enriched and extracellular proteomes

Fixation of selenium byLactobacillus reuteriLb2 BM DSM 16143 into secreted proteins as selenocysteine.

An enzymatic route to selenazolines

Ringing the changes: Selenazolines have applications in medicinal chemistry, but their synthesis is challenging. We report a new convenient and less toxic route to these heterocycles that starts from commercially available selenocysteine. The new route depends on a heterocyclase enzyme that creates oxazolines and thiazolines from serines/threonines and cysteines.

TIGRFAMs and Genome Properties in 2013

TIGRFAMs, available online at http://www.jcvi.org/tigrfams is a database of protein family definitions. Each entry features a seed alignment of trusted representative sequences, a hidden Markov model (HMM) built from that alignment, cutoff scores that let automated annotation pipelines decide which proteins are members, and annotations for transfer onto member proteins. Most TIGRFAMs models are designated equivalog, meaning they assign a specific name to proteins conserved in function from a common ancestral sequence. Models describing more functionally heterogeneous families are designated subfamily or domain, and assign less specific but more widely applicable annotations. The Genome Properties database, available at http://www.jcvi.org/genome-properties, specifies how computed evidence, including TIGRFAMs HMM results, should be used to judge whether an enzymatic pathway, a protein complex or another type of molecular subsystem is encoded in a genome. TIGRFAMs and Genome Properties content are developed in concert because subsystems reconstruction for large numbers of genomes guides selection of seed alignment sequences and cutoff values during protein family construction. Both databases specialize heavily in bacterial and archaeal subsystems. At present, 4284 models appear in TIGRFAMs, while 628 systems are described by Genome Properties. Content derives both from subsystem discovery work and from biocuration of the scientific literature.

ProPhylo: partial phylogenetic profiling to guide protein family construction and assignment of biological process

Background

Phylogenetic profiling is a technique of scoring co-occurrence between a protein family and some other trait, usually another protein family, across a set of taxonomic groups. In spite of several refinements in recent years, the technique still invites significant improvement. To be its most effective, a phylogenetic profiling algorithm must be able to examine co-occurrences among protein families whose boundaries are uncertain within large homologous protein superfamilies.

Results

Partial Phylogenetic Profiling (PPP) is an iterative algorithm that scores a given taxonomic profile against the taxonomic distribution of families for all proteins in a genome. The method works through optimizing the boundary of each protein family, rather than by relying on prebuilt protein families or fixed sequence similarity thresholds. Double Partial Phylogenetic Profiling (DPPP) is a related procedure that begins with a single sequence and searches for optimal granularities for its surrounding protein family in order to generate the best query profiles for PPP. We present ProPhylo, a high-performance software package for phylogenetic profiling studies through creating individually optimized protein family boundaries. ProPhylo provides precomputed databases for immediate use and tools for manipulating the taxonomic profiles used as queries.

Conclusion

ProPhylo results show universal markers of methanogenesis, a new DNA phosphorothioation-dependent restriction enzyme, and efficacy in guiding protein family construction. The software and the associated databases are freely available under the open source Perl Artistic License from ftp://ftp.jcvi.org/pub/data/ppp/.

More than 200 genes required for methane formation from H₂ and CO₂ and energy conservation are present in Methanothermobacter marburgensis and Methanothermobacter thermautotrophicus

The hydrogenotrophic methanogens Methanothermobacter marburgensis and Methanothermobacter thermautotrophicus can easily be mass cultured. They have therefore been used almost exclusively to study the biochemistry of methanogenesis from H₂ and CO₂, and the genomes of these two model organisms have been sequenced. The close relationship of the two organisms is reflected in their genomic architecture and coding potential. Within the 1,607 protein coding sequences (CDS) in common, we identified approximately 200 CDS required for the synthesis of the enzymes, coenzymes, and prosthetic groups involved in CO₂ reduction to methane and in coupling this process with the phosphorylation of ADP. Approximately 20 additional genes, such as those for the biosynthesis of F₄₃₀ and methanofuran and for the posttranslational modifications of the two methyl-coenzyme M reductases, remain to be identified.

An ancient family of SelB elongation factor-like proteins with a broad but disjunct distribution across archaea

Background

SelB is the dedicated elongation factor for delivery of selenocysteinyl-tRNA to the ribosome. In archaea, only a subset of methanogens utilizes selenocysteine and encodes archaeal SelB (aSelB). A SelB-like (aSelBL) homolog has previously been identified in an archaeon that does not encode selenosysteine, and has been proposed to be a pyrrolysyl-tRNA-specific elongation factor (EF-Pyl). However, elongation factor EF-Tu is capable of binding archaeal Pyl-tRNA in bacteria, suggesting the archaeal ortholog EF1A may also be capable of delivering Pyl-tRNA to the ribosome without the need of a specialized factor.

Results

We have phylogenetically characterized the aSelB and aSelBL families in archaea. We find the distribution of aSelBL to be wider than both selenocysteine and pyrrolysine usage. The aSelBLs also lack the carboxy terminal domain usually involved in recognition of the selenocysteine insertion sequence in the target mRNA. While most aSelBL-encoding archaea are methanogenic Euryarchaea, we also find aSelBL representatives in Sulfolobales and Thermoproteales of Crenarchaea, and in the recently identified phylum Thaumarchaea, suggesting that aSelBL evolution has involved horizontal gene transfer and/or parallel loss. Severe disruption of the GTPase domain suggests that some family members may employ a hitherto unknown mechanism of nucleotide hydrolysis, or have lost their GTPase ability altogether. However, patterns of sequence conservation indicate that aSelBL is still capable of binding the ribosome and aminoacyl-tRNA.

Conclusions

Although it is closely related to SelB, aSelBL appears unlikely to either bind selenocysteinyl-tRNA or function as a classical GTP hydrolyzing elongation factor. We propose that following duplication of aSelB, the resultant aSelBL was recruited for binding another aminoacyl-tRNA. In bacteria, aminoacylation with selenocysteine is essential for efficient thermodynamic coupling of SelB binding to tRNA and GTP. Therefore, change in tRNA specificity of aSelBL could have disrupted its GTPase cycle, leading to relaxation of selective pressure on the GTPase domain and explaining its apparent degradation. While the specific role of aSelBL is yet to be experimentally tested, its broad phylogenetic distribution, surpassing that of aSelB, indicates its importance.

A selenium-dependent xanthine dehydrogenase triggers biofilm proliferation in Enterococcus faecalis through oxidant production

Selenium has been shown to be present as a labile cofactor in a small class of molybdenum hydroxylase enzymes in several species of clostridia that specialize in the fermentation of purines and pyrimidines. This labile cofactor is poorly understood, yet recent bioinformatic studies have suggested that Enterococcus faecalis could serve as a model system to better understand the way in which this enzyme cofactor is built and the role of these metalloenzymes in the physiology of the organism. An mRNA that encodes a predicted selenium-dependent molybdenum hydroxylase (SDMH) has also been shown to be specifically increased during the transition from planktonic growth to biofilm growth. Based on these studies, we examined whether this organism produces an SDMH and probed whether selenoproteins may play a role in biofilm physiology. We observed a substantial increase in biofilm density upon the addition of uric acid to cells grown in a defined culture medium, but only when molybdate (Mo) and selenite (Se) were also added. We also observed a significant increase in biofilm density in cells cultured in tryptic soy broth with 1% glucose (TSBG) when selenite was added. In-frame deletion of selD, which encodes selenophosphate synthetase, also blocked biofilm formation that occurred upon addition of selenium. Moreover, mutation in the gene encoding the molybdoenzyme (xdh) prevented the induction of biofilm proliferation upon supplementation with selenium. Tungstate or auranofin addition also blocked this enhanced biofilm density, likely through inhibition of molybdenum or selenium cofactor synthesis. A large protein complex labeled with (75)Se is present in higher concentrations in biofilms than in planktonic cells, and the same complex is formed in TSBG. Xanthine dehydrogenase activity correlates with the presence of this labile selenoprotein complex and is absent in a selD or an xdh mutant. Enhanced biofilm density correlates strongly with higher levels of extracellular peroxide, which is produced upon the addition of selenite to TSBG. Peroxide levels are not increased in either the selD or the xdh mutant upon addition of selenite. Extracellular superoxide production, a phenomenon well established to be linked to clinical isolates, is abolished in both mutant strains. Taken together, these data provide evidence that an SDMH is involved in biofilm formation in Enterococcus faecalis, contributing to oxidant production either directly or alternatively through its involvement in redox-dependent processes linked to oxidant production.

In silico identification of genes involved in selenium metabolism: evidence for a third selenium utilization trait

Background

Selenium (Se) is a trace element that occurs in proteins in the form of selenocysteine (Sec) and in tRNAs in the form of selenouridine (SeU). Selenophosphate synthetase (SelD) is required for both utilization traits. However, previous research also revealed SelDs in two organisms lacking Sec and SeU, suggesting a possible additional use of Se that is dependent on SelD.

Results

In this study, we conducted comparative genomics and phylogenetic analyses to characterize genes involved in Se utilization. Candidate genes identified included SelA/SelB and YbbB that define Sec and SeU pathways, respectively, and NADH oxidoreductase that is predicted to generate a SelD substrate. In addition, among 227 organisms containing SelD, 10 prokaryotes were identified that lacked SelA/SelB and YbbB. Investigation of selD neighboring genes in these organisms revealed a SirA-like protein and two hypothetical proteins HP1 and HP2 that were strongly linked to a novel Se utilization. With these new signature proteins, 32 bacteria and archaea were found that utilized these proteins, likely as part of the new Se utilization trait. Metabolic labeling of one organism containing an orphan SelD, Enterococcus faecalis, with ⁷⁵Se revealed a protein containing labile Se species that could be released by treatment with reducing agents, suggesting non-Sec utilization of Se in this organism.

Conclusion

These studies suggest the occurrence of a third Se utilization trait in bacteria and archaea.