Aberrant splicing isoforms detected by full-length transcriptome sequencing as transcripts of potential neoantigens in non-small cell lung cancer

BACKGROUND

Long-read sequencing of full-length cDNAs enables the detection of structures of aberrant splicing isoforms in cancer cells. These isoforms are occasionally translated, presented by HLA molecules, and recognized as neoantigens. This study used a long-read sequencer (MinION) to construct a comprehensive catalog of aberrant splicing isoforms in non-small-cell lung cancers, by which novel isoforms and potential neoantigens are identified.

RESULTS

Full-length cDNA sequencing is performed using 22 cell lines, and a total of 2021 novel splicing isoforms are identified. The protein expression of some of these isoforms is then validated by proteome analysis. Ablations of a nonsense-mediated mRNA decay (NMD) factor, UPF1, and a splicing factor, SF3B1, are found to increase the proportion of aberrant transcripts. NetMHC evaluation of the binding affinities to each type of HLA molecule reveals that some of the isoforms potentially generate neoantigen candidates. We also identify aberrant splicing isoforms in seven non-small-cell lung cancer specimens. An enzyme-linked immune absorbent spot assay indicates that approximately half the peptide candidates have the potential to activate T cell responses through their interaction with HLA molecules. Finally, we estimate the number of isoforms in The Cancer Genome Atlas (TCGA) datasets by referring to the constructed catalog and found that disruption of NMD factors is significantly correlated with the number of splicing isoforms found in the TCGA-Lung Adenocarcinoma data collection.

CONCLUSIONS

Our results indicate that long-read sequencing of full-length cDNAs is essential for the precise identification of aberrant transcript structures in cancer cells.

The jPOST Repository as a Public Data Repository for Shotgun Proteomics

In recent years, mass spectrometry-based proteomics approach has made significant progress and the number of datasets related to various proteomics projects has increased worldwide. To promote the sharing and reuse of promising datasets, it is important to build an appropriate, high-quality public data repository. For this purpose, several repositories have already been created. The jPOST repository that we developed in 2016 has successfully implemented several unique features, such as fast file upload, flexible file management, and an easy-to-use interface. In addition, this repository is an official member of the ProteomeXchange Consortium established to facilitate standard data submission and global dissemination of mass spectrometry proteomics data. Our repository contributes to the global partnership for sharing and storing all the datasets related to various proteomics experiments.

The jPOST environment: an integrated proteomics data repository and database

Rapid progress is being made in mass spectrometry (MS)-based proteomics, yielding an increasing number of larger datasets with higher quality and higher throughput. To integrate proteomics datasets generated from various projects and institutions, we launched a project named jPOST (Japan ProteOme STandard Repository/Database, https://jpostdb.org/) in 2015. Its proteomics data repository, jPOSTrepo, began operations in 2016 and has accepted more than 10 TB of MS-based proteomics datasets in the past two years. In addition, we have developed a new proteomics database named jPOSTdb in which the published raw datasets in jPOSTrepo are reanalyzed using standardized protocol. jPOSTdb provides viewers showing the frequency of detected post-translational modifications, the co-occurrence of phosphorylation sites on a peptide and peptide sharing among proteoforms. jPOSTdb also provides basic statistical analysis tools to compare proteomics datasets.

BioHackathon 2015: Semantics of data for life sciences and reproducible research

We report on the activities of the 2015 edition of the BioHackathon, an annual event that brings together researchers and developers from around the world to develop tools and technologies that promote the reusability of biological data. We discuss issues surrounding the representation, publication, integration, mining and reuse of biological data and metadata across a wide range of biomedical data types of relevance for the life sciences, including chemistry, genotypes and phenotypes, orthology and phylogeny, proteomics, genomics, glycomics, and metabolomics. We describe our progress to address ongoing challenges to the reusability and reproducibility of research results, and identify outstanding issues that continue to impede the progress of bioinformatics research. We share our perspective on the state of the art, continued challenges, and goals for future research and development for the life sciences Semantic Web.

jPOSTrepo: an international standard data repository for proteomes

Major advancements have recently been made in mass spectrometry-based proteomics, yielding an increasing number of datasets from various proteomics projects worldwide. In order to facilitate the sharing and reuse of promising datasets, it is important to construct appropriate, high-quality public data repositories. jPOSTrepo (https://repository.jpostdb.org/) has successfully implemented several unique features, including high-speed file uploading, flexible file management and easy-to-use interfaces. This repository has been launched as a public repository containing various proteomic datasets and is available for researchers worldwide. In addition, our repository has joined the ProteomeXchange consortium, which includes the most popular public repositories such as PRIDE in Europe for MS/MS datasets and PASSEL for SRM datasets in the USA. Later MassIVE was introduced in the USA and accepted into the ProteomeXchange, as was our repository in July 2016, providing important datasets from Asia/Oceania. Accordingly, this repository thus contributes to a global alliance to share and store all datasets from a wide variety of proteomics experiments. Thus, the repository is expected to become a major repository, particularly for data collected in the Asia/Oceania region.

Identification of human basic fetoprotein as glucose-6-phosphate isomerase by using N- and C-terminal sequence tags and terminal tag database

Human basic fetoprotein (BFP), found in fetal serum and tissue extracts as well as in extracts of various cancer tissues, has long been known as a marker protein for cancers; however, the primary sequence has not yet been reported. This paper describes the identification of BFP using the N- and C-terminal amino acid sequence tags (Ac-AALTRDPQFQ and QQREARVQ, respectively) clarified by mass spectrometry-based methods, and a terminal tag database (ProteinCarta). In this study, BFP was identified as glucose-6-phosphate isomerase (G6PI_HUMAN).

Computational survey of sequence specificity for protein terminal tags covering nine organisms and its application to protein identification

In 1998, Wilkins et al. (J. Mol. Biol. 1998, 278, 599-608) reported high specificity in terminal regions (terminal tags) of 15 519 proteins from five organisms and proposed a methodology for identifying proteins by terminal tags. However, their examined sequence data were not based on complete genome sequences. Here, we examined current proteome data (217 249 entries from UniProt 2013_6 complete/reference proteome for nine organisms including human) in terms of the specificity of terminal tags and their computational annotation. One example from the results indicated that the specificity of N-terminal tags plateaued at 28% at a length of six residues for human; even when using both N- and C-terminal tags, specificity was merely 66%. In order to determine the cause of these low specificities, the annotation of proteins sharing terminal tags with other proteins was examined. The results suggested that a large majority were phylogenetically or functionally related, whereas nonrelated proteins sharing terminal tags made up less than 1% of human proteome data. On the basis of these findings, we constructed the terminal tag sequence database ProteinCarta (http://ms3d.jp/software/proteincarta/), which includes all terminal tags of proteomes from the nine organisms analyzed here, in order to confirm the specificity of terminal tags and to identify the parent protein.

A simple peak detection and label-free quantitation algorithm for chromatography-mass spectrometry

Background

Label-free quantitation of mass spectrometric data is one of the simplest and least expensive methods for differential expression profiling of proteins and metabolites. The need for high accuracy and performance computational label-free quantitation methods is still high in the biomarker and drug discovery research field. However, recent most advanced types of LC-MS generate huge amounts of analytical data with high scan speed, high accuracy and resolution, which is often impossible to interpret manually. Moreover, there are still issues to be improved for recent label-free methods, such as how to reduce false positive/negatives of the candidate peaks, how to expand scalability and how to enhance and automate data processing. AB3D (A simple label-free quantitation algorithm for Biomarker Discovery in Diagnostics and Drug discovery using LC-MS) has addressed these issues and has the capability to perform label-free quantitation using MS1 for proteomics study.

Results

We developed an algorithm called AB3D, a label free peak detection and quantitative algorithm using MS1 spectral data. To test our algorithm, practical applications of AB3D for LC-MS data sets were evaluated using 3 datasets. Comparisons were then carried out between widely used software tools such as MZmine 2, MSight, SuperHirn, OpenMS and our algorithm AB3D, using the same LC-MS datasets. All quantitative results were confirmed manually, and we found that AB3D could properly identify and quantify known peptides with fewer false positives and false negatives compared to four other existing software tools using either the standard peptide mixture or the real complex biological samples of Bartonella quintana (strain JK31). Moreover, AB3D showed the best reliability by comparing the variability between two technical replicates using a complex peptide mixture of HeLa and BSA samples. For performance, the AB3D algorithm is about 1.2 - 15 times faster than the four other existing software tools.

Conclusions

AB3D is a simple and fast algorithm for label-free quantitation using MS1 mass spectrometry data for large scale LC-MS data analysis with higher true positive and reasonable false positive rates. Furthermore, AB3D demonstrated the best reproducibility and is about 1.2- 15 times faster than those of existing 4 software tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0376-0) contains supplementary material, which is available to authorized users.

Integrated proteomics identified novel activation of dynein IC2-GR-COX-1 signaling in neurofibromatosis type I (NF1) disease model cells

Neurofibromatosis type 1 (NF1) tumor suppressor gene product, neurofibromin, functions in part as a Ras-GAP, and though its loss is implicated in the neuronal abnormality of NF1 patients, its precise cellular function remains unclear. To study the molecular mechanism of NF1 pathogenesis, we prepared NF1 gene knockdown (KD) PC12 cells, as a NF1 disease model, and analyzed their molecular (gene and protein) expression profiles with a unique integrated proteomics approach, comprising iTRAQ, 2D-DIGE, and DNA microarrays, using an integrated protein and gene expression analysis chart (iPEACH). In NF1-KD PC12 cells showing abnormal neuronal differentiation after NGF treatment, of 3198 molecules quantitatively identified and listed in iPEACH, 97 molecules continuously up- or down-regulated over time were extracted. Pathway and network analysis further revealed overrepresentation of calcium signaling and transcriptional regulation by glucocorticoid receptor (GR) in the up-regulated protein set, whereas nerve system development was overrepresented in the down-regulated protein set. The novel up-regulated network we discovered, "dynein IC2-GR-COX-1 signaling," was then examined in NF1-KD cells. Validation studies confirmed that NF1 knockdown induces altered splicing and phosphorylation patterns of dynein IC2 isomers, up-regulation and accumulation of nuclear GR, and increased COX-1 expression in NGF-treated cells. Moreover, the neurite retraction phenotype observed in NF1-KD cells was significantly recovered by knockdown of the dynein IC2-C isoform and COX-1. In addition, dynein IC2 siRNA significantly inhibited nuclear translocation and accumulation of GR and up-regulation of COX-1 expression. These results suggest that dynein IC2 up-regulates GR nuclear translocation and accumulation, and subsequently causes increased COX-1 expression, in this NF1 disease model. Our integrated proteomics strategy, which combines multiple approaches, demonstrates that NF1-related neural abnormalities are, in part, caused by up-regulation of dynein IC2-GR-COX-1 signaling, which may be a novel therapeutic target for NF1.

KEGG OC: a large-scale automatic construction of taxonomy-based ortholog clusters

The identification of orthologous genes in an increasing number of fully sequenced genomes is a challenging issue in recent genome science. Here we present KEGG OC (http://www.genome.jp/tools/oc/), a novel database of ortholog clusters (OCs). The current version of KEGG OC contains 1 176 030 OCs, obtained by clustering 8 357 175 genes in 2112 complete genomes (153 eukaryotes, 1830 bacteria and 129 archaea). The OCs were constructed by applying the quasi-clique-based clustering method to all possible protein coding genes in all complete genomes, based on their amino acid sequence similarities. It is computationally efficient to calculate OCs, which enables to regularly update the contents. KEGG OC has the following two features: (i) It consists of all complete genomes of a wide variety of organisms from three domains of life, and the number of organisms is the largest among the existing databases; and (ii) It is compatible with the KEGG database by sharing the same sets of genes and identifiers, which leads to seamless integration of OCs with useful components in KEGG such as biological pathways, pathway modules, functional hierarchy, diseases and drugs. The KEGG OC resources are accessible via OC Viewer that provides an interactive visualization of OCs at different taxonomic levels.

CIPRO 2.5: Ciona intestinalis protein database, a unique integrated repository of large-scale omics data, bioinformatic analyses and curated annotation, with user rating and reviewing functionality

The Ciona intestinalis protein database (CIPRO) is an integrated protein database for the tunicate species C. intestinalis. The database is unique in two respects: first, because of its phylogenetic position, Ciona is suitable model for understanding vertebrate evolution; and second, the database includes original large-scale transcriptomic and proteomic data. Ciona intestinalis has also been a favorite of developmental biologists. Therefore, large amounts of data exist on its development and morphology, along with a recent genome sequence and gene expression data. The CIPRO database is aimed at collecting those published data as well as providing unique information from unpublished experimental data, such as 3D expression profiling, 2D-PAGE and mass spectrometry-based large-scale analyses at various developmental stages, curated annotation data and various bioinformatic data, to facilitate research in diverse areas, including developmental, comparative and evolutionary biology. For medical and evolutionary research, homologs in humans and major model organisms are intentionally included. The current database is based on a recently developed KH model containing 36 034 unique sequences, but for higher usability it covers 89 683 all known and predicted proteins from all gene models for this species. Of these sequences, more than 10 000 proteins have been manually annotated. Furthermore, to establish a community-supported protein database, these annotations are open to evaluation by users through the CIPRO website. CIPRO 2.5 is freely accessible at http://cipro.ibio.jp/2.5.

ODB: a database for operon organizations, 2011 update

ODB (Operon DataBase) aims to collect data of all known and conserved operons in completely sequenced genomes. Three newly updated features of this database have been added as follows: (i) Data from included operons were updated. The genome-wide analysis of transcription and transcriptional units has become popular recently and ODB successfully integrates these high-throughput operon data, including genome-wide transcriptional units of five prokaryotes and two eukaryotes. The current version of our database contains information from about 10 000 known operons in more than 50 genomes, and more than 400 000 conserved operons obtained from more than 1000 bacterial genomes. (ii) ODB proposes the idea of reference operons as a new operon prediction tool. A reference operon, a set of possible orthologous genes that organize operons, is defined by clustering all known operons. A large number of known operons, including the recently added genome-wide analysis of operons, allowed us to define more reliable reference operons. (iii) ODB also provides new graphical interfaces. One is for comparative analyses of operon structures in multiple genomes. The other is for visualization of possible operons in multiple genomes obtained from the reference operons. The 2011 updated version of ODB is now available at http://operondb.jp/.

Crystal structure of aminomethyltransferase in complex with dihydrolipoyl-H-protein of the glycine cleavage system: implications for recognition of lipoyl protein substrate, disease-related mutations, and reaction mechanism

Aminomethyltransferase, a component of the glycine cleavage system termed T-protein, reversibly catalyzes the degradation of the aminomethyl moiety of glycine attached to the lipoate cofactor of H-protein, resulting in the production of ammonia, 5,10-methylenetetrahydrofolate, and dihydrolipoate-bearing H-protein in the presence of tetrahydrofolate. Several mutations in the human T-protein gene are known to cause nonketotic hyperglycinemia. Here, we report the crystal structure of Escherichia coli T-protein in complex with dihydrolipoate-bearing H-protein and 5-methyltetrahydrofolate, a complex mimicking the ternary complex in the reverse reaction. The structure of the complex shows a highly interacting intermolecular interface limited to a small area and the protein-bound dihydrolipoyllysine arm inserted into the active site cavity of the T-protein. Invariant Arg(292) of the T-protein is essential for complex assembly. The structure also provides novel insights in understanding the disease-causing mutations, in addition to the disease-related impairment in the cofactor-enzyme interactions reported previously. Furthermore, structural and mutational analyses suggest that the reversible transfer of the methylene group between the lipoate and tetrahydrofolate should proceed through the electron relay-assisted iminium intermediate formation.

KAAS: an automatic genome annotation and pathway reconstruction server

The number of complete and draft genomes is rapidly growing in recent years, and it has become increasingly important to automate the identification of functional properties and biological roles of genes in these genomes. In the KEGG database, genes in complete genomes are annotated with the KEGG orthology (KO) identifiers, or the K numbers, based on the best hit information using Smith-Waterman scores as well as by the manual curation. Each K number represents an ortholog group of genes, and it is directly linked to an object in the KEGG pathway map or the BRITE functional hierarchy. Here, we have developed a web-based server called KAAS (KEGG Automatic Annotation Server: http://www.genome.jp/kegg/kaas/) i.e. an implementation of a rapid method to automatically assign K numbers to genes in the genome, enabling reconstruction of KEGG pathways and BRITE hierarchies. The method is based on sequence similarities, bi-directional best hit information and some heuristics, and has achieved a high degree of accuracy when compared with the manually curated KEGG GENES database.

Extracting Sequence Motifs and the Phylogenetic Features of SNARE-Dependent Membrane Traffic

The SNARE proteins are required for membrane fusion during intracellular vesicular transport and for its specificity. Only the unique combination of SNARE proteins (cognates) can be bound and can lead to membrane fusion, although the characteristics of the possible specificity of the binding combinations encoded in the SNARE sequences have not yet been determined. We discovered by whole genome sequence analysis that sequence motifs (conserved sequences) in the SNARE motif domains for each protein group correspond to localization sites or transport pathways. We claim that these motifs reflect the specificity of the binding combinations of SNARE motif domains. Using these motifs, we could classify SNARE proteins from 48 organisms into their localization sites or transport pathways. The classification result shows that more than 10 SNARE subgroups are kingdom specific and that the SNARE paralogs involved in the plasma membrane-related transport pathways have developed greater variations in higher animals and higher plants than those involved in the endoplasmic reticulum-related transport pathways throughout eukaryotic evolution.

The repertoire of desaturases for unsaturated fatty acid synthesis in 397 genomes

Fatty acids are the major components of membrane molecules. The fact that unsaturated fatty acids play multiple important roles, physically and biologically, means that the ratio of unsaturated to saturated fatty acids in the membrane needs to be strictly regulated to maintain cellular homeostasis. The most ubiquitous and widespread modification to fatty acids, which results in a great diversity of different structures, is the insertion of double bonds. Fatty acid desaturases directly introduce regioselective double bonds into fatty acids. A phylogenetic analysis of desaturases suggests that the sequences of these proteins include a highly conserved domain, which determines the differences in specificity and regioselectivity found in these enzymes. In this study, we performed a systematic analysis of fatty acid desaturases found in the genomic data of 397 organisms. We obtained a set of desaturases clustered by regioselectivity using a hierarchal clustering analysis.

Development of a streak-camera-based time-resolved microscope fluorimeter and its application to studies of membrane fusion in single cells

A time-resolved microscope fluorimeter based on a synchroscan streak camera and a fast pulsed laser system has been developed to measure the fluorescence lifetime decay under the fluorescence microscope. This system allows one to measure the nanosecond fluorescence lifetimes of fluorophores in a small spot (0.8-6.3 microns diameter) in single cultured cells under a fluorescence microscope, while the cells are being viewed under a high-power objective lens. A signal acquisition time between a second and a minute was usually sufficient to obtain fluorescence decay curves with good quality for 10(3)-10(5) fluorophores localized in 1 microns 2 domain. A signal-to-noise ratio better than 30 was obtained for approximately 30,000 fluorescein-labeled band 3 molecules in a 2 microns 2 region in a single human erythrocyte ghost after signal accumulation for 30 s. The measured lifetimes for a variety of fluorescent probes attached to proteins in solution and lipids in liposomes showed a good agreement with those measured in a cuvette under standard conditions by time-correlated single photon counting. With the development of this instrument, microscope fluorimetry has become a practical, straightforward, quantitative technique for investigation of molecular processes in single cells in culture. Time-resolved microscope fluorimetry has been applied to observe fusion of liposomes in vitro and that of endosomes in single cells by monitoring resonance energy transfer. Inspection of individual liposomes and endosomes revealed the extent of fusion for each vesicle. Since the use of time-resolved microscope fluorimetry eliminates the need for subcellular fractionation or the complex correction procedures in steady-state microfluorimetry, it greatly simplifies the assay for endosome fusion in vivo. The results showed that extensive fusion of sequentially formed endosomes takes place all over the cell matrix in cultured cells. This suggests that extensive fusion with incoming endosomes takes place in many endosomal compartments, possibly sorting organelles, or that the early endosomes fuse with the preexisting network of tubular cisternae of the endosomal compartment at many points in the network. It is concluded that time-resolved microscope fluorimetry is a powerful noninvasive technique for studies of in situ biochemistry and biophysics using cells and tissues.