Proteomics-based refinement of Deinococcus deserti genome annotation reveals an unwonted use of non-canonical translation initiation codons

Novel model organisms and proteomics for a better biological understanding

The concept of « model organisms » is being revisited in the light of the latest advances in multi-omics technologies that can now capture the full range of molecular events that occur over time, regardless of the organism studied. Classic, well-studied models, such as Escherichia coli, Saccharomyces cerevisiae, to name a few, have long been valuable for hypothesis testing, reproducibility, and sharing common platforms among researchers. However, they are not suitable for all types of research. The complexity of unanswered questions in biology demands more elaborated systems, particularly to study plant and animal biodiversity, microbial ecosystems and their interactions with their hosts if any. More integrated systems, known as « holobionts », are emerging to describe and unify host organisms and associated microorganisms, providing an overview of all their possible interactions and trajectories. Comparative evolutionary proteomics offers interesting prospects for extrapolating knowledge from a few selected model organisms to others. This approach enables a deeper characterization of the diversity of proteins and proteoforms across the three branches of the tree of life, i.e. Bacteria, Archaea, and Eukarya. It also provides a powerful means to address remaining biological questions, such as identifying the key molecular players in organisms when they are confronted to environmental challenges, like anthropogenic toxicants, pathogens, dietary shifts or climate stressors, and proposing long-term sustainable solutions. SIGNIFICANCE: In this commentary, we reevaluated the concept of "model organisms" in light of advancements in multi-omics technologies. Traditional models have proven invaluable for hypothesis testing, reproducibility, and fostering shared research frameworks. However, we discussed that they are not universally applicable. To address complexities such as biodiversity and understand microbial ecosystems and their host interactions, integrated systems like "holobionts," which encompass host organisms and their associated microbes, are gaining prominence. Comparative evolutionary proteomics further enhances our understanding by enabling detailed exploration of protein diversity across organisms. This approach also facilitates the identification of critical molecular players in organisms facing environmental challenges, such as pollutants, pathogens, dietary changes, or climate stress, and contributes to developing sustainable long-term solutions.

A Proteogenomic Approach for the Identification of Virulence Factors in Leishmania Parasites

Identifying new genes involved in virulence and drug resistance may hold the key to a better understanding of parasitic diseases. The proteogenomic profiling of various Leishmania species, the causative agents of leishmaniasis, has identified several novel genes, N- and C-terminal extensions of proteins, and corrections of existing gene models. Various virulence factors (VFs) responsible for leishmaniasis have been previously annotated through a proteogenomic approach, including the C-terminal extension of heat shock protein 70 (HSP70). Furthermore, the diversity of VFs across Leishmania donovani, L. infantum, L. major, and L. mexicana was determined using phylogenetic analysis. Moreover, protein-protein interaction networks (PPINs) of VFs with HSPs aid in making significant biological interpretations. Overall, an integrated omics approach involving proteogenomics was used to identify and study the relationship among VFs with other interacting proteins, including HSPs. This chapter provides a step-by-step guide to the identification of new genes in Leishmania using a proteogenomic approach and their functional assignment using a bioinformatics-based approach.

The Construction of an Extreme Radiation-Resistant Perchlorate-Reducing Bacterium Using Deinococcus deserti Promoters

Perchlorate is one of the major inorganic pollutants in the natural environment and the living environment, which is toxic to organisms and difficult to degrade due to its special structure. As previously reported, the Phoenix Mars lander detected approximately 0.6% perchlorate in the Martian soil, indicating challenges for Earth-based life to survive there. Currently, biological approaches using dissimilatory perchlorate-reducing bacteria (DPRB) are the most promising methods for perchlorate degradation. However, the majority of DPRB exhibit limited radiation resistance, rendering them unsuitable for survival on Mars. In this study, we obtained the transcriptome data of Deinococcus deserti, and predicted and identified multiple constitutive expression promoters of D. deserti with varying activities. The top-five most active promoters were separately fused to specific genes involved in the degradation of perchlorate from DPRB Dechloromonas agitata CKB, and transformed into Deinococcus radiodurans R1, forming a novel dissimilatory perchlorate-reducing bacterium, R1-CKB. It exhibited both efficient perchlorate degradation capability and strong radiation resistance, potentially offering a valuable tool for the further enhancement of the Martian atmosphere in the future.

Taxonomical composition and functional analysis of biofilms sampled from a nuclear storage pool

Sampling small amounts of biofilm from harsh environments such as the biofilm present on the walls of a radioactive material storage pool offers few analytical options if taxonomic characterization and estimation of the different biomass contributions are the objectives. Although 16S/18S rRNA amplification on extracted DNA and sequencing is the most widely applied method, its reliability in terms of quantitation has been questioned as yields can be species-dependent. Here, we propose a tandem-mass spectrometry proteotyping approach consisting of acquiring peptide data and interpreting then against a generalist database without any a priori. The peptide sequence information is transformed into useful taxonomical information that allows to obtain the different biomass contributions at different taxonomical ranks. This new methodology is applied for the first time to analyze the composition of biofilms from minute quantities of material collected from a pool used to store radioactive sources in a nuclear facility. For these biofilms, we report the identification of three genera, namely Sphingomonas, Caulobacter, and Acidovorax, and their functional characterization by metaproteomics which shows that these organisms are metabolic active. Differential expression of Gene Ontology GOslim terms between the two main microorganisms highlights their metabolic specialization.

Experimental validation that human microbiome phages use alternative genetic coding

Previous bioinformatic analyses of metagenomic data have indicated that bacteriophages can use genetic codes different from those of their host bacteria. In particular, reassignment of stop codon TAG to glutamine (a variation known as 'genetic code 15') has been predicted. Here, we use LC-MS/MS-based metaproteomics of human fecal samples to provide experimental evidence of the use of genetic code 15 in two crAss-like phages. Furthermore, the proteomic data from several phage structural proteins supports the reassignment of the TAG stop codon to glutamine late in the phage infection cycle. Thus, our work experimentally validates the expression of genetic code 15 in human microbiome phages.

Deep N-terminomics of Mycobacterium tuberculosis H37Rv extensively correct annotated encoding genes

Mycobacterium tuberculosis (MTB) is a severe causing agent of tuberculosis (TB). Although H37Rv, the type strain of M. tuberculosis was sequenced in 1998, annotation errors of encoding genes have been frequently reported in hundreds of papers. This phenomenon is particularly severe at the 5' end of the genes. Here, we applied a TMPP [(N-Succinimidyloxycarbonylmethyl) tris (2,4,6-trimethoxyphenyl) phosphonium bromide] labeling combined with StageTip separating strategy on M. tuberculosis H37Rv to characterize the N-terminal start sites of its annotated encoding genes. Totally, 1047 proteins were identified with 2058 TMPP labeled N-terminal peptides from all the 2625 mass spectrometer (MS) sequenced proteins. Comparative genomics analysis allowed the re-annotation of 43 proteins' N-termini in H37Rv and 762 proteins in Mycobacteriaceae. All revised N-termini start sites were distributed in 5'-UTR of annotated genes due to over-annotation of previous N-terminal initiation codon, especially the ATG. In addition, we identified and verified a novel gene Rv1078A in +3 frame different from the annotated gene Rv1078 in +2 frame. Altogether, our findings contribute to the better understanding of N-terminal of H37Rv and other species from Mycobacteriaceae that can assist future studies on biological study.

A new era of radiation resistance bacteria in bioremediation and production of bioactive compounds with therapeutic potential and other aspects: An in-perspective review

Microorganisms that survive in extreme environmental conditions are known as 'extremophiles'. Recently, extremophiles draw an impression in biotechnology/pharmaceutical researches/industries because of their novel molecules, known as 'extremolytes'. The intriguing phenomenon of microbial radiation resistance probably arose independently throughout their evolution of selective pressures (e.g. UV, X-ray, Gamma radiation etc.). Radiation produces multiple types of damage/oxidation to nucleic acids, proteins and other crucial cellular components. Most of the literature on microbial radiation resistance is based on acute γ-irradiation experiments performed in the laboratory, typically involving pure cultures isolation and their application on bioremediation/therapeutic field. There is much less information other than bioremediation and therapeutic application of such promising microbes we called as 'new era'. Here we discus origin and diversity of radiation resistance bacteria as well as selective mechanisms by which microorganisms can sustain in radiation rich environment. Potential uses of these radiations resistant microbes in the field of bioremediation, bioactive compounds and therapeutic industry. Last but not the least, which is the new aspect of radiation resistance microbes. Our review suggest that resistance to chronic radiation is not limited to rare specialized strains from extreme environments, but can occur among common microbial taxa, perhaps due to overlap molecular mechanisms of resistance to radiation and other stressors. These stress tolerance potential make them potential for radionuclides remediation, their extremolytes can be useful as anti-oxidant and anti-proliferative agents. In current scenario they can be useful in various fields from natural dye synthesis to nanoparticles production and anti-cancer treatment.

Identification of Proteolysis Products in Protein Therapeutics through TMPP N-Terminal Tagging and Electron Transfer Dissociation Product Triggered Collisional Induced Dissociation Fragmentation

Thorough characterization of protein therapeutics is often challenging due to the heterogeneity arising from primary sequence variants, post-translational modifications, proteolytic clipping, or incomplete processing of the signal peptide. Modern mass spectrometry (MS) techniques are now routinely used to characterize such heterogeneous protein populations. Here, we present an LC-MS/MS method using (N-succinimidyloxycarbonylmethyl)-tris (2,4,6-trimethoxyphenyl) phosphonium bromide (TMPP-Ac-OSu) to label any free N-terminal α-amines to rapidly and selectively identify proteolytic clipping events. Electron transfer dissociation (ETD) fragmentation of these chemically tagged peptides generates two unique TMPP product ions, TMPP⁺ and TMPP-Ac-NH₂/c0. The presence of these signature ions following ETD is used to trigger subsequent collisional induced dissociation (CID) fragmentation of the precursor ion. This results in a small subset of CID tandem MS spectra that are used in a customized database search. Using a purified fusion monoclonal antibody (mAb) as an example, we demonstrate how TMPP labeling followed by ETD product ion triggered CID fragmentation is used to accurately identify two undesired clipping sites.

Bacterial riboproteogenomics: the era of N-terminal proteoform existence revealed

With the rapid increase in the number of sequenced prokaryotic genomes, relying on automated gene annotation became a necessity. Multiple lines of evidence, however, suggest that current bacterial genome annotations may contain inconsistencies and are incomplete, even for so-called well-annotated genomes. We here discuss underexplored sources of protein diversity and new methodologies for high-throughput genome reannotation. The expression of multiple molecular forms of proteins (proteoforms) from a single gene, particularly driven by alternative translation initiation, is gaining interest as a prominent contributor to bacterial protein diversity. In consequence, riboproteogenomic pipelines were proposed to comprehensively capture proteoform expression in prokaryotes by the complementary use of (positional) proteomics and the direct readout of translated genomic regions using ribosome profiling. To complement these discoveries, tailored strategies are required for the functional characterization of newly discovered bacterial proteoforms.

TaxIt: An Iterative Computational Pipeline for Untargeted Strain-Level Identification Using MS/MS Spectra from Pathogenic Single-Organism Samples

Untargeted accurate strain-level classification of a priori unidentified organisms using tandem mass spectrometry is a challenging task. Reference databases often lack taxonomic depth, limiting peptide assignments to the species level. However, the extension with detailed strain information increases runtime and decreases statistical power. In addition, larger databases contain a higher number of similar proteomes. We present TaxIt, an iterative workflow to address the increasing search space required for MS/MS-based strain-level classification of samples with unknown taxonomic origin. TaxIt first applies reference sequence data for initial identification of species candidates, followed by automated acquisition of relevant strain sequences for low level classification. Furthermore, proteome similarities resulting in ambiguous taxonomic assignments are addressed with an abundance weighting strategy to increase the confidence in candidate taxa. For benchmarking the performance of our method, we apply our iterative workflow on several samples of bacterial and viral origin. In comparison to noniterative approaches using unique peptides or advanced abundance correction, TaxIt identifies microbial strains correctly in all examples presented (with one tie), thereby demonstrating the potential for untargeted and deeper taxonomic classification. TaxIt makes extensive use of public, unrestricted, and continuously growing sequence resources such as the NCBI databases and is available under open-source BSD license at https://gitlab.com/rki_bioinformatics/TaxIt.

In Vino Veritas: An Invitation for Ambitious, Collaborative Proteogenomics Campaigns on Plant and Animal Models

Vitis vinifera has been an emblematic plant for humans since the Neolithic period. Human civilization has been shaped by its domestication as both its medicinal and nutritional values were exploited. It is now cultivated on all habitable continents, and more than 5000 varieties have been developed. A global passion for the art of wine fuels innovation and a profound desire for knowledge on this plant. The genome sequence of a homozygotic cultivar and several RNA-seq datasets on other varieties have been released paving the way to gaining further insight into its biology and tailoring improvements to varieties. However, its genome annotation remains unpolished. In this issue of Proteomics, Chapman and Bellgard (Proteomics 2017, 17, 1700197) discuss how proteogenomics can help improve genome annotation. By mining shotgun proteomics data, they defined new protein-coding genes, refined gene structures, and corrected numerous mRNA splicing events. This stimulating study shows how large international consortia could work together to improve plant and animal genome annotation on a large scale. To achieve this aim, time should be invested to generate comprehensive, high-quality experimental datasets for a wide range of well-defined lineages and exploit them with pipelines capable of handling giant datasets.

The Novel Anaerobiosis-Responsive Overlapping Gene ano Is Overlapping Antisense to the Annotated Gene ECs2385 of Escherichia coli O157:H7 Sakai

Current notion presumes that only one protein is encoded at a given bacterial genetic locus. However, transcription and translation of an overlapping open reading frame (ORF) of 186 bp length were discovered by RNAseq and RIBOseq experiments. This ORF is almost completely embedded in the annotated L,D-transpeptidase gene ECs2385 of Escherichia coli O157:H7 Sakai in the antisense reading frame -3. The ORF is transcribed as part of a bicistronic mRNA, which includes the annotated upstream gene ECs2384, encoding a murein lipoprotein. The transcriptional start site of the operon resides 38 bp upstream of the ECs2384 start codon and is driven by a predicted σ⁷⁰ promoter, which is constitutively active under different growth conditions. The bicistronic operon contains a ρ-independent terminator just upstream of the novel gene, significantly decreasing its transcription. The novel gene can be stably expressed as an EGFP-fusion protein and a translationally arrested mutant of ano, unable to produce the protein, shows a growth advantage in competitive growth experiments compared to the wild type under anaerobiosis. Therefore, the novel antisense overlapping gene is named ano (anaerobiosis responsive overlapping gene). A phylostratigraphic analysis indicates that ano originated very recently de novo by overprinting after the Escherichia/Shigella clade separated from other enterobacteria. Therefore, ano is one of the very rare cases of overlapping genes known in the genus Escherichia.

Proteogenomic insights into uranium tolerance of a Chernobyl's Microbacterium bacterial isolate

Microbacterium oleivorans A9 is a uranium-tolerant actinobacteria isolated from the trench T22 located near the Chernobyl nuclear power plant. This site is contaminated with different radionuclides including uranium. To observe the molecular changes at the proteome level occurring in this strain upon uranyl exposure and understand molecular mechanisms explaining its uranium tolerance, we established its draft genome and used this raw information to perform an in-depth proteogenomics study. High-throughput proteomics were performed on cells exposed or not to 10μM uranyl nitrate sampled at three previously identified phases of uranyl tolerance. We experimentally detected and annotated 1532 proteins and highlighted a total of 591 proteins for which abundances were significantly differing between conditions. Notably, proteins involved in phosphate and iron metabolisms show high dynamics. A large ratio of proteins more abundant upon uranyl stress, are distant from functionally-annotated known proteins, highlighting the lack of fundamental knowledge regarding numerous key molecular players from soil bacteria.

BIOLOGICAL SIGNIFICANCE

Microbacterium oleivorans A9 is an interesting environmental model to understand biological processes engaged in tolerance to radionuclides. Using an innovative proteogenomics approach, we explored its molecular mechanisms involved in uranium tolerance. We sequenced its genome, interpreted high-throughput proteomic data against a six-reading frame ORF database deduced from the draft genome, annotated the identified proteins and compared protein abundances from cells exposed or not to uranyl stress after a cascade search. These data show that a complex cellular response to uranium occurs in Microbacterium oleivorans A9, where one third of the experimental proteome is modified. In particular, the uranyl stress perturbed the phosphate and iron metabolic pathways. Furthermore, several transporters have been identified to be specifically associated to uranyl stress, paving the way to the development of biotechnological tools for uranium decontamination.

Identification of new protein-coding genes with a potential role in the virulence of the plant pathogen Xanthomonas euvesicatoria

Background

Bacteria of the genus Xanthomonas are economically important plant pathogens. Pathogenicity of Xanthomonas spp. depends on the type III-secretion system and additional virulence determinants. The number of sequenced Xanthomonas genomes increases rapidly, however, accurate annotation of these genomes is difficult, because it relies on gene prediction programs. In this study, we used a mass-spectrometry (MS)-based approach to identify the proteome of Xanthomonas euvesicatoria (Xe) strain 85–10 also known as X. campestris pv. vesicatoria, a well-studied member of plant-pathogenic Xanthomonadaceae.

Results

Using different culture conditions, MS-datasets were searched against a six-frame-translated genome database of Xe. In total, we identified 2588 proteins covering 55% of the Xe genome, including 764 hitherto hypothetical proteins. Our proteogenomic approach identified 30 new protein-coding genes and allowed correction of the N-termini of 50 protein-coding genes. For five novel and two N-terminally corrected genes the corresponding proteins were confirmed by immunoblot. Furthermore, our data indicate that two putative type VI-secretion systems encoded in Xe play no role in bacterial virulence which was experimentally confirmed.

Conclusions

The discovery and re-annotation of numerous genes in the genome of Xe shows that also a well-annotated genome can be improved. Additionally, our proteogenomic analyses validates “hypothetical” proteins and will improve annotation of Xanthomonadaceae genomes, providing a solid basis for further studies.

Electronic supplementary material

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

N-terminal Proteomics Assisted Profiling of the Unexplored Translation Initiation Landscape in Arabidopsis thaliana

Proteogenomics is an emerging research field yet lacking a uniform method of analysis. Proteogenomic studies in which N-terminal proteomics and ribosome profiling are combined, suggest that a high number of protein start sites are currently missing in genome annotations. We constructed a proteogenomic pipeline specific for the analysis of N-terminal proteomics data, with the aim of discovering novel translational start sites outside annotated protein coding regions. In summary, unidentified MS/MS spectra were matched to a specific N-terminal peptide library encompassing protein N termini encoded in the Arabidopsis thaliana genome. After a stringent false discovery rate filtering, 117 protein N termini compliant with N-terminal methionine excision specificity and indicative of translation initiation were found. These include N-terminal protein extensions and translation from transposable elements and pseudogenes. Gene prediction provided supporting protein-coding models for approximately half of the protein N termini. Besides the prediction of functional domains (partially) contained within the newly predicted ORFs, further supporting evidence of translation was found in the recently released Araport11 genome re-annotation of Arabidopsis and computational translations of sequences stored in public repositories. Most interestingly, complementary evidence by ribosome profiling was found for 23 protein N termini. Finally, by analyzing protein N-terminal peptides, an in silico analysis demonstrates the applicability of our N-terminal proteogenomics strategy in revealing protein-coding potential in species with well- and poorly-annotated genomes.

Conservation and diversity of the IrrE/DdrO-controlled radiation response in radiation-resistant Deinococcus bacteria

The extreme radiation resistance of Deinococcus bacteria requires the radiation‐stimulated cleavage of protein DdrO by a specific metalloprotease called IrrE. DdrO is the repressor of a predicted radiation/desiccation response (RDR) regulon, composed of radiation‐induced genes having a conserved DNA motif (RDRM) in their promoter regions. Here, we showed that addition of zinc ions to purified apo‐IrrE, and short exposure of Deinococcus cells to zinc ions, resulted in cleavage of DdrO in vitro and in vivo, respectively. Binding of IrrE to RDRM‐containing DNA or interaction of IrrE with DNA‐bound DdrO was not observed. The data are in line with IrrE being a zinc peptidase, and indicate that increased zinc availability, caused by oxidative stress, triggers the in vivo cleavage of DdrO unbound to DNA. Transcriptomics and proteomics of Deinococcus deserti confirmed the IrrE‐dependent regulation of predicted RDR regulon genes and also revealed additional members of this regulon. Comparative analysis showed that the RDR regulon is largely well conserved in Deinococcus species, but also showed diversity in the regulon composition. Notably, several RDR genes with an important role in radiation resistance in Deinococcus radiodurans, for example pprA, are not conserved in some other radiation‐resistant Deinococcus species.

Sharpening Host Defenses during Infection: Proteases Cut to the Chase

The human immune system consists of an intricate network of tightly controlled pathways, where proteases are essential instigators and executioners at multiple levels. Invading microbial pathogens also encode proteases that have evolved to manipulate and dysregulate host proteins, including host proteases during the course of disease. The identification of pathogen proteases as well as their substrates and mechanisms of action have empowered significant developments in therapeutics for infectious diseases. Yet for many pathogens, there remains a great deal to be discovered. Recently, proteomic techniques have been developed that can identify proteolytically processed proteins across the proteome. These “degradomics” approaches can identify human substrates of microbial proteases during infection in vivo and expose the molecular-level changes that occur in the human proteome during infection as an operational network to develop hypotheses for further research as well as new therapeutics. This Perspective Article reviews how proteases are utilized during infection by both the human host and invading bacterial pathogens, including archetypal virulence-associated microbial proteases, such as the Clostridia spp. botulinum and tetanus neurotoxins. We highlight the potential knowledge that degradomics studies of host–pathogen interactions would uncover, as well as how degradomics has been successfully applied in similar contexts, including use with a viral protease. We review how microbial proteases have been targeted in current therapeutic approaches and how microbial proteases have shaped and even contributed to human therapeutics beyond infectious disease. Finally, we discuss how, moving forward, degradomics research can greatly contribute to our understanding of how microbial pathogens cause disease in vivo and lead to the identification of novel substrates in vivo, and the development of improved therapeutics to counter these pathogens.

Natural mitochondrial proteolysis confirms transcription systematically exchanging/deleting nucleotides, peptides coded by expanded codons

Protein sequences have higher linguistic complexities than human languages. This indicates undeciphered multilayered, overprinted information/genetic codes. Some superimposed genetic information is revealed by detections of transcripts systematically (a) exchanging nucleotides (nine symmetric, e.g. A<->C, fourteen asymmetric, e.g. A->C->G->A, swinger RNAs) translated according to tri-, tetra- and pentacodons, and (b) deleting mono-, dinucleotides after each trinucleotide (delRNAs). Here analyses of two independent proteomic datasets considering natural proteolysis confirm independently translation of these non-canonical RNAs, also along tetra- and pentacodons, increasing coverage of putative, cryptically encoded proteins. Analyses assuming endoproteinase GluC and elastase digestions (cleavages after residues D, E, and A, L, I, V, respectively) detect additional peptides colocalizing with detected non-canonical RNAs. Analyses detect fewer peptides matching GluC-, elastase- than trypsin-digestions: artificial trypsin-digestion outweighs natural proteolysis. Results suggest occurrences of complete proteins entirely matching non-canonical, superimposed encoding(s). Protein-coding after bijective transformations could explain genetic code symmetries, such as along Rumer's transformation.

Reannotation of Genomes by Means of Proteomics Data

Omics approaches have become popular in biology as powerful discovery tools, and currently gain in interest for diagnostic applications. Establishing the accurate genome sequence of any organism is easy, but the outcome of its annotation by means of automatic pipelines remains imprecise. Some protein-encoding genes may be missed as soon as they are specific and poorly conserved in a given taxon, while important to explain the specific traits of the organism. Translational starts are also poorly predicted in a relatively important number of cases, thus impacting the protein sequence database used in proteomics, comparative genomics, and systems biology. The use of high-throughput proteomics data to improve genome annotation is an attractive option to obtain a more comprehensive molecular picture of a given organism. Here, protocols for reannotating prokaryote genomes are described based on shotgun proteomics and derivatization of protein N-termini with a positively charged reagent coupled to high-resolution tandem mass spectrometry.

GAPP: A Proteogenomic Software for Genome Annotation and Global Profiling of Post-translational Modifications in Prokaryotes

Although the number of sequenced prokaryotic genomes is growing rapidly, experimentally verified annotation of prokaryotic genome remains patchy and challenging. To facilitate genome annotation efforts for prokaryotes, we developed an open source software called GAPP for genome annotation and global profiling of post-translational modifications (PTMs) in prokaryotes. With a single command, it provides a standard workflow to validate and refine predicted genetic models and discover diverse PTM events. We demonstrated the utility of GAPP using proteomic data from Helicobacter pylori, one of the major human pathogens that is responsible for many gastric diseases. Our results confirmed 84.9% of the existing predicted H. pylori proteins, identified 20 novel protein coding genes, and corrected four existing gene models with regard to translation initiation sites. In particular, GAPP revealed a large repertoire of PTMs using the same proteomic data and provided a rich resource that can be used to examine the functions of reversible modifications in this human pathogen. This software is a powerful tool for genome annotation and global discovery of PTMs and is applicable to any sequenced prokaryotic organism; we expect that it will become an integral part of ongoing genome annotation efforts for prokaryotes. GAPP is freely available at https://sourceforge.net/projects/gappproteogenomic/.

Power of positive thinking in quantitative proteomics

Derivatization of proteins with specific isotope reagents has been widely explored for quantitative proteomics where the relative abundances of proteins present in different complex samples are compared by MS. This represents an interesting arena for innovation, where protein chemistry and MS are associated for the best of both worlds. Among the numerous reagents developed, those that introduce a permanent positive charge, such as (N-succinimidyloxycarbonylmethyl)-tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP), increase the ionizability of their targets and thus improve the sensitivity of the approach. TMPP labeling also modifies the hydrophobicity and changes the peptide fragmentation pattern. Because TMPP reacts preferably with the N-termini of proteins and peptides, its use has been explored for proteogenomics and de novo protein sequencing. In this issue of Proteomics, Shen et al. (Proteomics 2015, 15, 2903-2909) show that accurate quantitation of proteins can be obtained with light/heavy TMPP-labeling of peptides, which can be easily prepared and desalted in a homemade C8-SCX-C8 stagetip, and then monitored by nano-LC-MS/MS analysis. Their results demonstrate enhanced sequence coverage compared with other approaches. Combined with an efficient enrichment procedure, the higher sensitivity of this "positive attitude" reagent may facilitate much deeper investigations into the quantitative proteomics of complex samples.

Tissue-specific Proteogenomic Analysis of Plutella xylostella Larval Midgut Using a Multialgorithm Pipeline

The diamondback moth, Plutella xylostella (L.), is the major cosmopolitan pest of brassica and other cruciferous crops. Its larval midgut is a dynamic tissue that interfaces with a wide variety of toxicological and physiological processes. The draft sequence of the P. xylostella genome was recently released, but its annotation remains challenging because of the low sequence coverage of this branch of life and the poor description of exon/intron splicing rules for these insects. Peptide sequencing by computational assignment of tandem mass spectra to genome sequence information provides an experimental independent approach for confirming or refuting protein predictions, a concept that has been termed proteogenomics. In this study, we carried out an in-depth proteogenomic analysis to complement genome annotation of P. xylostella larval midgut based on shotgun HPLC-ESI-MS/MS data by means of a multialgorithm pipeline. A total of 876,341 tandem mass spectra were searched against the predicted P. xylostella protein sequences and a whole-genome six-frame translation database. Based on a data set comprising 2694 novel genome search specific peptides, we discovered 439 novel protein-coding genes and corrected 128 existing gene models. To get the most accurate data to seed further insect genome annotation, more than half of the novel protein-coding genes, i.e. 235 over 439, were further validated after RT-PCR amplification and sequencing of the corresponding transcripts. Furthermore, we validated 53 novel alternative splicings. Finally, a total of 6764 proteins were identified, resulting in one of the most comprehensive proteogenomic study of a nonmodel animal. As the first tissue-specific proteogenomics analysis of P. xylostella, this study provides the fundamental basis for high-throughput proteomics and functional genomics approaches aimed at deciphering the molecular mechanisms of resistance and controlling this pest.

Identification of Small Novel Coding Sequences, a Proteogenomics Endeavor

The identification of small proteins and peptides has consistently proven to be challenging. However, technological advances as well as multi-omics endeavors facilitate the identification of novel small coding sequences, leading to new insights. Specifically, the application of next generation sequencing technologies (NGS), providing accurate and sample specific transcriptome / translatome information, into the proteomics field led to more comprehensive results and new discoveries. This book chapter focuses on the inclusion of RNA-Seq and RIBO-Seq also known as ribosome profiling, an RNA-Seq based technique sequencing the +/- 30 bp long fragments captured by translating ribosomes. We emphasize the identification of micropeptides and neo-antigens, two distinct classes of small translation products, triggering our current understanding of biology. RNA-Seq is capable of capturing sample specific genomic variations, enabling focused neo-antigen identification. RIBO-Seq can identify translation events in small open reading frames which are considered to be non-coding, leading to the discovery of micropeptides. The identification of small translation products requires the integration of multi-omics data, stressing the importance of proteogenomics in this novel research area.

Proteogenomics of rare taxonomic phyla: A prospective treasure trove of protein coding genes

Sustainable innovations in sequencing technologies have resulted in a torrent of microbial genome sequencing projects. However, the prokaryotic genomes sequenced so far are unequally distributed along their phylogenetic tree; few phyla contain the majority, the rest only a few representatives. Accurate genome annotation lags far behind genome sequencing. While automated computational prediction, aided by comparative genomics, remains a popular choice for genome annotation, substantial fraction of these annotations are erroneous. Proteogenomics utilizes protein level experimental observations to annotate protein coding genes on a genome wide scale. Benefits of proteogenomics include discovery and correction of gene annotations regardless of their phylogenetic conservation. This not only allows detection of common, conserved proteins but also the discovery of protein products of rare genes that may be horizontally transferred or taxonomy specific. Chances of encountering such genes are more in rare phyla that comprise a small number of complete genome sequences. We collated all bacterial and archaeal proteogenomic studies carried out to date and reviewed them in the context of genome sequencing projects. Here, we present a comprehensive list of microbial proteogenomic studies, their taxonomic distribution, and also urge for targeted proteogenomics of underexplored taxa to build an extensive reference of protein coding genes.

The abundant and essential HU proteins in Deinococcus deserti and Deinococcus radiodurans are translated from leaderless mRNA

HU proteins have an important architectural role in nucleoid organization in bacteria. Compared with HU of many bacteria, HU proteins from Deinococcus species possess an N-terminal lysine-rich extension similar to the eukaryotic histone H1 C-terminal domain involved in DNA compaction. The single HU gene in Deinococcus radiodurans, encoding DrHU, is required for nucleoid compaction and cell viability. Deinococcus deserti contains three expressed HU genes, encoding DdHU1, DdHU2 and DdHU3. Here, we show that either DdHU1 or DdHU2 is essential in D. deserti. DdHU1 and DdHU2, but not DdHU3, can substitute for DrHU in D. radiodurans, indicating that DdHU3 may have a non-essential function different from DdHU1, DdHU2 and DrHU. Interestingly, the highly abundant DrHU and DdHU1 proteins, and also the less expressed DdHU2, are translated in Deinococcus from leaderless mRNAs, which lack a 5'-untranslated region and, hence, the Shine-Dalgarno sequence. Unexpectedly, cloning the DrHU or DdHU1 gene under control of a strong promoter in an expression plasmid, which results in leadered transcripts, strongly reduced the DrHU and DdHU1 protein level in D. radiodurans compared with that obtained from the natural leaderless gene. We also show that the start codon position for DrHU and DdHU1 should be reannotated, resulting in proteins that are 15 and 4 aa residues shorter than initially reported. The expression level and start codon correction were crucial for functional characterization of HU in Deinococcus.

Toward a standard in structural genome annotation for prokaryotes

Background

In an effort to identify the best practice for finding genes in prokaryotic genomes and propose it as a standard for automated annotation pipelines, 1,004,576 peptides were collected from various publicly available resources, and were used as a basis to evaluate various gene-calling methods. The peptides came from 45 bacterial replicons with an average GC content from 31 % to 74 %, biased toward higher GC content genomes. Automated, manual, and semi-manual methods were used to tally errors in three widely used gene calling methods, as evidenced by peptides mapped outside the boundaries of called genes.

Results

We found that the consensus set of identical genes predicted by the three methods constitutes only about 70 % of the genes predicted by each individual method (with start and stop required to coincide). Peptide data was useful for evaluating some of the differences between gene callers, but not reliable enough to make the results conclusive, due to limitations inherent in any proteogenomic study.

Conclusions

A single, unambiguous, unanimous best practice did not emerge from this analysis, since the available proteomics data were not adequate to provide an objective measurement of differences in the accuracy between these methods. However, as a result of this study, software, reference data, and procedures have been better matched among participants, representing a step toward a much-needed standard. In the absence of sufficient amount of exprimental data to achieve a universal standard, our recommendation is that any of these methods can be used by the community, as long as a single method is employed across all datasets to be compared.

Electronic supplementary material

The online version of this article (doi:10.1186/s40793-015-0034-9) contains supplementary material, which is available to authorized users.

Structural and functional characterization of two unusual endonuclease III enzymes from Deinococcus radiodurans

While most bacteria possess a single gene encoding the bifunctional DNA glycosylase Endonuclease III (EndoIII) in their genomes, Deinococcus radiodurans possesses three: DR2438 (DrEndoIII1), DR0289 (DrEndoIII2) and DR0982 (DrEndoIII3). Here we have determined the crystal structures of DrEndoIII1 and an N-terminally truncated form of DrEndoIII3 (DrEndoIII3Δ76). We have also generated a homology model of DrEndoIII2 and measured activity of the three enzymes. All three structures consist of two all α-helical domains, one of which exhibits a [4Fe-4S] cluster and the other a HhH-motif, separated by a DNA binding cleft, similar to previously determined structures of endonuclease III from Escherichia coli and Geobacillus stearothermophilus. However, both DrEndoIII1 and DrEndoIII3 possess an extended HhH motif with extra helical features and an altered electrostatic surface potential. In addition, the DNA binding cleft of DrEndoIII3 seems to be less accessible for DNA interactions, while in DrEndoIII1 it seems to be more open. Analysis of the enzyme activities shows that DrEndoIII2 is most similar to the previously studied enzymes, while DrEndoIII1 seems to be more distant with a weaker activity towards substrate DNA containing either thymine glycol or an abasic site. DrEndoIII3 is the most distantly related enzyme and displays no detectable activity towards these substrates even though the suggested catalytic residues are conserved. Based on a comparative structural analysis, we suggest that the altered surface potential, shape of the substrate-binding pockets and specific amino acid substitutions close to the active site and in the DNA interacting loops may underlie the unexpected differences in activity.

Genome and proteome analysis of Pseudomonas chloritidismutans AW-1T that grows on n-decane with chlorate or oxygen as electron acceptor

Growth of Pseudomonas chloritidismutans AW-1^T on C7 to C12 n-alkanes with oxygen or chlorate as electron acceptor was studied by genome and proteome analysis. Whole genome shotgun sequencing resulted in a 5 Mbp assembled sequence with a G + C content of 62.5%. The automatic annotation identified 4767 protein-encoding genes and a putative function could be assigned to almost 80% of the predicted proteins. The distinct phylogenetic position of P. chloritidismutans AW-1^T within the Pseudomonas stutzeri cluster became clear by comparison of average nucleotide identity values of sequenced genomes. Analysis of the proteome of P. chloritidismutans AW-1^T showed the versatility of this bacterium to adapt to aerobic and anaerobic growth conditions with acetate or n-decane as substrates. All enzymes involved in the alkane oxidation pathway were identified. An alkane monooxygenase was detected in n-decane-grown cells, but not in acetate-grown cells. The enzyme was found when grown in the presence of oxygen or chlorate, indicating that under both conditions an oxygenase-mediated pathway is employed for alkane degradation. Proteomic and biochemical data also showed that both chlorate reductase and chlorite dismutase are constitutively present, but most abundant under chlorate-reducing conditions.

Comparative genomics of Cluster O mycobacteriophages

Mycobacteriophages--viruses of mycobacterial hosts--are genetically diverse but morphologically are all classified in the Caudovirales with double-stranded DNA and tails. We describe here a group of five closely related mycobacteriophages--Corndog, Catdawg, Dylan, Firecracker, and YungJamal--designated as Cluster O with long flexible tails but with unusual prolate capsids. Proteomic analysis of phage Corndog particles, Catdawg particles, and Corndog-infected cells confirms expression of half of the predicted gene products and indicates a non-canonical mechanism for translation of the Corndog tape measure protein. Bioinformatic analysis identifies 8-9 strongly predicted SigA promoters and all five Cluster O genomes contain more than 30 copies of a 17 bp repeat sequence with dyad symmetry located throughout the genomes. Comparison of the Cluster O phages provides insights into phage genome evolution including the processes of gene flux by horizontal genetic exchange.

Crystal structure of the DNA polymerase III β subunit (β-clamp) from the extremophile Deinococcus radiodurans

BACKGROUND

Deinococcus radiodurans is an extremely radiation and desiccation resistant bacterium which can tolerate radiation doses up to 5,000 Grays without losing viability. We are studying the role of DNA repair and replication proteins for this unusual phenotype by a structural biology approach. The DNA polymerase III β subunit (β-clamp) acts as a sliding clamp on DNA, promoting the binding and processivity of many DNA-acting proteins, and here we report the crystal structure of D. radiodurans β-clamp (Drβ-clamp) at 2.0 Å resolution.

RESULTS

The sequence verification process revealed that at the time of the study the gene encoding Drβ-clamp was wrongly annotated in the genome database, encoding a protein of 393 instead of 362 amino acids. The short protein was successfully expressed, purified and used for crystallisation purposes in complex with Cy5-labeled DNA. The structure, which was obtained from blue crystals, shows a typical ring-shaped bacterial β-clamp formed of two monomers, each with three domains of identical topology, but with no visible DNA in electron density. A visualisation of the electrostatic surface potential reveals a highly negatively charged outer surface while the inner surface and the dimer forming interface have a more even charge distribution.

CONCLUSIONS

The structure of Drβ-clamp was determined to 2.0 Å resolution and shows an evenly distributed electrostatic surface charge on the DNA interacting side. We hypothesise that this charge distribution may facilitate efficient movement on encircled DNA and help ensure efficient DNA metabolism in D. radiodurans upon exposure to high doses of ionizing irradiation or desiccation.

RNA sequencing and proteogenomics reveal the importance of leaderless mRNAs in the radiation-tolerant bacterium Deinococcus deserti

Deinococcus deserti is a desiccation- and radiation-tolerant desert bacterium. Differential RNA sequencing (RNA-seq) was performed to explore the specificities of its transcriptome. Strikingly, for 1,174 (60%) mRNAs, the transcription start site was found exactly at (916 cases, 47%) or very close to the translation initiation codon AUG or GUG. Such proportion of leaderless mRNAs, which may resemble ancestral mRNAs, is unprecedented for a bacterial species. Proteomics showed that leaderless mRNAs are efficiently translated in D. deserti. Interestingly, we also found 173 additional transcripts with a 5′-AUG or 5′-GUG that would make them competent for ribosome binding and translation into novel small polypeptides. Fourteen of these are predicted to be leader peptides involved in transcription attenuation. Another 30 correlated with new gene predictions and/or showed conservation with annotated and nonannotated genes in other Deinococcus species, and five of these novel polypeptides were indeed detected by mass spectrometry. The data also allowed reannotation of the start codon position of 257 genes, including several DNA repair genes. Moreover, several novel highly radiation-induced genes were found, and their potential roles are discussed. On the basis of our RNA-seq and proteogenomics data, we propose that translation of many of the novel leaderless transcripts, which may have resulted from single-nucleotide changes and maintained by selective pressure, provides a new explanation for the generation of a cellular pool of small peptides important for protection of proteins against oxidation and thus for radiation/desiccation tolerance and adaptation to harsh environmental conditions.

Prioritizing targets for structural biology through the lens of proteomics: the archaeal protein TGAM_1934 from Thermococcus gammatolerans

ORFans are hypothetical proteins lacking any significant sequence similarity with other proteins. Here, we highlighted by quantitative proteomics the TGAM_1934 ORFan from the hyperradioresistant Thermococcus gammatolerans archaeon as one of the most abundant hypothetical proteins. This protein has been selected as a priority target for structure determination on the basis of its abundance in three cellular conditions. Its solution structure has been determined using multidimensional heteronuclear NMR spectroscopy. TGAM_1934 displays an original fold, although sharing some similarities with the 3D structure of the bacterial ortholog of frataxin, CyaY, a protein conserved in bacteria and eukaryotes and involved in iron-sulfur cluster biogenesis. These results highlight the potential of structural proteomics in prioritizing ORFan targets for structure determination based on quantitative proteomics data. The proteomic data and structure coordinates have been deposited to the ProteomeXchange with identifier PXD000402 (http://proteomecentral.proteomexchange.org/dataset/PXD000402) and Protein Data Bank under the accession number 2mcf, respectively.

Identification of the over alkylation sites of a protein by IAM in MALDI-TOF/TOF tandem mass spectrometry

Overalkylation often appears during the proteolytic digestion process when using iodoacetamide (IAM) to protect the produced side chain thiol of Cys from disulfide bonds.

Proteogenomics in microbiology: taking the right turn at the junction of genomics and proteomics

High-accuracy and high-throughput proteomic methods have completely changed the way we can identify and characterize proteins. MS-based proteomics can now provide a unique supplement to genomic data and add a new level of information to the interpretation of genomic sequences. Proteomics-driven genome annotation has become especially relevant in microbiology where genomes are sequenced on a daily basis and limitations of an in silico driven annotation process are well recognized. In this review paper, we outline different strategies on how one can design a proteogenomic experiment, for example on genome-sequenced (synonymous proteogenomics) versus unsequenced organisms (ortho-proteogenomics) or with the aid of other "omic" data such as RNA-seq. We touch upon many challenges that are encountered during a typical proteogenomic study, mostly concerning bioinformatics methods and downstream data analysis, but also related to creation and use of sequence databases. A large list of proteogenomic case studies of different microorganisms is provided to illustrate the mapping of MS/MS-derived peptide spectra to genomic DNA sequences. These investigations have led to accurate determination of translational initiation sites, pointed out eventual read-throughs or programmed frameshifts, detected signal peptide processing or other protein maturation events, removed questionable annotation assignments, and provided evidence for predicted hypothetical proteins.

N-terminomics and proteogenomics, getting off to a good start

Proteogenomics consists of the annotation or reannotation of protein-coding nucleic acid sequences based on the empirical observation of their gene products. While functional annotation of predicted genes is increasingly feasible given the multiplicity of genomes available for many branches of the tree of life, the accurate annotation of the translational start sites is still a point of contention. Extensive coverage of the proteome, including specifically the N-termini, is now possible, thanks to next-generation mass spectrometers able to record data from thousands of proteins at once. Efforts to increase the peptide coverage of protein sequences and to detect low abundance proteins are important to make proteomic and proteogenomic studies more comprehensive. In this review, we present the panoply of N-terminus-oriented strategies that have been developed over the last decade.

Integrating genomic, transcriptomic, and interactome data to improve Peptide and protein identification in shotgun proteomics

Mass spectrometry (MS)-based shotgun proteomics is an effective technology for global proteome profiling. The ultimate goal is to assign tandem MS spectra to peptides and subsequently infer proteins and their abundance. In addition to database searching and protein assembly algorithms, computational approaches have been developed to integrate genomic, transcriptomic, and interactome information to improve peptide and protein identification. Earlier efforts focus primarily on making databases more comprehensive using publicly available genomic and transcriptomic data. More recently, with the increasing affordability of the Next Generation Sequencing (NGS) technologies, personalized protein databases derived from sample-specific genomic and transcriptomic data have emerged as an attractive strategy. In addition, incorporating interactome data not only improves protein identification but also puts identified proteins into their functional context and thus facilitates data interpretation. In this paper, we survey the major integrative bioinformatics approaches that have been developed during the past decade and discuss their merits and demerits.

Assessing the exoproteome of marine bacteria, lesson from a RTX-toxin abundantly secreted by Phaeobacter strain DSM 17395

Bacteria from the Roseobacter clade are abundant in surface marine ecosystems as over 10% of bacterial cells in the open ocean and 20% in coastal waters belong to this group. In order to document how these marine bacteria interact with their environment, we analyzed the exoproteome of Phaeobacter strain DSM 17395. We grew the strain in marine medium, collected the exoproteome and catalogued its content with high-throughput nanoLC-MS/MS shotgun proteomics. The major component represented 60% of the total protein content but was refractory to either classical proteomic identification or proteogenomics. We de novo sequenced this abundant protein with high-resolution tandem mass spectra which turned out being the 53 kDa RTX-toxin ZP_02147451. It comprised a peptidase M10 serralysin domain. We explained its recalcitrance to trypsin proteolysis and proteomic identification by its unusual low number of basic residues. We found this is a conserved trait in RTX-toxins from Roseobacter strains which probably explains their persistence in the harsh conditions around bacteria. Comprehensive analysis of exoproteomes from environmental bacteria should take into account this proteolytic recalcitrance.

N-Terminal-oriented proteogenomics of the marine bacterium roseobacter denitrificans Och114 using N-Succinimidyloxycarbonylmethyl)tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP) labeling and diagonal chromatography

Given the ease of whole genome sequencing with next-generation sequencers, structural and functional gene annotation is now purely based on automated prediction. However, errors in gene structure are frequent, the correct determination of start codons being one of the main concerns. Here, we combine protein N termini derivatization using (N-Succinimidyloxycarbonylmethyl)tris(2,4,6-trimethoxyphenyl)phosphonium bromide (TMPP Ac-OSu) as a labeling reagent with the COmbined FRActional DIagonal Chromatography (COFRADIC) sorting method to enrich labeled N-terminal peptides for mass spectrometry detection. Protein digestion was performed in parallel with three proteases to obtain a reliable automatic validation of protein N termini. The analysis of these N-terminal enriched fractions by high-resolution tandem mass spectrometry allowed the annotation refinement of 534 proteins of the model marine bacterium Roseobacter denitrificans OCh114. This study is especially efficient regarding mass spectrometry analytical time. From the 534 validated N termini, 480 confirmed existing gene annotations, 41 highlighted erroneous start codon annotations, five revealed totally new mis-annotated genes; the mass spectrometry data also suggested the existence of multiple start sites for eight different genes, a result that challenges the current view of protein translation initiation. Finally, we identified several proteins for which classical genome homology-driven annotation was inconsistent, questioning the validity of automatic annotation pipelines and emphasizing the need for complementary proteomic data. All data have been deposited to the ProteomeXchange with identifier PXD000337.

Moving from unsequenced to sequenced genome: reanalysis of the proteome of Leishmania donovani

The kinetoplastid protozoan parasite, Leishmania donovani, is the causative agent of kala azar or visceral leishmaniasis. Kala azar is a severe form of leishmaniasis that is fatal in the majority of untreated cases. Studies on proteomic analysis of L. donovani thus far have been carried out using homology-based identification based on related Leishmania species (L. infantum, L. major and L. braziliensis) whose genomes have been sequenced. Recently, the genome of L. donovani was fully sequenced and the data became publicly available. We took advantage of the availability of its genomic sequence to carry out a more accurate proteogenomic analysis of L. donovani proteome using our previously generated dataset. This resulted in identification of 17,504 unique peptides upon database-dependent search against the annotated proteins in L. donovani. These peptides were assigned to 3999 unique proteins in L. donovani. 2296 proteins were identified in both the life stages of L. donovani, while 613 and 1090 proteins were identified only from amastigote and promastigote stages, respectively. The proteomic data was also searched against six-frame translated L. donovani genome, which led to 255 genome search-specific peptides (GSSPs) resulting in identification of 20 novel genes and correction of 40 existing gene models in L. donovani.

BIOLOGICAL SIGNIFICANCE

Leishmania donovani genome sequencing was recently completed, which permitted us to use a proteogenomic approach to map its proteome and to carry out annotation of it genome. This resulted in mapping of 50% (3999 proteins) of L. donovani proteome. Our study identified 20 novel genes previously not predicted from the L. donovani genome in addition to correcting annotations of 40 existing gene models. The identified proteins may help in better understanding of stage-specific protein expression profiles in L. donovani and to identify novel stage-specific drug targets in L. donovani which could be used in the treatment of leishmaniasis. This article is part of a Special Issue entitled: Trends in Microbial Proteomics.