Dynamic evolution of selenocysteine utilization in bacteria: a balance between selenoprotein loss and evolution of selenocysteine from redox active cysteine residues

An RNA modification enzyme directly senses reactive oxygen species for translational regulation in Enterococcus faecalis

Bacteria possess elaborate systems to manage reactive oxygen and nitrogen species (ROS) arising from exposure to the mammalian immune system and environmental stresses. Here we report the discovery of an ROS-sensing RNA-modifying enzyme that regulates translation of stress-response proteins in the gut commensal and opportunistic pathogen Enterococcus faecalis. We analyze the tRNA epitranscriptome of E. faecalis in response to reactive oxygen species (ROS) or sublethal doses of ROS-inducing antibiotics and identify large decreases in N²-methyladenosine (m²A) in both 23 S ribosomal RNA and transfer RNA. This we determine to be due to ROS-mediated inactivation of the Fe-S cluster-containing methyltransferase, RlmN. Genetic knockout of RlmN gives rise to a proteome that mimics the oxidative stress response, with an increase in levels of superoxide dismutase and decrease in virulence proteins. While tRNA modifications were established to be dynamic for fine-tuning translation, here we report the discovery of a dynamically regulated, environmentally responsive rRNA modification. These studies lead to a model in which RlmN serves as a redox-sensitive molecular switch, directly relaying oxidative stress to modulating translation through the rRNA and the tRNA epitranscriptome, adding a different paradigm in which RNA modifications can directly regulate the proteome.

Selenium and Selenoproteins in Health

Selenium is a trace mineral that is essential for health. After being obtained from food and taken up by the liver, selenium performs various physiological functions in the body in the form of selenoproteins, which are best known for their redox activity and anti-inflammatory properties. Selenium stimulates the activation of immune cells and is important for the activation of the immune system. Selenium is also essential for the maintenance of brain function. Selenium supplements can regulate lipid metabolism, cell apoptosis, and autophagy, and have displayed significant alleviating effects in most cardiovascular diseases. However, the effect of increased selenium intake on the risk of cancer remains unclear. Elevated serum selenium levels are associated with an increased risk of type 2 diabetes, and this relationship is complex and nonlinear. Selenium supplementation seems beneficial to some extent; however, existing studies have not fully explained the influence of selenium on various diseases. Further, more intervention trials are needed to verify the beneficial or harmful effects of selenium supplementation in various diseases.

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.

Decoding the Virtual 2D Map of the Chloroplast Proteomes

BACKGROUND

The chloroplast is a semi-autonomous organelle having its own genome and corresponding proteome. Although chloroplast genomes have been reported, no reports exist on their corresponding proteomes. Therefore, a proteome-wide analysis of the chloroplast proteomes of 2893 species was conducted, and a virtual 2D map was constructed.

RESULTS

The resulting virtual 2D map of the chloroplast proteome exhibited a bimodal distribution. The molecular mass of the chloroplast proteome ranged from 0.448 to 616.334 kDa, and the isoelectric point (pI) ranged from 2.854 to 12.954. Chloroplast proteomes were dominated by basic pI proteins with an average pI of 7.852. The molecular weight and isoelectric point of chloroplast proteome were found to show bimodal distribution. Leu was the most abundant and Cys the least abundant amino acid in the chloroplast proteome. Notably, Trp amino acid was absent in the chloroplast protein sequences of Pilostyles aethiopica. In addition, Selenocysteine (Sec) and Pyrrolysine (Pyl) amino acids were also found to be lacking in the chloroplast proteomes.

CONCLUSION

The virtual 2D map and amino acid composition of chloroplast proteome will enable the researchers to understand the biochemistry of chloroplast protein in detail. Further, the amino acid composition of the chloroplast proteome will also allow us to understand the codon usage bias. The codon usage bias and amino acid usage bias of chloroplast will be crucial to understanding their relationship.

SEPHS1: Its evolution, function and roles in development and diseases

Selenophosphate synthetase (SEPHS) was originally discovered in prokaryotes as an enzyme that catalyzes selenophosphate synthesis using inorganic selenium and ATP as substrates. However, in contrast to prokaryotes, two paralogs, SEPHS1 and SEPHS2, occur in many eukaryotes. Prokaryotic SEPHS, also known as SelD, contains either cysteine (Cys) or selenocysteine (Sec) in the catalytic domain. In eukaryotes, only SEPHS2 carries out selenophosphate synthesis and contains Sec at the active site. However, SEPHS1 contains amino acids other than Sec or Cys at the catalytic position. Phylogenetic analysis of SEPHSs reveals that the ancestral SEPHS contains both selenophosphate synthesis and another unknown activity, and that SEPHS1 lost the selenophosphate synthesis activity. The three-dimensional structure of SEPHS1 suggests that its homodimer is unable to form selenophosphate, but retains ATPase activity to produce ADP and inorganic phosphate. The most prominent function of SEPHS1 is that it is implicated in the regulation of cellular redox homeostasis. Deficiency of SEPHS1 leads to the disturbance in the expression of genes involved in redox homeostasis. Different types of reactive oxygen species (ROS) are accumulated in response to SEPHS deficiency depending on cell or tissue types. The accumulation of ROS causes pleiotropic effects such as growth retardation, apoptosis, DNA damage, and embryonic lethality. SEPHS1 deficiency in mouse embryos affects retinoic signaling and other related signaling pathways depending on the embryonal stage until the embryo dies at E11.5. Dysregulated SEPHS1 is associated with the pathogenesis of various diseases including cancer, Crohn's disease, and osteoarthritis.

Thioredoxin reductase selenoproteins from different organisms as potential drug targets for treatment of human diseases

Human thioredoxin reductase (TrxR) is a selenoprotein with a central role in cellular redox homeostasis, utilizing a highly reactive and solvent-exposed selenocysteine (Sec) residue in its active site. Pharmacological modulation of TrxR can be obtained with several classes of small compounds showing different mechanisms of action, but most often dependent upon interactions with its Sec residue. The clinical implications of TrxR modulation as mediated by small compounds have been studied in diverse diseases, from rheumatoid arthritis and ischemia to cancer and parasitic infections. The possible involvement of TrxR in these diseases was in some cases serendipitously discovered, by finding that existing clinically used drugs are also TrxR inhibitors. Inhibiting isoforms of human TrxR is, however, not the only strategy for human disease treatment, as some pathogenic parasites also depend upon Sec-containing TrxR variants, including S. mansoni, B. malayi or O. volvulus. Inhibiting parasite TrxR has been shown to selectively kill parasites and can thus become a promising treatment strategy, especially in the context of quickly emerging resistance towards other drugs. Here we have summarized the basis for the targeting of selenoprotein TrxR variants with small molecules for therapeutic purposes in different human disease contexts. We discuss how Sec engagement appears to be an indispensable part of treatment efficacy and how some therapeutically promising compounds have been evaluated in preclinical or clinical studies. Several research questions remain before a wider application of selenoprotein TrxR inhibition as a first-line treatment strategy might be developed. These include further mechanistic studies of downstream effects that may mediate treatment efficacy, identification of isoform-specific enzyme inhibition patterns for some given therapeutic compounds, and the further elucidation of cell-specific effects in disease contexts such as in the tumor microenvironment or in host-parasite interactions, and which of these effects may be dependent upon the specific targeting of Sec in distinct TrxR isoforms.

Functions and Regulation of Translation Elongation Factors

Translation elongation is a key step of protein synthesis, during which the nascent polypeptide chain extends by one amino acid residue during one elongation cycle. More and more data revealed that the elongation is a key regulatory node for translational control in health and disease. During elongation, elongation factor Tu (EF-Tu, eEF1A in eukaryotes) is used to deliver aminoacyl-tRNA (aa-tRNA) to the A-site of the ribosome, and elongation factor G (EF-G, EF2 in eukaryotes and archaea) is used to facilitate the translocation of the tRNA₂-mRNA complex on the ribosome. Other elongation factors, such as EF-Ts/eEF1B, EF-P/eIF5A, EF4, eEF3, SelB/EFsec, TetO/Tet(M), RelA and BipA, have been found to affect the overall rate of elongation. Here, we made a systematic review on the canonical and non-canonical functions and regulation of these elongation factors. In particular, we discussed the close link between translational factors and human diseases, and clarified how post-translational modifications control the activity of translational factors in tumors.

Fungal genomes: suffering with functional annotation errors

Background

The genome sequence data of more than 65985 species are publicly available as of October 2021 within the National Center for Biotechnology Information (NCBI) database alone and additional genome sequences are available in other databases and also continue to accumulate at a rapid pace. However, an error-free functional annotation of these genome is essential for the research communities to fully utilize these data in an optimum and efficient manner.

Results

An analysis of proteome sequence data of 689 fungal species (7.15 million protein sequences) was conducted to identify the presence of functional annotation errors. Proteins associated with calcium signaling events, including calcium dependent protein kinases (CDPKs), calmodulins (CaM), calmodulin-like (CML) proteins, WRKY transcription factors, selenoproteins, and proteins associated with the terpene biosynthesis pathway, were targeted in the analysis. Gene associated with CDPKs and selenoproteins are known to be absent in fungal genomes. Our analysis, however, revealed the presence of proteins that were functionally annotated as CDPK proteins. However, InterproScan analysis indicated that none of the protein sequences annotated as “calcium dependent protein kinase” were found to encode calcium binding EF-hands at the regulatory domain. Similarly, none of a protein sequences annotated as a “selenocysteine” were found to contain a Sec (U) amino acid. Proteins annotated as CaM and CMLs also had significant discrepancies. CaM proteins should contain four calcium binding EF-hands, however, a range of 2–4 calcium binding EF-hands were present in the fungal proteins that were annotated as CaM proteins. Similarly, CMLs should possess four calcium binding EF-hands, but some of the CML annotated fungal proteins possessed either three or four calcium binding EF-hands. WRKY transcription factors are characterized by the presence of a WRKY domain and are confined to the plant kingdom. Several fungal proteins, however, were annotated as WRKY transcription factors, even though they did not contain a WRKY domain.

Conclusion

The presence of functional annotation errors in fungal genome and proteome databases is of considerable concern and needs to be addressed in a timely manner.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43008-021-00083-x.

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.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

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.

Programmed Deviations of Ribosomes From Standard Decoding in Archaea

Genetic code decoding, initially considered to be universal and immutable, is now known to be flexible. In fact, in specific genes, ribosomes deviate from the standard translational rules in a programmed way, a phenomenon globally termed recoding. Translational recoding, which has been found in all domains of life, includes a group of events occurring during gene translation, namely stop codon readthrough, programmed ± 1 frameshifting, and ribosome bypassing. These events regulate protein expression at translational level and their mechanisms are well known and characterized in viruses, bacteria and eukaryotes. In this review we summarize the current state-of-the-art of recoding in the third domain of life. In Archaea, it was demonstrated and extensively studied that translational recoding regulates the decoding of the 21st and the 22nd amino acids selenocysteine and pyrrolysine, respectively, and only one case of programmed -1 frameshifting has been reported so far in Saccharolobus solfataricus P2. However, further putative events of translational recoding have been hypothesized in other archaeal species, but not extensively studied and confirmed yet. Although this phenomenon could have some implication for the physiology and adaptation of life in extreme environments, this field is still underexplored and genes whose expression could be regulated by recoding are still poorly characterized. The study of these recoding episodes in Archaea is urgently needed.

coupleCoC+: An information-theoretic co-clustering-based transfer learning framework for the integrative analysis of single-cell genomic data

Technological advances have enabled us to profile multiple molecular layers at unprecedented single-cell resolution and the available datasets from multiple samples or domains are growing. These datasets, including scRNA-seq data, scATAC-seq data and sc-methylation data, usually have different powers in identifying the unknown cell types through clustering. So, methods that integrate multiple datasets can potentially lead to a better clustering performance. Here we propose coupleCoC+ for the integrative analysis of single-cell genomic data. coupleCoC+ is a transfer learning method based on the information-theoretic co-clustering framework. In coupleCoC+, we utilize the information in one dataset, the source data, to facilitate the analysis of another dataset, the target data. coupleCoC+ uses the linked features in the two datasets for effective knowledge transfer, and it also uses the information of the features in the target data that are unlinked with the source data. In addition, coupleCoC+ matches similar cell types across the source data and the target data. By applying coupleCoC+ to the integrative clustering of mouse cortex scATAC-seq data and scRNA-seq data, mouse and human scRNA-seq data, mouse cortex sc-methylation and scRNA-seq data, and human blood dendritic cells scRNA-seq data from two batches, we demonstrate that coupleCoC+ improves the overall clustering performance and matches the cell subpopulations across multimodal single-cell genomic datasets. coupleCoC+ has fast convergence and it is computationally efficient. The software is available at https://github.com/cuhklinlab/coupleCoC_plus.

The recent advances in single-cell technologies have enabled multiple biological layers to be probed and provides unprecedented opportunities to assay cellular heterogeneity. To analyze the complex biological processes varying across cells, we need to obtain and integrate different types of genomic features through flexible but rigorous computational methods. The most important challenge for data integration is to link data from different sources in a way that is biologically meaningful. In this work, we have developed a transfer learning method based on the information-theoretic co-clustering framework for the integrative analysis of single-cell genomic data. This method utilizes the information from one dataset to boost the analysis of another dataset, and it also uses the information of the features that are unlinked in the two datasets. We demonstrate that our transfer learning-based clustering method significantly improves clustering performance in single-cell genomic datasets. Our results show that transfer learning is promising for the integrative analysis of single-cell genomic data.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Pan-Genome Analysis Reveals Host-Specific Functional Divergences in Burkholderia gladioli

Burkholderia gladioli has high versatility and adaptability to various ecological niches. Here, we constructed a pan-genome using 14 genome sequences of B. gladioli, which originate from different niches, including gladiolus, rice, humans, and nature. Functional roles of core and niche-associated genomes were investigated by pathway enrichment analyses. Consequently, we inferred the uniquely important role of niche-associated genomes in (1) selenium availability during competition with gladiolus host; (2) aromatic compound degradation in seed-borne and crude oil-accumulated environments, and (3) stress-induced DNA repair system/recombination in the cystic fibrosis-niche. We also identified the conservation of the rhizomide biosynthetic gene cluster in all the B. gladioli strains and the concentrated distribution of this cluster in human isolates. It was confirmed the absence of complete CRISPR/Cas system in both plant and human pathogenic B. gladioli and the presence of the system in B. gladioli living in nature, possibly reflecting the inverse relationship between CRISPR/Cas system and virulence.

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.

Bioinformatic Prediction of an tRNASec Gene Nested inside an Elongation Factor SelB Gene in Alphaproteobacteria

In bacteria, selenocysteine (Sec) is incorporated into proteins via the recoding of a particular codon, the UGA stop codon in most cases. Sec-tRNASec is delivered to the ribosome by the Sec-dedicated elongation factor SelB that also recognizes a Sec-insertion sequence element following the codon on the mRNA. Since the excess of SelB may lead to sequestration of Sec-tRNASec under selenium deficiency or oxidative stress, the expression levels of SelB and tRNASec should be regulated. In this bioinformatic study, I analyzed the Rhizobiales SelB species because they were annotated to have a non-canonical C-terminal extension. I found that the open reading frame (ORF) of diverse Alphaproteobacteria selB genes includes an entire tRNASec sequence (selC) and overlaps with the start codon of the downstream ORF. A remnant tRNASec sequence was found in the Sinorhizobium melilotiselB genes whose products have a shorter C-terminal extension. Similar overlapping traits were found in Gammaproteobacteria and Nitrospirae. I hypothesized that once the tRNASec moiety is folded and processed, the expression of the full-length SelB may be repressed. This is the first report on a nested tRNA gene inside a protein ORF in bacteria.

Selective cysteine-to-selenocysteine changes in a [NiFe]-hydrogenase confirm a special position for catalysis and oxygen tolerance

In [NiFe]-hydrogenases, the active-site Ni is coordinated by four cysteine-S ligands (Cys; C), two of which are bridging to the Fe(CO)(CN)2 fragment. Substitution of a single Cys residue by selenocysteine (Sec; U) occurs occasionally in nature. Using a recent method for site-specific Sec incorporation into proteins, each of the four Ni-coordinating cysteine residues in the oxygen-tolerant Escherichia coli [NiFe]-hydrogenase-1 (Hyd-1) has been replaced by U to identify its importance for enzyme function. Steady-state solution activity of each Sec-substituted enzyme (on a per-milligram basis) is lowered, although this may reflect the unquantified presence of recalcitrant inactive/immature/misfolded forms. Protein film electrochemistry, however, reveals detailed kinetic data that are independent of absolute activities. Like native Hyd-1, the variants have low apparent KMH2 values, do not produce H2 at pH 6, and display the same onset overpotential for H2 oxidation. Mechanistically important differences were identified for the C576U variant bearing the equivalent replacement found in native [NiFeSe]-hydrogenases, its extreme O2 tolerance (apparent KMH2 and Vmax [solution] values relative to native Hyd-1 of 0.13 and 0.04, respectively) implying the importance of a selenium atom in the position cis to the site where exogenous ligands (H-, H2, O2) bind. Observation of the same unusual electrocatalytic signature seen earlier for the proton transfer-defective E28Q variant highlights the direct role of the chalcogen atom (S/Se) at position 576 close to E28, with the caveat that Se is less effective than S in facilitating proton transfer away from the Ni during H2 oxidation by this enzyme.

Biogeographic and Evolutionary Patterns of Trace Element Utilization in Marine Microbial World

Trace elements are required by all organisms, which are key components of many enzymes catalyzing important biological reactions. Many trace element-dependent proteins have been characterized; however, little is known about their occurrence in microbial communities in diverse environments, especially the global marine ecosystem. Moreover, the relationships between trace element utilization and different types of environmental stressors are unclear. In this study, we used metagenomic data from the Global Ocean Sampling expedition project to identify the biogeographic distribution of genes encoding trace element-dependent proteins (for copper, molybdenum, cobalt, nickel, and selenium) in a variety of marine and non-marine aquatic samples. More than 56,000 metalloprotein and selenoprotein genes corresponding to nearly 100 families were predicted, becoming the largest dataset of marine metalloprotein and selenoprotein genes reported to date. In addition, samples with enriched or depleted metalloprotein/selenoprotein genes were identified, suggesting an active or inactive usage of these micronutrients in various sites. Further analysis of interactions among the elements showed significant correlations between some of them, especially those between nickel and selenium/copper. Finally, investigation of the relationships between environmental conditions and metalloprotein/selenoprotein families revealed that many environmental factors might contribute to the evolution of different metalloprotein and/or selenoprotein genes in the marine microbial world. Our data provide new insights into the utilization and biological roles of these trace elements in extant marine microbes, and might also be helpful for the understanding of how these organisms have adapted to their local environments.

Construction of anti-codon table of the plant kingdom and evolution of tRNA selenocysteine (tRNASec)

Background

The tRNAs act as a bridge between the coding mRNA and incoming amino acids during protein translation. The anti-codon of tRNA recognizes the codon of the mRNA and deliver the amino acid into the protein translation chain. However, we did not know about the exact abundance of anti-codons in the genome and whether the frequency of abundance remains same across the plant lineage or not.

Results

Therefore, we analysed the tRNAnome of 128 plant species and reported an anti-codon table of the plant kingdom. We found that CAU anti-codon of tRNA^Met has highest (5.039%) whereas GCG anti-codon of tRNA^Arg has lowest (0.004%) abundance. However, when we compared the anti-codon frequencies according to the tRNA isotypes, we found tRNA^Leu (7.808%) has highest abundance followed by tRNA^Ser (7.668%) and tRNA^Gly (7.523%). Similarly, suppressor tRNA (0.036%) has lowest abundance followed by tRNA^Sec (0.066%) and tRNA^His (2.109). The genome of Ipomoea nil, Papaver somniferum, and Zea mays encoded the highest number of anti-codons (isoacceptor) at 59 each whereas the genome of Ostreococcus tauri was found to encode only 18 isoacceptors. The tRNA^Sec genes undergone losses more frequently than duplication and we found that tRNA^Sec showed anti-codon switch during the course of evolution.

Conclusion

The anti-codon table of the plant tRNA will enable us to understand the synonymous codon usage of the plant kingdom and can be very helpful to understand which codon is preferred over other during the translation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07216-3.

Microbial selenium metabolism: a brief history, biogeochemistry and ecophysiology

Selenium is an essential trace element for organisms from all three domains of life. Microorganisms, in particular, mediate reductive transformations of selenium that govern the element's mobility and bioavailability in terrestrial and aquatic environments. Selenium metabolism is not just ubiquitous but an ancient feature of life likely extending back to the universal common ancestor of all cellular lineages. As with the sulfur biogeochemical cycle, reductive transformations of selenium serve two metabolic functions: assimilation into macromolecules and dissimilatory reduction during anaerobic respiration. This review begins with a historical overview of how research in both aspects of selenium metabolism has developed. We then provide an overview of the global selenium biogeochemical cycle, emphasizing the central role of microorganisms in the cycle. This serves as a basis for a robust discussion of current models for the evolution of the selenium biogeochemical cycle over geologic time, and how knowledge of the evolution and ecophysiology of selenium metabolism can enrich and refine these models. We conclude with a discussion of the ecophysiological function of selenium-respiring prokaryotes within the cycle, and the tantalizing possibility of oxidative selenium transformations during chemolithoautotrophic growth.

Bioinformatics of Selenoproteins

Significance: Bioinformatics has brought important insights into the field of selenium research. The progress made in the development of computational tools in the last two decades, coordinated with growing genome resources, provided new opportunities to study selenoproteins. The present review discusses existing tools for selenoprotein gene finding and other bioinformatic approaches to study the biology of selenium. Recent Advances: The availability of complete selenoproteomes allowed assessing a global distribution of the use of selenocysteine (Sec) across the tree of life, as well as studying the evolution of selenoproteins and their biosynthetic pathway. Beyond gene identification and characterization, human genetic variants in selenoprotein genes were used to examine adaptations to selenium levels in diverse human populations and to estimate selective constraints against gene loss. Critical Issues: The synthesis of selenoproteins is essential for development in mice. In humans, several mutations in selenoprotein genes have been linked to rare congenital disorders. And yet, the mechanism of Sec insertion and the regulation of selenoprotein synthesis in mammalian cells are not completely understood. Future Directions: Omics technologies offer new possibilities to study selenoproteins and mechanisms of Sec incorporation in cells, tissues, and organisms.

Bioinformatics of Metalloproteins and Metalloproteomes

Trace metals are inorganic elements that are required for all organisms in very low quantities. They serve as cofactors and activators of metalloproteins involved in a variety of key cellular processes. While substantial effort has been made in experimental characterization of metalloproteins and their functions, the application of bioinformatics in the research of metalloproteins and metalloproteomes is still limited. In the last few years, computational prediction and comparative genomics of metalloprotein genes have arisen, which provide significant insights into their distribution, function, and evolution in nature. This review aims to offer an overview of recent advances in bioinformatic analysis of metalloproteins, mainly focusing on metalloprotein prediction and the use of different metals across the tree of life. We describe current computational approaches for the identification of metalloprotein genes and metal-binding sites/patterns in proteins, and then introduce a set of related databases. Furthermore, we discuss the latest research progress in comparative genomics of several important metals in both prokaryotes and eukaryotes, which demonstrates divergent and dynamic evolutionary patterns of different metalloprotein families and metalloproteomes. Overall, bioinformatic studies of metalloproteins provide a foundation for systematic understanding of trace metal utilization in all three domains of life.

Tetrahymena Glutathione Peroxidase Family: A Comparative Analysis of These Antioxidant Enzymes and Differential Gene Expression to Metals and Oxidizing Agents

In the present work, an extensive analysis of the putative glutathione peroxidases (GPx) of the eukaryotic microorganism model Tetrahymena thermophila is carried out. A comparative analysis with GPx present in other Tetrahymena species and other very taxonomically diverse ciliates is also performed. A majority of ciliate GPx have replaced the selenocysteine (Sec) by Cys in its catalytic center, so they can be considered as phospholipid hydroperoxide glutathione peroxidases (PHGPx). Selenocysteine insertion sequence (SECIS) elements have been detected in several ciliate GPx that do not incorporate Sec in their amino acid sequences, and conversely, in other ciliate GPx with Sec, no SECIS elements are detected. These anomalies are analyzed and discussed. From the phylogenetic analysis using the ciliate GPx amino acid sequences, the existence of extensive intra- and interspecific gene duplications that produced multiple GPx isoforms in each species is inferred. The ancestral character of the selenoproteins is also corroborated. The analysis by qRT-PCR of six selected T. thermophila GPx genes has shown a quantitative differential expression between them, depending on the stressor (oxidizing agents, apoptotic inducer or metals) and the time of exposure.

Methane, arsenic, selenium and the origins of the DMSO reductase family

Mononuclear molybdoenzymes of the dimethyl sulfoxide reductase (DMSOR) family catalyze a number of reactions essential to the carbon, nitrogen, sulfur, arsenic, and selenium biogeochemical cycles. These enzymes are also ancient, with many lineages likely predating the divergence of the last universal common ancestor into the Bacteria and Archaea domains. We have constructed rooted phylogenies for over 1,550 representatives of the DMSOR family using maximum likelihood methods to investigate the evolution of the arsenic biogeochemical cycle. The phylogenetic analysis provides compelling evidence that formylmethanofuran dehydrogenase B subunits, which catalyze the reduction of CO₂ to formate during hydrogenotrophic methanogenesis, constitutes the most ancient lineage. Our analysis also provides robust support for selenocysteine as the ancestral ligand for the Mo/W atom. Finally, we demonstrate that anaerobic arsenite oxidase and respiratory arsenate reductase catalytic subunits represent a more ancient lineage of DMSORs compared to aerobic arsenite oxidase catalytic subunits, which evolved from the assimilatory nitrate reductase lineage. This provides substantial support for an active arsenic biogeochemical cycle on the anoxic Archean Earth. Our work emphasizes that the use of chalcophilic elements as substrates as well as the Mo/W ligand in DMSORs has indelibly shaped the diversification of these enzymes through deep time.

Studying O2 pathways in [NiFe]- and [NiFeSe]-hydrogenases

Hydrogenases are efficient biocatalysts for H₂ production and oxidation with various potential biotechnological applications.[NiFe]-class hydrogenases are highly active in both production and oxidation processes—albeit primarily biased to the latter—but suffer from being sensitive to O₂.[NiFeSe] hydrogenases are a subclass of [NiFe] hydrogenases with, usually, an increased insensitivity to aerobic environments. In this study we aim to understand the structural causes of the low sensitivity of a [NiFeSe]-hydrogenase, when compared with a [NiFe] class enzyme, by studying the diffusion of O₂. To unravel the differences between the two enzymes, we used computational methods comprising Molecular Dynamics simulations with explicit O₂ and Implicit Ligand Sampling methodologies. With the latter, we were able to map the free energy landscapes for O₂ permeation in both enzymes. We derived pathways from these energy landscapes and selected the kinetically more relevant ones with reactive flux analysis using transition path theory. These studies evidence the existence of quite different pathways in both enzymes and predict a lower permeation efficiency for O₂ in the case of the [NiFeSe]-hydrogenase when compared with the [NiFe] enzyme. These differences can explain the experimentally observed lower inhibition by O₂ on [NiFeSe]-hydrogenases, when compared with [NiFe]-hydrogenases. A comprehensive map of the residues lining the most important O₂ pathways in both enzymes is also presented.

Metabolic flux and functional potential of microbial community in an acidogenic dynamic membrane bioreactor

In this study, a laboratory-scale acidogenic dynamic membrane bioreactor (DMBR) was operated with or without alginate immobilization under various hydraulic retention times. The interactions between biofilm formation, metabolic flux and predicted functional genes were investigated in the DMBR. Regardless of the alginate immobilization, hydrogen yield was increased after biofilm was formed on a supporting material. Metabolic flux shifts from homoacetogenic to hydrogenic pathways along with biofilm formation was confirmed by flux balance analysis for metabolites. KEGG ortholog count prediction of functional genes from microbial consortia using bioinformatics package (PICRUSt). Predicted enzymes showed similar levels of functional gene expression regardless of the immobilization status. PICRUSt result showed strong evidence of the functional relationship among hydrogen production, biofilm formation, and metabolic pathway. Taken together, this study provides to estimate characteristics of biofilm formation changes and to better understand the microbial metabolism of DMBR.

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.

Role of Selenoproteins in Bacterial Pathogenesis

The trace element selenium is an essential micronutrient that plays an important role in maintaining homeostasis of several tissues including the immune system of mammals. The vast majority of the biological functions of selenium are mediated via selenoproteins, proteins which incorporate the selenium-containing amino acid selenocysteine. Several bacterial infections of humans and animals are associated with decreased levels of selenium in the blood and an adjunct therapy with selenium often leads to favorable outcomes. Many pathogenic bacteria are also capable of synthesizing selenocysteine suggesting that selenoproteins may have a role in bacterial physiology. Interestingly, the composition of host microbiota is also regulated by dietary selenium levels. Therefore, bacterial pathogens, microbiome, and host immune cells may be competing for a limited supply of selenium. Elucidating how selenium, in particular selenoproteins, may regulate pathogen virulence, microbiome diversity, and host immune response during a bacterial infection is critical for clinical management of infectious diseases.

Selenocysteine β-Lyase: Biochemistry, Regulation and Physiological Role of the Selenocysteine Decomposition Enzyme

The enzyme selenocysteine β-lyase (SCLY) was first isolated in 1982 from pig livers, followed by its identification in bacteria. SCLY works as a homodimer, utilizing pyridoxal 5'-phosphate as a cofactor, and catalyzing the specific decomposition of the amino acid selenocysteine into alanine and selenide. The enzyme is thought to deliver its selenide as a substrate for selenophosphate synthetases, which will ultimately be reutilized in selenoprotein synthesis. SCLY subcellular localization is unresolved, as it has been observed both in the cytosol and in the nucleus depending on the technical approach used. The highest SCLY expression and activity in mammals is found in the liver and kidneys. Disruption of the Scly gene in mice led to obesity, hyperinsulinemia, glucose intolerance, and hepatic steatosis, with SCLY being suggested as a participant in the regulation of energy metabolism in a sex-dependent manner. With the physiological role of SCLY still not fully understood, this review attempts to discuss the available literature regarding SCLY in animals and provides avenues for possible future investigation.

RefSeq: an update on prokaryotic genome annotation and curation

The Reference Sequence (RefSeq) project at the National Center for Biotechnology Information (NCBI) provides annotation for over 95 000 prokaryotic genomes that meet standards for sequence quality, completeness, and freedom from contamination. Genomes are annotated by a single Prokaryotic Genome Annotation Pipeline (PGAP) to provide users with a resource that is as consistent and accurate as possible. Notable recent changes include the development of a hierarchical evidence scheme, a new focus on curating annotation evidence sources, the addition and curation of protein profile hidden Markov models (HMMs), release of an updated pipeline (PGAP-4), and comprehensive re-annotation of RefSeq prokaryotic genomes. Antimicrobial resistance proteins have been reannotated comprehensively, improved structural annotation of insertion sequence transposases and selenoproteins is provided, curated complex domain architectures have given upgraded names to millions of multidomain proteins, and we introduce a new kind of annotation rule—BlastRules. Continual curation of supporting evidence, and propagation of improved names onto RefSeq proteins ensures that the functional annotation of genomes is kept current. An increasing share of our annotation now derives from HMMs and other sets of annotation rules that are portable by nature, and available for download and for reuse by other investigators. RefSeq is found at https://www.ncbi.nlm.nih.gov/refseq/.

The molecular mass and isoelectric point of plant proteomes

Background

Cell contain diverse array of proteins with different molecular weight and isoelectric point (pI). The molecular weight and pI of protein play important role in determining the molecular biochemical function. Therefore, it was important to understand the detail regarding the molecular weight and pI of the plant proteins.

Results

A proteome-wide analysis of plant proteomes from 145 species revealed a pI range of 1.99 (epsin) to 13.96 (hypothetical protein). The spectrum of molecular mass of the plant proteins varied from 0.54 to 2236.8 kDa. A putative Type-I polyketide synthase (22244 amino acids) in Volvox carteri was found to be the largest protein in the plant kingdom. However, Type-I polyketide synthase was not found in higher plant species. Titin (806.46 kDa) and misin/midasin (730.02 kDa) were the largest proteins identified in higher plant species. The pI and molecular weight of the plant proteins showed a trimodal distribution. An acidic pI (56.44% of proteins) was found to be predominant over a basic pI (43.34% of proteins) and the abundance of acidic pI proteins was higher in unicellular algae species relative to multicellular higher plants. In contrast, the seaweed, Porphyra umbilicalis, possesses a higher proportion of basic pI proteins (70.09%). Plant proteomes were also found to contain selenocysteine (Sec), amino acid that was found only in lower eukaryotic aquatic plant lineage. Amino acid composition analysis showed Leu was high and Trp was low abundant amino acids in the plant proteome. Additionally, the plant proteomes also possess ambiguous amino acids Xaa (unknown), Asx (asparagine or aspartic acid), Glx (glutamine or glutamic acid), and Xle (leucine or isoleucine) as well.

Conclusion

The diverse molecular weight and isoelectric point range of plant proteome will be helpful to understand their biochemical and functional aspects. The presence of selenocysteine proteins in lower eukaryotic organism is of interest and their expression in higher plant system can help us to understand their functional role.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5983-8) contains supplementary material, which is available to authorized users.

Phylogenomics Provides New Insights into Gains and Losses of Selenoproteins among Archaeplastida

Selenoproteins that contain selenocysteine (Sec) are found in all kingdoms of life. Although they constitute a small proportion of the proteome, selenoproteins play essential roles in many organisms. In photosynthetic eukaryotes, selenoproteins have been found in algae but are missing in land plants (embryophytes). In this study, we explored the evolutionary dynamics of Sec incorporation by conveying a genomic search for the Sec machinery and selenoproteins across Archaeplastida. We identified a complete Sec machinery and variable sizes of selenoproteomes in the main algal lineages. However, the entire Sec machinery was missing in the Bangiophyceae-Florideophyceae clade (BV) of Rhodoplantae (red algae) and only partial machinery was found in three species of Archaeplastida, indicating parallel loss of Sec incorporation in different groups of algae. Further analysis of genome and transcriptome data suggests that all major lineages of streptophyte algae display a complete Sec machinery, although the number of selenoproteins is low in this group, especially in subaerial taxa. We conclude that selenoproteins tend to be lost in Archaeplastida upon adaptation to a subaerial or acidic environment. The high number of redox-active selenoproteins found in some bloom-forming marine microalgae may be related to defense against viral infections. Some of the selenoproteins in these organisms may have been gained by horizontal gene transfer from bacteria.

Prioritized brain selenium retention and selenoprotein expression: Nutritional insights into Parkinson's disease

Selenium (Se), an essential trace mineral, confers its physiological functions mainly through selenoproteins, most of which are oxidoreductases. Results from animal, epidemiological, and human genetic studies link Parkinson's disease to Se and certain selenoproteins. Parkinson's disease is characterized by multiple motor and non-motor symptoms that are difficult to diagnose at early stages of the pathogenesis. While irreversible, degenerative and age-related, the onset of Parkinson's disease may be delayed through proper dietary and environmental controls. One particular attribute of Se biology is that brain has the highest priority to receive and retain this nutrient even in Se deficiency. Thus, brain Se deficiency is rare; however, a strong body of recent evidence implicates selenoprotein dysfunction in Parkinson's disease. Direct and indirect evidence from mouse models implicate selenoprotein T, glutathione peroxidase 1, selenoprotein P and glutathione peroxidase 4 in counteracting Parkinson's disease through Se transportation to the brain and reduced oxidative stress. It is of future interest to further characterize the full selenoproteomes in various types of brain cells and elucidate the mechanism of their actions in Parkinson's disease.

The Thioredoxin-Like Family of Selenoproteins: Implications in Aging and Age-Related Degeneration

The thioredoxin-like (Rdx) family proteins contain four selenoproteins (selenoprotein H, SELENOH; selenoprotein T, SELENOT; selenoprotein V, SELENOV; selenoprotein W, SELENOW) and a nonselenoprotein Rdx12. They share a CxxU or a CxxC (C, cysteine; x, any amino acid; U, selenocysteine) motif and a stretch of eGxFEI(V) sequence. From the evolutionary perspective, SELENOW and SELENOV are clustered together and SELENOH and SELENOT are in another branch. Selenoproteins in the Rdx family exhibit tissue- and organelle-specific distribution and are differentially influenced in response to selenium deficiency. While SELENOH is nucleus-exclusive, SELENOT resides mainly in endoplasmic reticulum and SELENOW in cytosol. SELENOV is expressed essentially only in the testes with unknown cellular localization. SELENOH and SELENOW are more sensitive than SELENOT and SELENOV to selenium deficiency. While physiological functions of the Rdx family of selenoproteins are not fully understand, results from animal models demonstrated that (1) brain-specific SELENOT knockout mice are susceptible to 1-methyl-4-phenylpyridinium-induced Parkinson's disease in association with redox imbalance and (2) adult zebrafishes with heterozygous SELENOH knockout are prone to dimethylbenzanthracene-induced tumorigenesis together with increased DNA damage and oxidative stress. Further animal and human studies are needed to fully understand physiological roles of the Rdx family of selenoproteins in redox regulation, genome maintenance, aging, and age-related degeneration.

Utilization of selenocysteine in early-branching fungal phyla

Selenoproteins are a diverse group of proteins containing selenocysteine (Sec)-the twenty-first amino acid-incorporated during translation via a unique recoding mechanism^1,2. Selenoproteins fulfil essential roles in many organisms¹, yet are not ubiquitous across the tree of life^3-7. In particular, fungi were deemed devoid of selenoproteins^4,5,8. However, we show here that Sec is utilized by nine species belonging to diverse early-branching fungal phyla, as evidenced by the genomic presence of both Sec machinery and selenoproteins. Most fungal selenoproteins lack consensus Sec recoding signals (SECIS elements⁹) but exhibit other RNA structures, suggesting altered mechanisms of Sec insertion in fungi. Phylogenetic analyses support a scenario of vertical inheritance of the Sec trait within eukaryotes and fungi. Sec was then lost in numerous independent events in various fungal lineages. Notably, Sec was lost at the base of Dikarya, resulting in the absence of selenoproteins in Saccharomyces cerevisiae and other well-studied fungi. Our results indicate that, despite scattered occurrence, selenoproteins are found in all kingdoms of life.

A Quantitative Chemoproteomic Platform to Monitor Selenocysteine Reactivity within a Complex Proteome

Mammalian selenocysteine (Sec)-containing proteins, selenoproteins, are important to (patho)physiological processes, including redox homeostasis. Sec residues have been recalcitrant to mass spectrometry-based chemoproteomic methods that enrich for reactive cysteine (Cys) residues with electrophilic chemical probes, despite confirmed reactivity of Sec with these electrophiles. Highly abundant Cys peptides likely suppress low-abundant Sec peptides. By exploiting the decreased pK_a of Sec relative to Cys, we have developed a chemoproteomic platform that relies on low pH (pH 5.75) electrophile labeling, reducing Cys reactivity and enhancing identification of Sec-containing peptides across mouse tissues and cell lines. The utility of this Sec-profiling platform is underscored by evaluation of the selectivity of auranofin, an inhibitor of the selenoprotein, thioredoxin reductase, against both reactive Cys- and Sec-containing proteins. Platform limitations pertain to the non-physiological low-pH conditions that could perturb protein structure and function. Future work necessitates the discovery of Sec-selective electrophiles that function at physiological pH.

Recoding of the selenocysteine UGA codon by cysteine in the presence of a non-canonical tRNACys and elongation factor SelB

In many organisms, the UGA stop codon is recoded to insert selenocysteine (Sec) into proteins. Sec incorporation in bacteria is directed by an mRNA element, known as the Sec-insertion sequence (SECIS), located downstream of the Sec codon. Unlike other aminoacyl-tRNAs, Sec-tRNASec is delivered to the ribosome by a dedicated elongation factor, SelB. We recently identified a series of tRNASec-like tRNA genes distributed across Bacteria that also encode a canonical tRNASec. These tRNAs contain sequence elements generally recognized by cysteinyl-tRNA synthetase (CysRS). While some of these tRNAs contain a UCA Sec anticodon, most have a GCA Cys anticodon. tRNASec with GCA anticodons are known to recode UGA codons. Here we investigate the clostridial Desulfotomaculum nigrificans tRNASec-like tRNACys, and show that this tRNA is acylated by CysRS, recognized by SelB, and capable of UGA recoding with Cys in Escherichia coli. We named this non-canonical group of tRNACys as 'tRNAReC' (Recoding with Cys). We performed a comprehensive survey of tRNAReC genes to establish their phylogenetic distribution, and found that, in a particular lineage of clostridial Pelotomaculum, the Cys identity elements of tRNAReC had mutated. This novel tRNA, which contains a UCA anticodon, is capable of Sec incorporation in E. coli, albeit with lower efficiency relative to Pelotomaculum tRNASec. We renamed this unusual tRNASec derived from tRNAReC as 'tRNAReU' (Recoding with Sec). Together, our results suggest that tRNAReC and tRNAReU may serve as safeguards in the production of selenoproteins and - to our knowledge - they provide the first example of programmed codon-anticodon mispairing in bacteria.

Computational identification of the selenocysteine tRNA (tRNASec) in genomes

Selenocysteine (Sec) is known as the 21st amino acid, a cysteine analogue with selenium replacing sulphur. Sec is inserted co-translationally in a small fraction of proteins called selenoproteins. In selenoprotein genes, the Sec specific tRNA (tRNA^Sec) drives the recoding of highly specific UGA codons from stop signals to Sec. Although found in organisms from the three domains of life, Sec is not universal. Many species are completely devoid of selenoprotein genes and lack the ability to synthesize Sec. Since tRNA^Sec is a key component in selenoprotein biosynthesis, its efficient identification in genomes is instrumental to characterize the utilization of Sec across lineages. Available tRNA prediction methods fail to accurately predict tRNA^Sec, due to its unusual structural fold. Here, we present Secmarker, a method based on manually curated covariance models capturing the specific tRNA^Sec structure in archaea, bacteria and eukaryotes. We exploited the non-universality of Sec to build a proper benchmark set for tRNA^Sec predictions, which is not possible for the predictions of other tRNAs. We show that Secmarker greatly improves the accuracy of previously existing methods constituting a valuable tool to identify tRNA^Sec genes, and to efficiently determine whether a genome contains selenoproteins. We used Secmarker to analyze a large set of fully sequenced genomes, and the results revealed new insights in the biology of tRNA^Sec, led to the discovery of a novel bacterial selenoprotein family, and shed additional light on the phylogenetic distribution of selenoprotein containing genomes. Secmarker is freely accessible for download, or online analysis through a web server at http://secmarker.crg.cat.

Most proteins are made of twenty amino acids. However, there is a small group of proteins that incorporate a 21st amino acid, Selenocysteine (Sec). These proteins are called selenoproteins and are present in some, but not all, species from the three domains of life. Sec is inserted in selenoproteins in response to the UGA codon, normally a stop codon. A Sec specific tRNA (tRNA^Sec), which only exists in the organisms that synthesize selenoproteins recognizes the UGA codon. tRNA^Sec is not only indispensable for Sec incorporation into selenoproteins, but also for Sec synthesis, since Sec is synthesized on its own tRNA. The structure of tRNA^Sec differs from that of canonical tRNAs, and general tRNA detection methods fail to accurately predict it. We developed Secmarker, a tRNA^Sec specific identification tool based on the characteristic structural features of the tRNA^Sec. Our benchmark shows that Secmarker produces nearly flawless tRNA^Sec predictions. We used Secmarker to scan all currently available genome sequences. The analysis of the highly accurate predictions obtained revealed new insights into the biology of tRNA^Sec.

Lokiarchaeota Marks the Transition between the Archaeal and Eukaryotic Selenocysteine Encoding Systems

Selenocysteine (Sec) is the 21st amino acid in the genetic code, inserted in response to UGA codons with the help of RNA structures, the SEC Insertion Sequence (SECIS) elements. The three domains of life feature distinct strategies for Sec insertion in proteins and its utilization. While bacteria and archaea possess similar sets of selenoproteins, Sec biosynthesis is more similar among archaea and eukaryotes. However, SECIS elements are completely different in the three domains of life. Here, we analyze the archaeon Lokiarchaeota that resolves the relationships among Sec insertion systems. This organism has selenoproteins representing five protein families, three of which have multiple Sec residues. Remarkably, these archaeal selenoprotein genes possess conserved RNA structures that strongly resemble the eukaryotic SECIS element, including key eukaryotic protein-binding sites. These structures also share similarity with the SECIS element in archaeal selenoprotein VhuD, suggesting a relation of direct descent. These results identify Lokiarchaeota as an intermediate form between the archaeal and eukaryotic Sec-encoding systems and clarify the evolution of the Sec insertion system.

Novel Two-Step Hierarchical Screening of Mutant Pools Reveals Mutants under Selection in Chicks

Contaminated chicken/egg products are major sources of human salmonellosis, yet the strategies used by Salmonella to colonize chickens are poorly understood. We applied a novel two-step hierarchical procedure to identify new genes important for colonization and persistence of Salmonella enterica serotype Typhimurium in chickens. A library of 182 S. Typhimurium mutants each containing a targeted deletion of a group of contiguous genes (for a total of 2,069 genes deleted) was used to identify regions under selection at 1, 3, and 9 days postinfection in chicks. Mutants in 11 regions were under selection at all assayed times (colonization mutants), and mutants in 15 regions were under selection only at day 9 (persistence mutants). We assembled a pool of 92 mutants, each deleted for a single gene, representing nearly all genes in nine regions under selection. Twelve single gene deletion mutants were under selection in this assay, and we confirmed 6 of 9 of these candidate mutants via competitive infections and complementation analysis in chicks. STM0580, STM1295, STM1297, STM3612, STM3615, and STM3734 are needed for Salmonella to colonize and persist in chicks and were not previously associated with this ability. One of these key genes, STM1297 (selD), is required for anaerobic growth and supports the ability to utilize formate under these conditions, suggesting that metabolism of formate is important during infection. We report a hierarchical screening strategy to interrogate large portions of the genome during infection of animals using pools of mutants of low complexity. Using this strategy, we identified six genes not previously known to be needed during infection in chicks, and one of these (STM1297) suggests an important role for formate metabolism during infection.

Genomic analysis of the nitrate-respiring Sphingopyxis granuli (formerly Sphingomonas macrogoltabida) strain TFA

Background

Sphingomonads are Alphaproteobacteria that belong to the Sphingomonas, Novosphingobium, Sphingopyxis or Sphingobium genera, They are physiologically diverse and broadly distributed in nature, playing important roles in oligotrophic environments and in the degradation of recalcitrant polyaromatic compounds, Sphingopyxis is a poorly studied genus of which only one representative (S. alaskensis RB2256) has been deeply characterized. In this paper we analyze the genomic features of S. granuli strain TFA (formerly Sphingomonas macrogoltabida) in comparison with the available Sphingopyxis sequenced genomes, to describe common characteristics of this genus and to highlight unique characteristics of strain TFA.

Results

The TFA genome has been assembled in a single circular chromosome of 4.7 Mb. Genomic sequence analysis and proteome comparison re-assigned the TFA strain to the Sphingopyxis genus and the S. granuli species. Some regions of the TFA genome show high similarity (ca. 100 %) to other bacteria and several genomic islands have been detected. Pathways for aromatic compound degradation have been predicted but no growth of TFA has been detected using these as carbon or nitrogen sources. Genes for nitrate respiration have been identified as TFA exclusive. Experimental data on anaerobic growth of TFA using nitrate as a terminal electron acceptor are also provided.

Conclusions

Sphingopyxis representatives form a compact phylogenetic group (with the exception of S. baekryungensis DSM 16222) that share several characteristics, such as being naturally resistant to streptomycin, having only one ribosomal operon, a low number of prophages and CRISPR sequences, absence of selenoproteins and presence of ectoin and other biosynthesis pathways for secondary metabolites. Moreover, the TFA genome organization shows evidence of the presence of putative integrative and conjugative elements (ICE) responsible for the acquisition of several characteristics by horizontal transfer mechanisms. Sphingopyxis representatives have been described as strict aerobes but anaerobic growth using nitrate as a terminal electron acceptor might confer an environmental advantage to the first S. granuli strain characterized at genomic level.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2411-1) contains supplementary material, which is available to authorized users.

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.

Probing the active site tryptophan of Staphylococcus aureus thioredoxin with an analog

Genetically encoded non-canonical amino acids are powerful tools of protein research and engineering; in particular they allow substitution of individual chemical groups or atoms in a protein of interest. One such amino acid is the tryptophan (Trp) analog 3-benzothienyl-l-alanine (Bta) with an imino-to-sulfur substitution in the five-membered ring. Unlike Trp, Bta is not capable of forming a hydrogen bond, but preserves other properties of a Trp residue. Here we present a pyrrolysyl-tRNA synthetase-derived, engineered enzyme BtaRS that enables efficient and site-specific Bta incorporation into proteins of interest in vivo. Furthermore, we report a 2.1 Å-resolution crystal structure of a BtaRS•Bta complex to show how BtaRS discriminates Bta from canonical amino acids, including Trp. To show utility in protein mutagenesis, we used BtaRS to introduce Bta to replace the Trp28 residue in the active site of Staphylococcus aureus thioredoxin. This experiment showed that not the hydrogen bond between residues Trp28 and Asp58, but the bulky aromatic side chain of Trp28 is important for active site maintenance. Collectively, our study provides a new and robust tool for checking the function of Trp in proteins.

Genetic code flexibility in microorganisms: novel mechanisms and impact on physiology

The genetic code, initially thought to be universal and immutable, is now known to contain many variations, including biased codon usage, codon reassignment, ambiguous decoding and recoding. As a result of recent advances in the areas of genome sequencing, biochemistry, bioinformatics and structural biology, our understanding of genetic code flexibility has advanced substantially in the past decade. In this Review, we highlight the prevalence, evolution and mechanistic basis of genetic code variations in microorganisms, and we discuss how this flexibility of the genetic code affects microbial physiology.

Evolution of the Selenoproteome in Helicobacter pylori and Epsilonproteobacteria

By competing for the acquisition of essential nutrients, Helicobacter pylori has the unique ability to persist in the human stomach, also causing nutritional insufficiencies in the host. Although the H. pylori genome apparently encodes selenocysteine synthase (SelA, HP1513), a key pyridoxal phosphate (PLP)-dependent enzyme for the incorporation of selenium into bacterial proteins, nothing is known about the use of this essential element in protein synthesis by this pathogen. We analyzed the evolution of the complete machinery for incorporation of selenium into proteins and the selenoproteome of several H. pylori strains and related Epsilonproteobacteria. Our searches identified the presence of selenoproteins-including the previously unknown DUF466 family-in various Epsilonproteobacteria, but not in H. pylori. We found that a complete system for selenocysteine incorporation was present in the Helicobacteriaceae ancestor and has been recently lost before the split of Helicobacter acinonychis and H. pylori. Our results indicate that H. pylori, at variance with other gastric and enterohepatic Helicobacter, does not use selenocysteine in protein synthesis and does not use selenium for tRNA wobble base modification. However, selA has survived as a functional gene, having lost the domain for the binding of selenocysteine tRNA, but maintaining the ability to bind the PLP cofactor. The evolutionary modifications described for the SelA protein of H. pylori find parallels in other bacterial and archaeal species, suggesting that an alternative enzymatic function is hidden in many proteins annotated as selenocysteinyl-tRNA synthase.