Proteogenomics: concepts, applications and computational strategies

Application of Proteomics Technology Based on LC-MS Combined with Western Blotting and Co-IP in Antiviral Innate Immunity

As an interferon-stimulating factor protein, STING plays a role in the response and downstream liaison in antiviral natural immunity. Upon viral invasion, the immediate response of STING protein leads to a series of changes in downstream proteins, which ultimately leads to an antiviral immune response in the form of proinflammatory cytokines and type I interferons, thus triggering an innate immune response, an adaptive immune response in vivo, and long-term protection of the host. In the field of antiviral natural immunity, it is particularly important to rigorously and sequentially probe the dynamic changes in the antiviral natural immunity connector protein STING caused by the entire anti-inflammatory and anti-pathway mechanism and the differences in upstream and downstream proteins. Traditionally, proteomics technology has been validated by detecting proteins in a 2D platform, for which it is difficult to sensitively identify changes in the nature and abundance of target proteins. With the development of mass spectrometry (MS) technology, MS-based proteomics has made important contributions to characterizing the dynamic changes in the natural immune proteome induced by viral infections. MS analytical techniques have several advantages, such as high throughput, rapidity, sensitivity, accuracy, and automation. The most common techniques for detecting complex proteomes are liquid chromatography (LC) and mass spectrometry (MS). LC-MS (Liquid Chromatography-Mass Spectrometry), which combines the physical separation capability of LC and the mass analysis capability of MS, is a powerful technique mainly used for analyzing the proteome of cells, tissues, and body fluids. To explore the combination of traditional proteomics techniques such as Western blotting, Co-IP (co-Immunoprecipitation), and the latest LC-MS methods to probe the anti-inflammatory pathway and the differential changes in upstream and downstream proteins induced by the antiviral natural immune junction protein STING.

Systematic data analysis pipeline for quantitative morphological cell phenotyping

Quantitative morphological phenotyping (QMP) is an image-based method used to capture morphological features at both the cellular and population level. Its interdisciplinary nature, spanning from data collection to result analysis and interpretation, can lead to uncertainties, particularly among those new to this actively growing field. High analytical specificity for a typical QMP is achieved through sophisticated approaches that can leverage subtle cellular morphological changes. Here, we outline a systematic workflow to refine the QMP methodology. For a practical review, we describe the main steps of a typical QMP; in each step, we discuss the available methods, their applications, advantages, and disadvantages, along with the R functions and packages for easy implementation. This review does not cover theoretical backgrounds, but provides several references for interested researchers. It aims to broaden the horizons for future phenome studies and demonstrate how to exploit years of endeavors to achieve more with less.

PgxSAVy: A tool for comprehensive evaluation of variant peptide quality in proteogenomics - catching the (un)usual suspects

Variant peptides resulting from single nucleotide polymorphisms (SNPs) can lead to aberrant protein functions and have translational potential for disease diagnosis and personalized therapy. Variant peptides detected by proteogenomics are fraught with high number of false positives, but there is no uniform and comprehensive approach to assess variant quality across analysis pipelines. Despite class-specific FDR along with ad-hoc filters, the problem is far from solved. These protocols are typically manual and tedious, and thus not uniform across labs. We demonstrate that variant peptide rescoring, integrated with intensity, variant event information and search result features, allows better discrimination of correct variant peptides. Implemented into PgxSAVy - a tool for quality control of variant peptides, this method can tackle the high rate of false positives. PgxSAVy provides a rigorous framework for quality control and annotations of variant peptides on the basis of (i) variant quality, (ii) isobaric masses, and (iii) disease annotation. PgxSAVy demonstrated high accuracy by identifying true variants with 98.43% accuracy on simulated data. Large-scale proteogenomic reanalysis of ∼2.8 million spectra (PXD004010 and PXD001468) resulted in 12,705 variant peptide spectrum matches (PSMs), of which PgxSAVy evaluated 3028 (23.8%), 1409 (11.1%) and 8268 (65.1%) as confident, semi-confident and doubtful respectively. PgxSAVy also annotates the variants based on their pathogenicity and provides support for assisted manual validation. The analysis of proteins carrying variants can provide fine granularity in discovering important pathways. PgxSAVy will advance personalized medicine by providing a comprehensive framework for quality control and prioritization of proteogenomics variants. PgxSAVy is freely available at https://pgxsavy.igib.res.in/ as a webserver and https://github.com/anuragraj/PgxSAVy as a stand-alone tool.

Myoglobin as a molecular biomarker for meat authentication and traceability

The global incidence of economically motivated meat adulteration represents a crucial issue for the food industry. Undeclared addition of cheaper or low-quality species to meat products of high commercial value has become a common practice that needs to be countered with specific measures. In this framework, myoglobin (Mb) is a sarcoplasmic haemoprotein, primarily responsible for meat colour and has been successfully used in meat fraud authentication. Mb is highly soluble in water, easily monitored at 409 nm and species-specific. Knowing that various analytical DNA-based and protein-based methods, as well as spectroscopic techniques have been developed over the years for the detection of meat fraud, the aim of the present review is to take stock of the situation regarding the possible use of Mb as a molecular biomarker for the easy and rapid detection of undeclared species in meat products, avoiding the need of sophisticated or expensive equipment and specialised operators.

Proteogenomics-guided functional venomics resolves the toxin arsenal and activity of Deinagkistrodon acutus venom

Snakebite primarily impacts rural communities of Africa, Asia, and Latin America. The sharp-nosed viper (Deinagkistrodon acutus) is among the snakes of highest medical importance in Asia. Despite various studies on its venom using modern venomics techniques, a comprehensive understanding of composition and function of this species' venom remains lacking. We combined proteogenomics with extensive bioactivity profiling to present the first genome-level catalogue of D. acutus venom proteins and their exochemistry. Our analysis identified an unusually simple venom containing 45 components from 20 distinct protein families. Relative toxin abundances indicate that C-type lectin and C-type lectin-related protein (CTL), snake venom metalloproteinase (svMP), snake venom serine protease (svSP), and phospholipase A2 (PLA2) constitute 90 % of the venom. Bioassays targeting key aspects of viperid envenomation showed considerable concentration-dependent cytotoxicity, particularly in kidney and lung cells, and potent protease and PLA2 activity. Factor Xa and thrombin activities were minor, and no plasmin activity was observed. Effects on haemolysis, intracellular calcium (Ca2+) release, and nitric oxide (NO) synthesis were negligible. Our analysis provides the first holistic genome-based overview of the toxin arsenal of D. acutus, predicting the molecular and functional basis of its life-threatening effects, and opens novel avenues for treating envenomation by this highly dangerous snake.

Phenotyping Tumor Heterogeneity through Proteogenomics: Study Models and Challenges

Tumor heterogeneity refers to the diversity observed among tumor cells: both between different tumors (inter-tumor heterogeneity) and within a single tumor (intra-tumor heterogeneity). These cells can display distinct morphological and phenotypic characteristics, including variations in cellular morphology, metastatic potential and variability treatment responses among patients. Therefore, a comprehensive understanding of such heterogeneity is necessary for deciphering tumor-specific mechanisms that may be diagnostically and therapeutically valuable. Innovative and multidisciplinary approaches are needed to understand this complex feature. In this context, proteogenomics has been emerging as a significant resource for integrating omics fields such as genomics and proteomics. By combining data obtained from both Next-Generation Sequencing (NGS) technologies and mass spectrometry (MS) analyses, proteogenomics aims to provide a comprehensive view of tumor heterogeneity. This approach reveals molecular alterations and phenotypic features related to tumor subtypes, potentially identifying therapeutic biomarkers. Many achievements have been made; however, despite continuous advances in proteogenomics-based methodologies, several challenges remain: in particular the limitations in sensitivity and specificity and the lack of optimal study models. This review highlights the impact of proteogenomics on characterizing tumor phenotypes, focusing on the critical challenges and current limitations of its use in different clinical and preclinical models for tumor phenotypic characterization.

Evidence for widespread translation of 5' untranslated regions

Ribosome profiling experiments support the translation of a range of novel human open reading frames. By contrast, most peptides from large-scale proteomics experiments derive from just one source, 5' untranslated regions. Across the human genome we find evidence for 192 translated upstream regions, most of which would produce protein isoforms with extended N-terminal ends. Almost all of these N-terminal extensions are from highly abundant genes, which suggests that the novel regions we detect are just the tip of the iceberg. These upstream regions have characteristics that are not typical of coding exons. Their GC-content is remarkably high, even higher than 5' regions in other genes, and a large majority have non-canonical start codons. Although some novel upstream regions have cross-species conservation - five have orthologues in invertebrates for example - the reading frames of two thirds are not conserved beyond simians. These non-conserved regions also have no evidence of purifying selection, which suggests that much of this translation is not functional. In addition, non-conserved upstream regions have significantly more peptides in cancer cell lines than would be expected, a strong indication that an aberrant or noisy translation initiation process may play an important role in translation from upstream regions.

Rp3: Ribosome profiling-assisted proteogenomics improves coverage and confidence during microprotein discovery

There has been a dramatic increase in the identification of non-canonical translation and a significant expansion of the protein-coding genome. Among the strategies used to identify unannotated small Open Reading Frames (smORFs) that encode microproteins, Ribosome profiling (Ribo-Seq) is the gold standard for the annotation of novel coding sequences by reporting on smORF translation. In Ribo-Seq, ribosome-protected footprints (RPFs) that map to multiple genomic sites are removed since they cannot be unambiguously assigned to a specific genomic location. Furthermore, RPFs necessarily result in short (25-34 nucleotides) reads, increasing the chance of multi-mapping alignments, such that smORFs residing in these regions cannot be identified by Ribo-Seq. Moreover, it has been challenging to identify protein evidence for Ribo-Seq. To solve this, we developed Rp3, a pipeline that integrates proteogenomics and Ribosome profiling to provide unambiguous evidence for a subset of microproteins missed by current Ribo-Seq pipelines. Here, we show that Rp3 maximizes proteomics detection and confidence of microprotein-encoding smORFs.

Community Resource: Large-Scale Proteogenomics to Refine Wheat Genome Annotations

Triticum aestivum is an important crop whose reference genome (International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v2.1) offers a valuable resource for understanding wheat genetic structure, improving agronomic traits, and developing new cultivars. A key aspect of gene model annotation is protein-level evidence of gene expression obtained from proteomics studies, followed up by proteogenomics to physically map proteins to the genome. In this research, we have retrieved the largest recent wheat proteomics datasets publicly available and applied the Basic Local Alignment Search Tool (tBLASTn) algorithm to map the 861,759 identified unique peptides against IWGSC RefSeq v2.1. Of the 92,719 hits, 83,015 unique peptides aligned along 33,612 High Confidence (HC) genes, thus validating 31.4% of all wheat HC gene models. Furthermore, 6685 unique peptides were mapped against 3702 Low Confidence (LC) gene models, and we argue that these gene models should be considered for HC status. The remaining 2934 orphan peptides can be used for novel gene discovery, as exemplified here on chromosome 4D. We demonstrated that tBLASTn could not map peptides exhibiting mid-sequence frame shift. We supply all our proteogenomics results, Galaxy workflow and Python code, as well as Browser Extensible Data (BED) files as a resource for the wheat community via the Apollo Jbrowse, and GitHub repositories. Our workflow could be applied to other proteomics datasets to expand this resource with proteins and peptides from biotically and abiotically stressed samples. This would help tease out wheat gene expression under various environmental conditions, both spatially and temporally.

Improvement of Tumor Neoantigen Detection by High-Field Asymmetric Waveform Ion Mobility Mass Spectrometry

Cancer neoantigens have been shown to elicit cancer-specific T-cell responses and have garnered much attention for their roles in both spontaneous and therapeutically induced antitumor responses. Mass spectrometry (MS) profiling of tumor immunopeptidomes has been used, in part, to identify MHC-bound mutant neoantigen ligands. However, under standard conditions, MS-based detection of such rare but clinically relevant neoantigens is relatively insensitive, requiring 300 million cells or more. Here, to quantitatively define the minimum detectable amounts of therapeutically relevant MHC-I and MHC-II neoantigen peptides, we analyzed different dilutions of immunopeptidomes isolated from the well-characterized T3 mouse methylcholanthrene (MCA)-induced cell line by MS. Using either data-dependent acquisition or parallel reaction monitoring (PRM), we established the minimum amount of material required to detect the major T3 neoantigens in the presence or absence of high field asymmetric waveform ion mobility spectrometry (FAIMS). This analysis yielded a 14-fold enhancement of sensitivity in detecting the major T3 MHC-I neoantigen (mLama4) with FAIMS-PRM compared with PRM without FAIMS, allowing ex vivo detection of this neoantigen from an individual 100 mg T3 tumor. These findings were then extended to two other independent MCA-sarcoma lines (1956 and F244). This study demonstrates that FAIMS substantially increases the sensitivity of MS-based characterization of validated neoantigens from tumors.

Buckwheat OMICS: present status and future prospects

Buckwheat (Fagopyrum spp.) is an underutilized resilient crop of North Western Himalayas belonging to the family Polygonaceae and is a source of essential nutrients and therapeutics. Common Buckwheat and Tatary Buckwheat are the two main cultivated species used as food. It is the only grain crop possessing rutin, an important metabolite with high nutraceutical potential. Due to its inherent tolerance to various biotic and abiotic stresses and a short life cycle, Buckwheat has been proposed as a model crop plant. Nutritional security is one of the major concerns, breeding for a nutrient-dense crop such as Buckwheat will provide a sustainable solution. Efforts toward improving Buckwheat for nutrition and yield are limited due to the lack of available: genetic resources, genomics, transcriptomics and metabolomics. In order to harness the agricultural importance of Buckwheat, an integrated breeding and OMICS platforms needs to be established that can pave the way for a better understanding of crop biology and developing commercial varieties. This, coupled with the availability of the genome sequences of both Buckwheat species in the public domain, should facilitate the identification of alleles/QTLs and candidate genes. There is a need to further our understanding of the molecular basis of the genetic regulation that controls various economically important traits. The present review focuses on: the food and nutritional importance of Buckwheat, its various omics resources, utilization of omics approaches in understanding Buckwheat biology and, finally, how an integrated platform of breeding and omics will help in developing commercially high yielding nutrient rich cultivars in Buckwheat.

Advances in the Clinical Application of High-throughput Proteomics

High-throughput proteomics has become an exciting field and a potential frontier of modern medicine since the early 2000s. While significant progress has been made in the technical aspects of the field, translating proteomics to clinical applications has been challenging. This review summarizes recent advances in clinical applications of high-throughput proteomics and discusses the associated challenges, advantages, and future directions. We focus on research progress and clinical applications of high-throughput proteomics in breast cancer, bladder cancer, laryngeal squamous cell carcinoma, gastric cancer, colorectal cancer, and coronavirus disease 2019. The future application of high-throughput proteomics will face challenges such as varying protein properties, limitations of statistical modeling, technical and logistical difficulties in data deposition, integration, and harmonization, as well as regulatory requirements for clinical validation and considerations. However, there are several noteworthy advantages of high-throughput proteomics, including the identification of novel global protein networks, the discovery of new proteins, and the synergistic incorporation with other omic data. We look forward to participating in and embracing future advances in high-throughput proteomics, such as proteomics-based single-cell biology and its clinical applications, individualized proteomics, pathology informatics, digital pathology, and deep learning models for high-throughput proteomics. Several new proteomic technologies are noteworthy, including data-independent acquisition mass spectrometry, nanopore-based proteomics, 4-D proteomics, and secondary ion mass spectrometry. In summary, we believe high-throughput proteomics will drastically shift the paradigm of translational research, clinical practice, and public health in the near future.

Understanding bacterial pathogen diversity: A proteogenomic analysis and use of an array of genome assemblies to identify novel virulence factors of the honey bee bacterial pathogen Paenibacillus larvae

Mass spectrometry proteomics data are typically evaluated against publicly available annotated sequences, but the proteogenomics approach is a useful alternative. A single genome is commonly utilized in custom proteomic and proteogenomic data analysis. We pose the question of whether utilizing numerous different genome assemblies in a search database would be beneficial. We reanalyzed raw data from the exoprotein fraction of four reference Enterobacterial Repetitive Intergenic Consensus (ERIC) I-IV genotypes of the honey bee bacterial pathogen Paenibacillus larvae and evaluated them against three reference databases (from NCBI-protein, RefSeq, and UniProt) together with an array of protein sequences generated by six-frame direct translation of 15 genome assemblies from GenBank. The wide search yielded 453 protein hits/groups, which UpSet analysis categorized into 50 groups based on the success of protein identification by the 18 database components. Nine hits that were not identified by a unique peptide were not considered for marker selection, which discarded the only protein that was not identified by the reference databases. We propose that the variability in successful identifications between genome assemblies is useful for marker mining. The results suggest that various strains of P. larvae can exhibit specific traits that set them apart from the established genotypes ERIC I-V.

Rescoring Peptide Spectrum Matches: Boosting Proteomics Performance by Integrating Peptide Property Predictors Into Peptide Identification

Rescoring of peptide spectrum matches originating from database search engines enabled by peptide property predictors is exceeding the performance of peptide identification from traditional database search engines. In contrast to the peptide spectrum match scores calculated by traditional database search engines, rescoring peptide spectrum matches generates scores based on comparing observed and predicted peptide properties, such as fragment ion intensities and retention times. These newly generated scores enable a more efficient discrimination between correct and incorrect peptide spectrum matches. This approach was shown to lead to substantial improvements in the number of confidently identified peptides, facilitating the analysis of challenging datasets in various fields such as immunopeptidomics, metaproteomics, proteogenomics, and single-cell proteomics. In this review, we summarize the key elements leading up to the recent introduction of multiple data-driven rescoring pipelines. We provide an overview of relevant post-processing rescoring tools, introduce prominent data-driven rescoring pipelines for various applications, and highlight limitations, opportunities, and future perspectives of this approach and its impact on mass spectrometry-based proteomics.

Data-Independent Acquisition: A Milestone and Prospect in Clinical Mass Spectrometry-Based Proteomics

Data-independent acquisition (DIA) has revolutionized the field of mass spectrometry (MS)-based proteomics over the past few years. DIA stands out for its ability to systematically sample all peptides in a given m/z range, allowing an unbiased acquisition of proteomics data. This greatly mitigates the issue of missing values and significantly enhances quantitative accuracy, precision, and reproducibility compared to many traditional methods. This review focuses on the critical role of DIA analysis software tools, primarily focusing on their capabilities and the challenges they address in proteomic research. Advances in MS technology, such as trapped ion mobility spectrometry, or high field asymmetric waveform ion mobility spectrometry require sophisticated analysis software capable of handling the increased data complexity and exploiting the full potential of DIA. We identify and critically evaluate leading software tools in the DIA landscape, discussing their unique features, and the reliability of their quantitative and qualitative outputs. We present the biological and clinical relevance of DIA-MS and discuss crucial publications that paved the way for in-depth proteomic characterization in patient-derived specimens. Furthermore, we provide a perspective on emerging trends in clinical applications and present upcoming challenges including standardization and certification of MS-based acquisition strategies in molecular diagnostics. While we emphasize the need for continuous development of software tools to keep pace with evolving technologies, we advise researchers against uncritically accepting the results from DIA software tools. Each tool may have its own biases, and some may not be as sensitive or reliable as others. Our overarching recommendation for both researchers and clinicians is to employ multiple DIA analysis tools, utilizing orthogonal analysis approaches to enhance the robustness and reliability of their findings.

Integration of proteomics in the molecular tumor board

Cancer remains one of the most complex and challenging diseases in mankind. To address the need for a personalized treatment approach for particularly complex tumor cases, molecular tumor boards (MTBs) have been initiated. MTBs are interdisciplinary teams that perform in-depth molecular diagnostics to cooperatively and interdisciplinarily advise on the best therapeutic strategy. Current molecular diagnostics are routinely performed on the transcriptomic and genomic levels, aiming to identify tumor-driving mutations. However, these approaches can only partially capture the actual phenotype and the molecular key players of tumor growth and progression. Thus, direct investigation of the expressed proteins and activated signaling pathways provide valuable complementary information on the tumor-driving molecular characteristics of the tissue. Technological advancements in mass spectrometry-based proteomics enable the robust, rapid, and sensitive detection of thousands of proteins in minimal sample amounts, paving the way for clinical proteomics and the probing of oncogenic signaling activity. Therefore, proteomics is currently being integrated into molecular diagnostics within MTBs and holds promising potential in aiding tumor classification and identifying personalized treatment strategies. This review introduces MTBs and describes current clinical proteomics, its potential in precision oncology, and highlights the benefits of multi-omic data integration.

Evaluation of the antitumor effect of neoantigen peptide vaccines derived from the translatome of lung cancer

Emerging evidence suggests that tumor-specific neoantigens are ideal targets for cancer immunotherapy. However, how to predict tumor neoantigens based on translatome data remains obscure. Through the extraction of ribosome-nascent chain complexes (RNCs) from LLC cells, followed by RNC-mRNA extraction, RNC-mRNA sequencing, and comprehensive bioinformatic analysis, we successfully identified proteins undergoing translatome and exhibiting mutations in the cells. Subsequently, novel antigens identification was analyzed by the interaction between their high affinity and the Major Histocompatibility Complex (MHC). Neoantigens immunogenicity was analyzed by enzyme-linked immunospot assay (ELISpot). Finally, in vivo experiments in mice were conducted to evaluate the antitumor effects of translatome-derived neoantigen peptides on lung cancer. The results showed that ten neoantigen peptides were identified and synthesized by translatome data from LLC cells; 8 out of the 10 neoantigens had strong immunogenicity. The neoantigen peptide vaccine group exhibited significant tumor growth inhibition effect. In conclusion, neoantigen peptide vaccine derived from the translatome of lung cancer exhibited significant tumor growth inhibition effect.

Proteogenomic characterization of difficult-to-treat breast cancer with tumor cells enriched through laser microdissection

BACKGROUND

Breast cancer (BC) is the most commonly diagnosed cancer and the leading cause of cancer death among women globally. Despite advances, there is considerable variation in clinical outcomes for patients with non-luminal A tumors, classified as difficult-to-treat breast cancers (DTBC). This study aims to delineate the proteogenomic landscape of DTBC tumors compared to luminal A (LumA) tumors.

METHODS

We retrospectively collected a total of 117 untreated primary breast tumor specimens, focusing on DTBC subtypes. Breast tumors were processed by laser microdissection (LMD) to enrich tumor cells. DNA, RNA, and protein were simultaneously extracted from each tumor preparation, followed by whole genome sequencing, paired-end RNA sequencing, global proteomics and phosphoproteomics. Differential feature analysis, pathway analysis and survival analysis were performed to better understand DTBC and investigate biomarkers.

RESULTS

We observed distinct variations in gene mutations, structural variations, and chromosomal alterations between DTBC and LumA breast tumors. DTBC tumors predominantly had more mutations in TP53, PLXNB3, Zinc finger genes, and fewer mutations in SDC2, CDH1, PIK3CA, SVIL, and PTEN. Notably, Cytoband 1q21, which contains numerous cell proliferation-related genes, was significantly amplified in the DTBC tumors. LMD successfully minimized stromal components and increased RNA-protein concordance, as evidenced by stromal score comparisons and proteomic analysis. Distinct DTBC and LumA-enriched clusters were observed by proteomic and phosphoproteomic clustering analysis, some with survival differences. Phosphoproteomics identified two distinct phosphoproteomic profiles for high relapse-risk and low relapse-risk basal-like tumors, involving several genes known to be associated with breast cancer oncogenesis and progression, including KIAA1522, DCK, FOXO3, MYO9B, ARID1A, EPRS, ZC3HAV1, and RBM14. Lastly, an integrated pathway analysis of multi-omics data highlighted a robust enrichment of proliferation pathways in DTBC tumors.

CONCLUSIONS

This study provides an integrated proteogenomic characterization of DTBC vs LumA with tumor cells enriched through laser microdissection. We identified many common features of DTBC tumors and the phosphopeptides that could serve as potential biomarkers for high/low relapse-risk basal-like BC and possibly guide treatment selections.

Plant genome information facilitates plant functional genomics

MAIN CONCLUSION

In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

pXg: Comprehensive Identification of Noncanonical MHC-I-Associated Peptides From De Novo Peptide Sequencing Using RNA-Seq Reads

Discovering noncanonical peptides has been a common application of proteogenomics. Recent studies suggest that certain noncanonical peptides, known as noncanonical major histocompatibility complex-I (MHC-I)-associated peptides (ncMAPs), that bind to MHC-I may make good immunotherapeutic targets. De novo peptide sequencing is a great way to find ncMAPs since it can detect peptide sequences from their tandem mass spectra without using any sequence databases. However, this strategy has not been widely applied for ncMAP identification because there is not a good way to estimate its false-positive rates. In order to completely and accurately identify immunopeptides using de novo peptide sequencing, we describe a unique pipeline called proteomics X genomics. In contrast to current pipelines, it makes use of genomic data, RNA-Seq abundance and sequencing quality, in addition to proteomic features to increase the sensitivity and specificity of peptide identification. We show that the peptide-spectrum match quality and genetic traits have a clear relationship, showing that they can be utilized to evaluate peptide-spectrum matches. From 10 samples, we found 24,449 canonical MHC-I-associated peptides and 956 ncMAPs by using a target-decoy competition. Three hundred eighty-seven ncMAPs and 1611 canonical MHC-I-associated peptides were new identifications that had not yet been published. We discovered 11 ncMAPs produced from a squirrel monkey retrovirus in human cell lines in addition to the two ncMAPs originating from a complementarity determining region 3 in an antibody thanks to the unrestricted search space assumed by de novo sequencing. These entirely new identifications show that proteomics X genomics can make the most of de novo peptide sequencing's advantages and its potential use in the search for new immunotherapeutic targets.

moPepGen: Rapid and Comprehensive Proteoform Identification

Gene expression is a multi-step transformation of biological information from its storage form (DNA) into functional forms (protein and some RNAs). Regulatory activities at each step of this transformation multiply a single gene into a myriad of proteoforms. Proteogenomics is the study of how genomic and transcriptomic variation creates this proteoform diversity, and is limited by the challenges of modeling the complexities of gene-expression. We therefore created moPepGen, a graph-based algorithm that comprehensively enumerates proteoforms in linear time. moPepGen works with multiple technologies, in multiple species and on all types of genetic and transcriptomic data. In human cancer proteomes, it detects and quantifies previously unobserved noncanonical peptides arising from germline and somatic genomic variants, noncoding open reading frames, RNA fusions and RNA circularization. By enabling efficient identification and quantitation of previously hidden proteins in both existing and new proteomic data, moPepGen facilitates all proteogenomics applications. It is available at: https://github.com/uclahs-cds/package-moPepGen.

Immunopeptidomics-based identification of naturally presented non-canonical circRNA-derived peptides

Circular RNAs (circRNAs) are covalently closed non-coding RNAs lacking the 5' cap and the poly-A tail. Nevertheless, it has been demonstrated that certain circRNAs can undergo active translation. Therefore, aberrantly expressed circRNAs in human cancers could be an unexplored source of tumor-specific antigens, potentially mediating anti-tumor T cell responses. This study presents an immunopeptidomics workflow with a specific focus on generating a circRNA-specific protein fasta reference. The main goal of this workflow is to streamline the process of identifying and validating human leukocyte antigen (HLA) bound peptides potentially originating from circRNAs. We increase the analytical stringency of our workflow by retaining peptides identified independently by two mass spectrometry search engines and/or by applying a group-specific FDR for canonical-derived and circRNA-derived peptides. A subset of circRNA-derived peptides specifically encoded by the region spanning the back-splice junction (BSJ) are validated with targeted MS, and with direct Sanger sequencing of the respective source transcripts. Our workflow identifies 54 unique BSJ-spanning circRNA-derived peptides in the immunopeptidome of melanoma and lung cancer samples. Our approach enlarges the catalog of source proteins that can be explored for immunotherapy.

Analysis and visualization of quantitative proteomics data using FragPipe-Analyst

The FragPipe computational proteomics platform is gaining widespread popularity among the proteomics research community because of its fast processing speed and user-friendly graphical interface. Although FragPipe produces well-formatted output tables that are ready for analysis, there is still a need for an easy-to-use and user-friendly downstream statistical analysis and visualization tool. FragPipe-Analyst addresses this need by providing an R shiny web server to assist FragPipe users in conducting downstream analyses of the resulting quantitative proteomics data. It supports major quantification workflows including label-free quantification, tandem mass tags, and data-independent acquisition. FragPipe-Analyst offers a range of useful functionalities, such as various missing value imputation options, data quality control, unsupervised clustering, differential expression (DE) analysis using Limma, and gene ontology and pathway enrichment analysis using Enrichr. To support advanced analysis and customized visualizations, we also developed FragPipeAnalystR, an R package encompassing all FragPipe-Analyst functionalities that is extended to support site-specific analysis of post-translational modifications (PTMs). FragPipe-Analyst and FragPipeAnalystR are both open-source and freely available.

AttnPep: A Self-Attention-Based Deep Learning Method for Peptide Identification in Shotgun Proteomics

In shotgun proteomics, the proteome search engine analyzes mass spectra obtained by experiments, and then a peptide-spectra match (PSM) is reported for each spectrum. However, most of the PSMs identified are incorrect, and therefore various postprocessing software have been developed for reranking the peptide identifications. Yet these methods suffer from issues such as dependency on distribution, reliance on shallow models, and limited effectiveness. In this work, we propose AttnPep, a deep learning model for rescoring PSM scores that utilizes the Self-Attention module. This module helps the neural network focus on features relevant to the classification of PSMs and ignore irrelevant features. This allows AttnPep to analyze the output of different search engines and improve PSM discrimination accuracy. We considered a PSM to be correct if it achieves a q-value <0.01 and compared AttnPep with existing mainstream software PeptideProphet, Percolator, and proteoTorch. The results indicated that AttnPep found an average increase in correct PSMs of 9.29% relative to the other methods. Additionally, AttnPep was able to better distinguish between correct and incorrect PSMs and found more synthetic peptides in the complex SWATH data set.

A Massive Proteogenomic Screen Identifies Thousands of Novel Peptides From the Human "Dark" Proteome

Although the human gene annotation has been continuously improved over the past 2 decades, numerous studies demonstrated the existence of a "dark proteome", consisting of proteins that were critical for biological processes but not included in widely used gene catalogs. The Genotype-Tissue Expression project generated more than 15,000 RNA-seq datasets from multiple tissues, which modeled 30 million transcripts in the human genome. To provide a resource of high-confidence novel proteins from the dark proteome, we screened 50,000 mass spectrometry runs from over 900 projects to identify proteins translated from the Genotype-Tissue Expression transcript model with proteomic support. We also integrated 3.8 million common genetic variants from the gnomAD database to improve peptide identification. As a result, we identified 170,529 novel peptides with proteomic evidence, of which 6048 passed the strictest standard we defined and were supported by PepQuery. We provided a user-friendly website (https://ncorf.genes.fun/) for researchers to check the evidence of novel peptides from their studies. The findings will improve our understanding of coding genes and facilitate genomic data interpretation in biomedical research.

Proteomic insights into mental health status: plasma markers in young adults

Global emphasis on enhancing prevention and treatment strategies necessitates an increased understanding of the biological mechanisms of psychopathology. Plasma proteomics is a powerful tool that has been applied in the context of specific mental disorders for biomarker identification. The p-factor, also known as the "general psychopathology factor", is a concept in psychopathology suggesting that there is a common underlying factor that contributes to the development of various forms of mental disorders. It has been proposed that the p-factor can be used to understand the overall mental health status of an individual. Here, we aimed to discover plasma proteins associated with the p-factor in 775 young adults in the FinnTwin12 cohort. Using liquid chromatography-tandem mass spectrometry, 13 proteins with a significant connection with the p-factor were identified, 8 of which were linked to epidermal growth factor receptor (EGFR) signaling. This exploratory study provides new insight into biological alterations associated with mental health status in young adults.

SpliceProt 2.0: A Sequence Repository of Human, Mouse, and Rat Proteoforms

SpliceProt 2.0 is a public proteogenomics database that aims to list the sequence of known proteins and potential new proteoforms in human, mouse, and rat proteomes. This updated repository provides an even broader range of computationally translated proteins and serves, for example, to aid with proteomic validation of splice variants absent from the reference UniProtKB/SwissProt database. We demonstrate the value of SpliceProt 2.0 to predict orthologous proteins between humans and murines based on transcript reconstruction, sequence annotation and detection at the transcriptome and proteome levels. In this release, the annotation data used in the reconstruction of transcripts based on the methodology of ternary matrices were acquired from new databases such as Ensembl, UniProt, and APPRIS. Another innovation implemented in the pipeline is the exclusion of transcripts predicted to be susceptible to degradation through the NMD pathway. Taken together, our repository and its applications represent a valuable resource for the proteogenomics community.

Intracellular and Extracellular Peptidomes of the Model Plant, Physcomitrium patens

Here, we report our approach to peptidomic analysis of the plant model Physcomitrium patens. Intracellular and extracellular peptides were extracted under conditions preventing proteolytic digestion by endogenous proteases. The extracts were fractionated on size exclusion columns to isolate intracellular peptides and on reversed-phase cartridges to isolate extracellular peptides, with the isolated peptides subjected to LC-MS/MS analysis. Mass spectrometry data were analyzed for the presence of peptides derived from the known proteins or microproteins encoded by small open reading frames (<100 aa, smORFs) predicted in the moss genome. Experimental details are provided for each step.

Tumour circular RNAs elicit anti-tumour immunity by encoding cryptic peptides

Emerging data have shown that previously defined noncoding genomes might encode peptides that bind human leukocyte antigen (HLA) as cryptic antigens to stimulate adaptive immunity1,2. However, the significance and mechanisms of action of cryptic antigens in anti-tumour immunity remain unclear. Here mass spectrometry of the HLA class I (HLA-I) peptidome coupled with ribosome sequencing of human breast cancer samples identified HLA-I-binding cryptic antigenic peptides that were noncanonically translated by a tumour-specific circular RNA (circRNA): circFAM53B. The cryptic peptides efficiently primed naive CD4+ and CD8+ T cells in an antigen-specific manner and induced anti-tumour immunity. Clinically, the expression of circFAM53B and its encoded peptides was associated with substantial infiltration of antigen-specific CD8+ T cells and better survival in patients with breast cancer and patients with melanoma. Mechanistically, circFAM53B-encoded peptides had strong binding affinity to both HLA-I and HLA-II molecules. In vivo, administration of vaccines consisting of tumour-specific circRNA or its encoded peptides in mice bearing breast cancer tumours or melanoma induced enhanced infiltration of tumour-antigen-specific cytotoxic T cells, which led to effective tumour control. Overall, our findings reveal that noncanonical translation of circRNAs can drive efficient anti-tumour immunity, which suggests that vaccination exploiting tumour-specific circRNAs may serve as an immunotherapeutic strategy against malignant tumours.

Exploring the Alternative Proteome with OpenProt and Mass Spectrometry

Proteogenomics has revealed the translation of unannotated open reading frames (ORFs) present in mRNAs and in noncoding RNAs (ncRNAs). OpenProt annotates all ORFs with a minimum of 30 codons in the transcriptome of several species and displays many functional features associated with the corresponding proteins. Two types of proteins are annotated: reference or canonical proteins which are proteins already annotated in UniProt, RefSeq, or Ensembl and noncanonical proteins. Noncanonical proteins form two groups: predicted novel isoforms that display a significant level of homology with a reference protein and alternative proteins that are new proteins with no significant homology to known proteins. This chapter describes how to check whether a gene and/or transcript contains multiple open reading frames and how to use OpenProt databases for the detection of alternative proteins and novel isoforms by mass spectrometry-based proteomics.

Bacterial tolerance and detoxification of cyanide, arsenic and heavy metals: Holistic approaches applied to bioremediation of industrial complex wastes

Cyanide is a highly toxic compound that is found in wastewaters generated from different industrial activities, such as mining or jewellery. These residues usually contain high concentrations of other toxic pollutants like arsenic and heavy metals that may form different complexes with cyanide. To develop bioremediation strategies, it is necessary to know the metabolic processes involved in the tolerance and detoxification of these pollutants, but most of the current studies are focused on the characterization of the microbial responses to each one of these environmental hazards individually, and the effect of co-contaminated wastes on microbial metabolism has been hardly addressed. This work summarizes the main strategies developed by bacteria to alleviate the effects of cyanide, arsenic and heavy metals, analysing interactions among these toxic chemicals. Additionally, it is discussed the role of systems biology and synthetic biology as tools for the development of bioremediation strategies of complex industrial wastes and co-contaminated sites, emphasizing the importance and progress derived from meta-omic studies.

Integrative Proteogenomics for Differential Expression and Splicing Variation in a DM1 Mouse Model

Dysregulated mRNA splicing is involved in the pathogenesis of many diseases including cancer, neurodegenerative diseases, and muscular dystrophies such as myotonic dystrophy type 1 (DM1). Comprehensive assessment of dysregulated splicing on the transcriptome and proteome level has been methodologically challenging, and thus investigations have often been targeting only few genes. Here, we performed a large-scale coordinated transcriptomic and proteomic analysis to characterize a DM1 mouse model (HSALR) in comparison to wild type. Our integrative proteogenomics approach comprised gene- and splicing-level assessments for mRNAs and proteins. It recapitulated many known instances of aberrant mRNA splicing in DM1 and identified new ones. It enabled the design and targeting of splicing-specific peptides and confirmed the translation of known instances of aberrantly spliced disease-related genes (e.g., Atp2a1, Bin1, Ryr1), complemented by novel findings (Flnc and Ywhae). Comparative analysis of large-scale mRNA and protein expression data showed quantitative agreement of differentially expressed genes and splicing patterns between disease and wild type. We hence propose this work as a suitable blueprint for a robust and scalable integrative proteogenomic strategy geared toward advancing our understanding of splicing-based disorders. With such a strategy, splicing-based biomarker candidates emerge as an attractive and accessible option, as they can be efficiently asserted on the mRNA and protein level in coordinated fashion.

Identification of proteoforms of short open reading frame-encoded peptides in Blautia producta under different cultivation conditions

IMPORTANCE

The identification of short open reading frame-encoded peptides (SEP) and different proteoforms in single cultures of gut microbes offers new insights into a largely neglected part of the microbial proteome landscape. This is of particular importance as SEP provide various predicted functions, such as acting as antimicrobial peptides, maintaining cell homeostasis under stress conditions, or even contributing to the virulence pattern. They are, thus, taking a poorly understood role in structure and function of microbial networks in the human body. A better understanding of SEP in the context of human health requires a precise understanding of the abundance of SEP both in commensal microbes as well as pathogens. For the gut beneficial B. producta, we demonstrate the importance of specific environmental conditions for biosynthesis of SEP expanding previous findings about their role in microbial interactions.

Multi-tissue proteogenomic analysis for mechanistic toxicology studies in non-model species

New approach methodologies (NAM), including omics and in vitro approaches, are contributing to the implementation of 3R (reduction, refinement and replacement) strategies in regulatory science and risk assessment. In this study, we present an integrative transcriptomics and proteomics analysis workflow for the validation and revision of complex fish genomes and demonstrate how proteogenomics expression matrices can be used to support multi-level omics data integration in non-model species in vivo and in vitro. Using Atlantic salmon as an example, we constructed proteogenomic databases from publicly available transcriptomic data and in-house generated RNA-Seq and LC-MS/MS data. Our analysis identified ∼80,000 peptides, providing direct evidence of translation for over 40,000 RefSeq structures. The data also highlighted 183 co-located peptide groups that supported a single transcript each, and in each case, either corrected a previous annotation, supported Ensembl annotations not present in RefSeq, or identified novel previously unannotated genes. Proteogenomics data-derived expression matrices revealed distinct profiles for the different tissue types analyzed. Focusing on proteins involved in defense against xenobiotics, we detected distinct expression patterns across different salmon tissues and observed homology in the expression of chemical defense proteins between in vivo and in vitro liver systems. Our study demonstrates the potential of proteogenomic analyses in extending our understanding of complex fish genomes and provides an advanced bioinformatic toolkit to support the further development of NAMs and their application in regulatory science and (eco)toxicological studies of non-model species.

Proteogenomics-based functional genome research: approaches, applications, and perspectives in plants

Proteogenomics (PG) integrates the proteome with the genome and transcriptome to refine gene models and annotation. Coupled with single-cell (SC) assays, PG effectively distinguishes heterogeneity among cell groups. Affiliating spatial information to PG reveals the high-resolution circuitry within SC atlases. Additionally, PG can investigate dynamic changes in protein-coding genes in plants across growth and development as well as stress and external stimulation, significantly contributing to the functional genome. Here we summarize existing PG research in plants and introduce the technical features of various methods. Combining PG with other omics, such as metabolomics and peptidomics, can offer even deeper insights into gene functions. We argue that the application of PG will represent an important font of foundational knowledge for plants.

Small proteins in bacteria - Big challenges in prediction and identification

Proteins with up to 100 amino acids have been largely overlooked due to the challenges associated with predicting and identifying them using traditional methods. Recent advances in bioinformatics and machine learning, DNA sequencing, RNA and Ribo-seq technologies, and mass spectrometry (MS) have greatly facilitated the detection and characterisation of these elusive proteins in recent years. This has revealed their crucial role in various cellular processes including regulation, signalling and transport, as toxins and as folding helpers for protein complexes. Consequently, the systematic identification and characterisation of these proteins in bacteria have emerged as a prominent field of interest within the microbial research community. This review provides an overview of different strategies for predicting and identifying these proteins on a large scale, leveraging the power of these advanced technologies. Furthermore, the review offers insights into the future developments that may be expected in this field.

Biological factors and statistical limitations prevent detection of most noncanonical proteins by mass spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry (MS) experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here, we leveraged recent advances in ribosome profiling and MS to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly expressed to be detected by shotgun MS at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for 4 noncanonical proteins in MS data, which were also supported by evolution and translation data. These results illustrate the power of MS to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly expressed proteins.

Recent advances in immunopeptidomic-based tumor neoantigen discovery

The role of aberrantly expressed proteins in tumors in driving immune-mediated control of cancer has been well documented for more than five decades. Today, we know that both aberrantly expressed normal proteins as well as mutant proteins (neoantigens) can function as tumor antigens in both humans and mice. Next-generation sequencing (NGS) and high-resolution mass spectrometry (MS) technologies have made significant advances since the early 2010s, enabling detection of rare but clinically relevant neoantigens recognized by T cells. MS profiling of tumor-specific immunopeptidomes remains the most direct method to identify mutant peptides bound to cellular MHC. However, the need for use of large numbers of cells or significant amounts of tumor tissue to achieve neoantigen detection has historically limited the application of MS. Newer, more sensitive MS technologies have recently demonstrated the capacities to detect neoantigens from fewer cells. Here, we highlight recent advancements in immunopeptidomics-based characterization of tumor-specific neoantigens. Various tumor antigen categories and neoantigen identification approaches are also discussed. Furthermore, we summarize recent reports that achieved successful tumor neoantigen detection by MS using a variety of starting materials, MS acquisition modes, and novel ion mobility devices.

Small proteins in Gram-positive bacteria

Small proteins comprising less than 100 amino acids have been often ignored in bacterial genome annotations. About 10 years ago, focused efforts started to investigate whole peptidomes, which resulted in the discovery of a multitude of small proteins, but only a number of them have been characterized in detail. Generally, small proteins can be either membrane or cytosolic proteins. The latter interact with larger proteins, RNA or even metal ions. Here, we summarize our current knowledge on small proteins from Gram-positive bacteria with a special emphasis on the model organism Bacillus subtilis. Our examples include membrane-bound toxins of type I toxin-antitoxin systems, proteins that block the assembly of higher order structures, regulate sporulation or modulate the RNA degradosome. We do not consider antimicrobial peptides. Furthermore, we present methods for the identification and investigation of small proteins.

Biological Factors and Statistical Limitations Prevent Detection of Most Noncanonical Proteins by Mass Spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here we leveraged recent advances in ribosome profiling and mass spectrometry to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly-expressed to be detected by shotgun mass spectrometry at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for four noncanonical proteins in mass spectrometry data, which were also supported by evolution and translation data. These results illustrate the power of mass spectrometry to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly-expressed proteins.

Proteogenomics in Cancer: Then and Now

For years, the paths of sequencing technologies and mass spectrometry have occurred in isolation, with each developing its own unique culture and expertise. These two technologies are crucial for inspecting complementary aspects of the molecular phenotype across the central dogma. Integrative multiomics strives to bridge the analysis gap among different fields to complete more comprehensive mechanisms of life events and diseases. Proteogenomics is one integrated multiomics field. Here in this review, we mainly summarize and discuss three aspects: workflow of proteogenomics, proteogenomics applications in cancer research, and the SWOT (Strengths, Weaknesses, Opportunities, Threats) analysis of proteogenomics in cancer research. In conclusion, proteogenomics has a promising future as it clarifies the functional consequences of many unannotated genomic abnormalities or noncanonical variants and identifies driver genes and novel therapeutic targets across cancers, which would substantially accelerate the development of precision oncology.

Retention Time and Fragmentation Predictors Increase Confidence in Identification of Common Variant Peptides

Precision medicine focuses on adapting care to the individual profile of patients, for example, accounting for their unique genetic makeup. Being able to account for the effect of genetic variation on the proteome holds great promise toward this goal. However, identifying the protein products of genetic variation using mass spectrometry has proven very challenging. Here we show that the identification of variant peptides can be improved by the integration of retention time and fragmentation predictors into a unified proteogenomic pipeline. By combining these intrinsic peptide characteristics using the search-engine post-processor Percolator, we demonstrate improved discrimination power between correct and incorrect peptide-spectrum matches. Our results demonstrate that the drop in performance that is induced when expanding a protein sequence database can be compensated, hence enabling efficient identification of genetic variation products in proteomics data. We anticipate that this enhancement of proteogenomic pipelines can provide a more refined picture of the unique proteome of patients and thereby contribute to improving patient care.

The Breast Cancer Proteome and Precision Oncology

The goal of precision oncology is to translate the molecular features of cancer into predictive and prognostic tests that can be used to individualize treatment leading to improved outcomes and decreased toxicity. Success for this strategy in breast cancer is exemplified by efficacy of trastuzumab in tumors overexpressing ERBB2 and endocrine therapy for tumors that are estrogen receptor positive. However, other effective treatments, including chemotherapy, immune checkpoint inhibitors, and CDK4/6 inhibitors are not associated with strong predictive biomarkers. Proteomics promises another tier of information that, when added to genomic and transcriptomic features (proteogenomics), may create new opportunities to improve both treatment precision and therapeutic hypotheses. Here, we review both mass spectrometry-based and antibody-dependent proteomics as complementary approaches. We highlight how these methods have contributed toward a more complete understanding of breast cancer and describe the potential to guide diagnosis and treatment more accurately.

The Integration of Proteogenomics and Ribosome Profiling Circumvents Key Limitations to Increase the Coverage and Confidence of Novel Microproteins

There has been a dramatic increase in the identification of non-conical translation and a significant expansion of the protein-coding genome and proteome. Among the strategies used to identify novel small ORFs (smORFs), Ribosome profiling (Ribo-Seq) is the gold standard for the annotation of novel coding sequences by reporting on smORF translation. In Ribo-Seq, ribosome-protected footprints (RPFs) that map to multiple sites in the genome are computationally removed since they cannot unambiguously be assigned to a specific genomic location, or to a specific transcript in the case of multiple isoforms. Furthermore, RPFs necessarily result in short (25-34 nucleotides) reads, increasing the chance of ambiguous and multi-mapping alignments, such that smORFs that reside in these regions cannot be identified by Ribo-Seq. Here, we show that the inclusion of proteogenomics to create a Ribosome Profiling and Proteogenomics Pipeline (RP3) bypasses this limitation to identify a group of microprotein-encoding smORFs that are missed by current Ribo-Seq pipelines. Moreover, we show that the microproteins identified by RP3 have different sequence compositions from the ones identified by Ribo-Seq-only pipelines, which can affect proteomics identification. In aggregate, the development of RP3 maximizes the detection and confidence of protein-encoding smORFs and microproteins.

Metaproteogenomic Profile of a Mesopelagic Adenylylsulfate Reductase: Course-Based Discovery Using the Ocean Protein Portal

Adenylylsulfate reductase (Apr) is a flavoprotein with a dissimilatory sulfate reductase function. Its ability to catalyze the reverse reaction in sulfur oxidizers has propelled a complex phylogenetic history of transfers with sulfate reducers and made this enzyme an important protein in ocean sulfur cycling. As part of a graduate course, we analyzed metaproteomic data from the Ocean Protein Portal and observed evidence of Apr alpha (AprA) and beta (AprB) subunits in the Central Pacific Ocean. The protein was originally taxonomically attributed toChlorobium tepidum TLS, a green sulfur bacterium. However, our phylogenomic and oceanographic contextual analysis contradicted this label, instead showing that this protein is consistent with the genomic material from the newly discovered Candidatus Lambdaproteobacteriaclass, implying that the ecological role of this lineage in oxygen minimum twilight zones is underappreciated. This study illustrates how metaproteogenomic analysis can contribute to more accurate metagenomic/proteomic annotations and comprehensive ocean biogeochemical processes conducive to course-based research experiences.

Multi-omic stratification of the missense variant cysteinome

Cancer genomes are rife with genetic variants; one key outcome of this variation is gain-ofcysteine, which is the most frequently acquired amino acid due to missense variants in COSMIC. Acquired cysteines are both driver mutations and sites targeted by precision therapies. However, despite their ubiquity, nearly all acquired cysteines remain uncharacterized. Here, we pair cysteine chemoproteomics-a technique that enables proteome-wide pinpointing of functional, redox sensitive, and potentially druggable residues-with genomics to reveal the hidden landscape of cysteine acquisition. For both cancer and healthy genomes, we find that cysteine acquisition is a ubiquitous consequence of genetic variation that is further elevated in the context of decreased DNA repair. Our chemoproteogenomics platform integrates chemoproteomic, whole exome, and RNA-seq data, with a customized 2-stage false discovery rate (FDR) error controlled proteomic search, further enhanced with a user-friendly FragPipe interface. Integration of CADD predictions of deleteriousness revealed marked enrichment for likely damaging variants that result in acquisition of cysteine. By deploying chemoproteogenomics across eleven cell lines, we identify 116 gain-of-cysteines, of which 10 were liganded by electrophilic druglike molecules. Reference cysteines proximal to missense variants were also found to be pervasive, 791 in total, supporting heretofore untapped opportunities for proteoform-specific chemical probe development campaigns. As chemoproteogenomics is further distinguished by sample-matched combinatorial variant databases and compatible with redox proteomics and small molecule screening, we expect widespread utility in guiding proteoform-specific biology and therapeutic discovery.

Identification of Splice Variants and Isoforms in Transcriptomics and Proteomics

Alternative splicing is pivotal to the regulation of gene expression and protein diversity in eukaryotic cells. The detection of alternative splicing events requires specific omics technologies. Although short-read RNA sequencing has successfully supported a plethora of investigations on alternative splicing, the emerging technologies of long-read RNA sequencing and top-down mass spectrometry open new opportunities to identify alternative splicing and protein isoforms with less ambiguity. Here, we summarize improvements in short-read RNA sequencing for alternative splicing analysis, including percent splicing index estimation and differential analysis. We also review the computational methods used in top-down proteomics analysis regarding proteoform identification, including the construction of databases of protein isoforms and statistical analyses of search results. While many improvements in sequencing and computational methods will result from emerging technologies, there should be future endeavors to increase the effectiveness, integration, and proteome coverage of alternative splicing events.

metaSpectraST: an unsupervised and database-independent analysis workflow for metaproteomic MS/MS data using spectrum clustering

BACKGROUND

The high diversity and complexity of the microbial community make it a formidable challenge to identify and quantify the large number of proteins expressed in the community. Conventional metaproteomics approaches largely rely on accurate identification of the MS/MS spectra to their corresponding short peptides in the digested samples, followed by protein inference and subsequent taxonomic and functional analysis of the detected proteins. These approaches are dependent on the availability of protein sequence databases derived either from sample-specific metagenomic data or from public repositories. Due to the incompleteness and imperfections of these protein sequence databases, and the preponderance of homologous proteins expressed by different bacterial species in the community, this computational process of peptide identification and protein inference is challenging and error-prone, which hinders the comparison of metaproteomes across multiple samples.

RESULTS

We developed metaSpectraST, an unsupervised and database-independent metaproteomics workflow, which quantitatively profiles and compares metaproteomics samples by clustering experimentally observed MS/MS spectra based on their spectral similarity. We applied metaSpectraST to fecal samples collected from littermates of two different mother mice right after weaning. Quantitative proteome profiles of the microbial communities of different mice were obtained without any peptide-spectrum identification and used to evaluate the overall similarity between samples and highlight any differentiating markers. Compared to the conventional database-dependent metaproteomics analysis, metaSpectraST is more successful in classifying the samples and detecting the subtle microbiome changes of mouse gut microbiomes post-weaning. metaSpectraST could also be used as a tool to select the suitable biological replicates from samples with wide inter-individual variation.

CONCLUSIONS

metaSpectraST enables rapid profiling of metaproteomic samples quantitatively, without the need for constructing the protein sequence database or identification of the MS/MS spectra. It maximally preserves information contained in the experimental MS/MS spectra by clustering all of them first and thus is able to better profile the complex microbial communities and highlight their functional changes, as compared with conventional approaches. tag the videobyte in this section as ESM4 Video Abstract.

Precision Medicine: Disease Subtyping and Tailored Treatment

The genomics-based concept of precision medicine began to emerge following the completion of the Human Genome Project. In contrast to evidence-based medicine, precision medicine will allow doctors and scientists to tailor the treatment of different subpopulations of patients who differ in their susceptibility to specific diseases or responsiveness to specific therapies. The current precision medicine model was proposed to precisely classify patients into subgroups sharing a common biological basis of diseases for more effective tailored treatment to achieve improved outcomes. Precision medicine has become a term that symbolizes the new age of medicine. In this review, we examine the history, development, and future perspective of precision medicine. We also discuss the concepts, principles, tools, and applications of precision medicine and related fields. In our view, for precision medicine to work, two essential objectives need to be achieved. First, diseases need to be classified into various subtypes. Second, targeted therapies must be available for each specific disease subtype. Therefore, we focused this review on the progress in meeting these two objectives.

MSBooster: improving peptide identification rates using deep learning-based features

Peptide identification in liquid chromatography-tandem mass spectrometry (LC-MS/MS) experiments relies on computational algorithms for matching acquired MS/MS spectra against sequences of candidate peptides using database search tools, such as MSFragger. Here, we present a new tool, MSBooster, for rescoring peptide-to-spectrum matches using additional features incorporating deep learning-based predictions of peptide properties, such as LC retention time, ion mobility, and MS/MS spectra. We demonstrate the utility of MSBooster, in tandem with MSFragger and Percolator, in several different workflows, including nonspecific searches (immunopeptidomics), direct identification of peptides from data independent acquisition data, single-cell proteomics, and data generated on an ion mobility separation-enabled timsTOF MS platform. MSBooster is fast, robust, and fully integrated into the widely used FragPipe computational platform.