An integrated model for predicting KRAS dependency.

The clinical approvals of KRAS G12C inhibitors have been a revolutionary advance in precision oncology, but response rates are often modest. To improve patient selection, we developed an integrated model to predict KRAS dependency. By integrating molecular profiles of a large panel of cell lines from the DEMETER2 dataset, we built a binary classifier to predict a tumor's KRAS dependency. Monte Carlo cross validation via ElasticNet within the training set was used to compare model performance and to tune parameters α and λ. The final model was then applied to the validation set. We validated the model with genetic depletion assays and an external dataset of lung cancer cells treated with a G12C inhibitor. We then applied the model to several Cancer Genome Atlas (TCGA) datasets. The final "K20" model contains 20 features, including expression of 19 genes and KRAS mutation status. In the validation cohort, K20 had an AUC of 0.94 and accurately predicted KRAS dependency in both mutant and KRAS wild-type cell lines following genetic depletion. It was also highly predictive across an external dataset of lung cancer lines treated with KRAS G12C inhibition. When applied to TCGA datasets, specific subpopulations such as the invasive subtype in colorectal cancer and copy number high pancreatic adenocarcinoma were predicted to have higher KRAS dependency. The K20 model has simple yet robust predictive capabilities that may provide a useful tool to select patients with KRAS mutant tumors that are most likely to respond to direct KRAS inhibitors.

Peripheral Blood IL5RA Gene Expression as a Diagnostic Biomarker for Eosinophilic Esophagitis.

Eosinophilic esophagitis (EoE) is an allergic inflammatory condition of the esophagus, often diagnosed late because of its challenging symptoms and costly and invasive diagnostic methods.1,2 To address the need for more accessible biomarkers in EoE,3 we aimed to investigate the potential of whole-blood RNA expression as a noninvasive biomarker for diagnosing and monitoring EoE, hypothesizing that genetic signatures in blood could distinguish EoE cases, correlate with disease activity, and predict treatment responses.

Pre-treatment differential correlation of gene expression and response to topical steroids in eosinophilic esophagitis.

Few predictors of response to topical corticosteroid (tCS) treatment have been identified in eosinophilic esophagitis (EoE). We aimed to determine whether baseline gene expression predicts histologic response to tCS treatment for EoE. We analyzed prospectively collected samples from incident EoE cases who were treated with tCS for 8 weeks in a development cohort (prospective study) or in an independent validation cohort (clinical trial). Whole transcriptome RNA expression was determined from a baseline (pre-treatment) RNA-later preserved esophageal biopsy. Baseline expression was compared between histologic responders (<15 eos/hpf) and non-responders (≥15 eos/hpf), and differential correlation was used to assess baseline gene expression by response status. In 87 EoE cases analyzed in the development set, there were no differentially expressed genes associated with treatment response (at false discovery rate = 0.1). However, differential correlation identified a module of 22 genes with statistically significantly high pairwise correlation in non-responders (mean correlation coefficient = 0.7) compared to low correlation in responders (coefficient = 0.3). When this 22-gene module was applied to the 89 EoE cases in the independent cohort, it was not validated to predict tCS response at the 15 eos/hpf threshold (mean correlation coefficient = 0.32 in responders and 0.25 in nonresponders). Exploration of other thresholds also did not validate any modules. Though we identified a 22 gene differential correlation module measured pre-treatment that was strongly associated with subsequent histologic response to tCS in EoE, this was not validated in an independent population. Alternative methods to predict steroid response should be explored.

Utility of TERT Promoter Mutations for Cutaneous Primary Melanoma Diagnosis.

Telomerase reverse transcriptase (TERT) promoter mutations are commonly found in malignant melanomas but rare in melanocytic nevi. To assess its potential diagnostic utility for the distinction of melanoma from nevus, we determined the TERT promoter mutation status of 86 primary melanomas, 72 melanocytic nevi, and 40 diagnostically problematic melanocytic proliferations. Of the 86 melanomas, 67 (77.9%) were TERT-positive, defined as harboring a hotspot TERT promoter mutation at positions -124C>T, -124_125CC>TT, -138_139CC>TT, or -146C>T. Of the 72 nevi, only 1 (1.4%) was TERT-positive. Of the 40 diagnostically uncertain melanocytic proliferations, 2 (5.0%) were TERT-positive. TERT positivity as a test for melanoma versus nevus had an accuracy of 87.3% [95% confidence interval (CI), 81.1-92.1], a sensitivity of 77.9% (95% CI, 68.9-85.4), a specificity of 98.6% (95% CI, 95.8-100), a positive predictive value of 98.5% (95% CI, 95.6-100), and a negative predictive value of 78.9% (95% CI, 72.6-85.4). Our results indicate that hotspot TERT promoter mutation status may be a useful ancillary parameter for the diagnosis of melanoma. In particular, the high specificity of these mutations for melanoma indicates the presence of a TERT promoter mutation in a melanocytic neoplasm associated with diagnostic controversy, or uncertainty should increase concern for a melanoma.

Prevalent RNA recognition motif duplication in the human genome.

The sequence-specific recognition of RNA by proteins is mediated through various RNA binding domains, with the RNA recognition motif (RRM) being the most frequent and present in >50% of RNA-binding proteins (RBPs). Many RBPs contain multiple RRMs, and it is unclear how each RRM contributes to the binding specificity of the entire protein. We found that RRMs within the same RBP (i.e., sibling RRMs) tend to have significantly higher similarity than expected by chance. Sibling RRM pairs from RBPs shared by multiple species tend to have lower similarity than those found only in a single species, suggesting that multiple RRMs within the same protein might arise from domain duplication followed by divergence through random mutations. This finding is exemplified by a recent RRM domain duplication in DAZ proteins and an ancient duplication in PABP proteins. Additionally, we found that different similarities between sibling RRMs are associated with distinct functions of an RBP and that the RBPs tend to contain repetitive sequences with low complexity. Taken together, this study suggests that the number of RBPs with multiple RRMs has expanded in mammals and that the multiple sibling RRMs may recognize similar target motifs in a cooperative manner.

Protein identification using top-down.

In the last two years, because of advances in protein separation and mass spectrometry, top-down mass spectrometry moved from analyzing single proteins to analyzing complex samples and identifying hundreds and even thousands of proteins. However, computational tools for database search of top-down spectra against protein databases are still in their infancy. We describe MS-Align+, a fast algorithm for top-down protein identification based on spectral alignment that enables searches for unexpected post-translational modifications. We also propose a method for evaluating statistical significance of top-down protein identifications and further benchmark various software tools on two top-down data sets from Saccharomyces cerevisiae and Salmonella typhimurium. We demonstrate that MS-Align+ significantly increases the number of identified spectra as compared with MASCOT and OMSSA on both data sets. Although MS-Align+ and ProSightPC have similar performance on the Salmonella typhimurium data set, MS-Align+ outperforms ProSightPC on the (more complex) Saccharomyces cerevisiae data set.

Spon1+ inflammatory monocytes promote collagen remodeling and lung cancer metastasis through lipoprotein receptor 8 signaling.

Lung cancer is the leading cause of cancer-related deaths in the world, and non-small cell lung cancer (NSCLC) is the most common subset. We previously found that infiltration of tumor inflammatory monocytes (TIMs) into lung squamous carcinoma (LUSC) tumors is associated with increased metastases and poor survival. To further understand how TIMs promote metastases, we compared RNA-Seq profiles of TIMs from several LUSC metastatic models with inflammatory monocytes (IMs) of non-tumor-bearing controls. We identified Spon1 as upregulated in TIMs and found that Spon1 expression in LUSC tumors corresponded with poor survival and enrichment of collagen extracellular matrix signatures. We observed SPON1+ TIMs mediate their effects directly through LRP8 on NSCLC cells, which resulted in TGF-ß1 activation and robust production of fibrillar collagens. Using several orthogonal approaches, we demonstrated that SPON1+ TIMs were sufficient to promote NSCLC metastases. Additionally, we found that Spon1 loss in the host, or Lrp8 loss in cancer cells, resulted in a significant decrease of both high-density collagen matrices and metastases. Finally, we confirmed the relevance of the SPON1/LRP8/TGF-ß1 axis with collagen production and survival in patients with NSCLC. Taken together, our study describes how SPON1+ TIMs promote collagen remodeling and NSCLC metastases through an LRP8/TGF-ß1 signaling axis.

Methylation of nonessential genes in cutaneous melanoma - Rule Out hypothesis.

Differential methylation plays an important role in melanoma development and is associated with survival, progression and response to treatment. However, the mechanisms by which methylation promotes melanoma development are poorly understood. The traditional explanation of selective advantage provided by differential methylation postulates that hypermethylation of regulatory 5'-cytosine-phosphate-guanine-3' dinucleotides (CpGs) downregulates the expression of tumor suppressor genes and therefore promotes tumorigenesis. We believe that other (not necessarily alternative) explanations of the selective advantages of methylation are also possible. Here, we hypothesize that melanoma cells use methylation to shut down transcription of nonessential genes - those not required for cell survival and proliferation. Suppression of nonessential genes allows tumor cells to be more efficient in terms of energy and resource usage, providing them with a selective advantage over the tumor cells that transcribe and subsequently translate genes they do not need. We named the hypothesis the Rule Out (RO) hypothesis. The RO hypothesis predicts higher methylation of CpGs located in regulatory regions (CpG islands) of nonessential genes. It also predicts the higher methylation of regulatory CpGs linked to nonessential genes in melanomas compared to nevi and lower expression of nonessential genes in malignant (derived from melanoma) versus normal (derived from nonaffected skin) melanocytes. The analyses conducted using in-house and publicly available data found that all predictions derived from the RO hypothesis hold, providing observational support for the hypothesis.

Urinary proteome analysis of irritable bowel syndrome (IBS) symptom subgroups.

Irritable bowel syndrome (IBS) is a functional gastrointestinal (GI) disorder characterized by chronic abdominal pain associated with alterations in bowel function. Given the heterogeneity of the symptoms, multiple pathophysiologic factors are suspected to play a role. We classified women with IBS into four subgroups based on distinct symptom profiles. In-depth shotgun proteomic analysis was carried out to profile the urinary proteomes to identify possible proteins associated with these subgroups. First void urine samples with urine creatinine level≥100 mg/dL were used after excluding samples that tested positive for blood. Urine from 10 subjects representing each symptom subgroup was pooled for proteomic analysis. The urine proteome was analyzed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) using a data-independent method known as Precursor Acquisition Independent From Ion Count (PAcIFIC) that allowed extended detectable dynamic range. Differences in protein quantities were determined by peptide spectral counting followed by validation of select proteins with ELISA or a targeted single reaction monitoring (LC-SRM/MS) approach. Four IBS symptom subgroups were selected: (1) Constipation, (2) Diarrhea+Low Pain, (3) Diarrhea+High Pain, and (4) High Pain+High Psychological Distress. A fifth group consisted of Healthy Control subjects. From comparisons of quantitative spectral counting data among the symptom subgroups and controls, a total of 18 proteins that showed quantitative differences in relative abundance and possible physiological relevance to IBS were selected for further investigation. Three of the 18 proteins were chosen for validation by either ELISA or SRM. An elevated expression of gelsolin (GSN) was associated with the high pain groups. Trefoil Factor 3 (TFF3) levels were higher in IBS groups compared to controls. In this study, the IBS patients subclassified by predominant symptoms showed differences in urine proteome levels. Proteins showing distinctive changes are involved in homeostasis of intestinal function and inflammatory response. These findings warrant future studies with larger, independent cohorts to enable more extensive assessment and validation of urinary protein markers as a diagnostic tool in adults with IBS.

Characterization of the CpG Island Hypermethylated Phenotype Subclass in Primary Melanomas.

Cutaneous melanoma can be lethal even if detected at an early stage. Epigenetic profiling may facilitate the identification of aggressive primary melanomas with unfavorable outcomes. We performed clustering of whole-genome methylation data to identify subclasses that were then assessed for survival, clinical features, methylation patterns, and biological pathways. Among 89 cutaneous primary invasive melanomas, we identified three methylation subclasses exhibiting low methylation, intermediate methylation, or hypermethylation of CpG islands, known as the CpG island methylator phenotype (CIMP). CIMP melanomas occurred as early as tumor stage 1b and, compared with low-methylation melanomas, were associated with age at diagnosis ≥65 years, lentigo maligna melanoma histologic subtype, presence of ulceration, higher American Joint Committee on Cancer stage and tumor stage, and lower tumor-infiltrating lymphocyte grade (all P < 0.05). Patients with CIMP melanomas had worse melanoma-specific survival (hazard ratio = 11.84; confidence interval = 4.65‒30.20) than those with low-methylation melanomas, adjusted for age, sex, American Joint Committee on Cancer stage, and tumor-infiltrating lymphocyte grade. Genes hypermethylated in CIMP compared with those in low-methylation melanomas included PTEN, VDR, PD-L1, TET2, and gene sets related to development/differentiation, the extracellular matrix, and immunity. CIMP melanomas exhibited hypermethylation of genes important in melanoma progression and tumor immunity, and although present in some early melanomas, CIMP was associated with worse survival independent of known prognostic factors.

Rapid idiosyncratic mechanisms of clinical resistance to KRAS G12C inhibition.

BACKGROUNDThe KRAS proto-oncogene is among the most frequently mutated genes in cancer, yet for 40 years it remained an elusive therapeutic target. Recently, allosteric inhibitors that covalently bind to KRAS G12C mutations have been approved for use in lung adenocarcinomas. Although responses are observed, they are often short-lived, thus making in-depth characterization of the mechanisms of resistance of paramount importance.METHODSHere, we present a rapid-autopsy case of a patient who had a KRASG12C-mutant lung adenocarcinoma who initially responded to a KRAS G12C inhibitor but then rapidly developed resistance. Using deep-RNA and whole-exome sequencing comparing pretreatment, posttreatment, and matched normal tissues, we uncover numerous mechanisms of resistance to direct KRAS inhibition.RESULTSIn addition to decreased KRAS G12C-mutant allele frequency in refractory tumors, we also found reactivation of the MAPK pathway despite no new mutations in KRAS or its downstream mediators. Tumor cell-intrinsic and non-cell autonomous mechanisms included increased complement activation, coagulation, and tumor angiogenesis, and several lines of evidence of immunologic evasion.CONCLUSIONTogether, our findings reveal numerous mechanisms of resistance to current KRAS G12C inhibitors through enrichment of clonal populations, KRAS-independent downstream signaling, and diverse remodeling of the tumor microenvironment.FUNDINGRichard and Fran Duley, Jimmy and Kay Mann, the NIH, and the North Carolina Biotechnology Center.

Landscape of the SOX2 protein-protein interactome.

SOX2 is a key gene implicated in maintaining the stemness of embryonic and adult stem cells that appears to re-activate in several human cancers including glioblastoma multiforme. Using immunoprecipitation (IP)/MS/MS, we identified 144 proteins that are putative SOX2 interacting proteins. Of note, SOX2 was found to interact with several heterogeneous nuclear ribonucleoprotein family proteins, including HNRNPA2B1, HNRNPA3, HNRNPC, HNRNPK, HNRNPL, HNRNPM, HNRNPR, HNRNPU, as well as other ribonucleoproteins, DNA repair proteins and helicases. Gene ontology (GO) analysis revealed that the SOX2 interactome was enriched for GO terms GO:0030529 ribonucleoprotein complex and GO:0004386 helicase activity. These findings indicate that SOX2 associates with the heterogeneous nuclear ribonucleoprotein complex, suggesting a possible role for SOX2 in post-transcriptional regulation in addition to its function as a transcription factor.

Urinary proteomics evaluation in interstitial cystitis/painful bladder syndrome: a pilot study.

PURPOSE: Interstitial cystitis/painful bladder syndrome (IC/PBS) is characterized by chronic pain, pressure and discomfort felt in the pelvis or bladder. An in-depth shotgun proteomics study was carried out to profile the urinary proteome of women with IC/PBS to identify possible specific proteins and networks associated with IC%PBS. MATERIALS AND METHODS: Urine samples from ten female IC/PBS patients and ten female asymptomatic, healthy control subjects were analyzed in quadruplicate by liquid chromatography-tandem mass spectrometry (LC-MS/MS) on a hybrid linear ion trap-orbitrap mass spectrometer. Gas-phase fractionation (GPF) was used to enhance protein identification. Differences in protein quantity were determined by peptide spectral counting. RESULTS: alpha-1B-glycoprotein (A1BG) and orosomucoid-1 (ORM1) were detected in all IC%PBS patients, and > or = 60% of these patients had elevated expression of these two proteins compared to control subjects. Transthyretin (TTR) and hemopexin (HPX) were detected in all control individuals, but > or = 60% of the IC/PBS patients had decreased expression levels of these two proteins. Enrichment functional analysis showed cell adhesion and response to stimuli were down-regulated whereas response to inflammation, wounding, and tissue degradation were up-regulated in IC/PBS. Activation of neurophysiological processes in synaptic inhibition, and lack of DNA damage repair may also be key components of IC%PBS. CONCLUSION: There are qualitative and quantitative differences between the urinary proteomes of women with and without IC%PBS. We identified a number of proteins as well as pathways%networks that might contribute to the pathology of IC%PBS or result from perturbations induced by this condition.

Endothelial miR-30c suppresses tumor growth via inhibition of TGF-ß-induced Serpine1.

In tumors, extravascular fibrin forms provisional scaffolds for endothelial cell (EC) growth and motility during angiogenesis. We report that fibrin-mediated angiogenesis was inhibited and tumor growth delayed following postnatal deletion of Tgfbr2 in the endothelium of Cdh5-CreERT2 Tgfbr2fl/fl mice (Tgfbr2iECKO mice). ECs from Tgfbr2iECKO mice failed to upregulate the fibrinolysis inhibitor plasminogen activator inhibitor 1 (Serpine1, also known as PAI-1), due in part to uncoupled TGF-ß-mediated suppression of miR-30c. Bypassing TGF-ß signaling with vascular tropic nanoparticles that deliver miR-30c antagomiRs promoted PAI-1-dependent tumor growth and increased fibrin abundance, whereas miR-30c mimics inhibited tumor growth and promoted vascular-directed fibrinolysis in vivo. Using single-cell RNA-Seq and a NanoString miRNA array, we also found that subtypes of ECs in tumors showed spectrums of Serpine1 and miR-30c expression levels, suggesting functional diversity in ECs at the level of individual cells; indeed, fresh EC isolates from lung and mammary tumor models had differential abilities to degrade fibrin and launch new vessel sprouts, a finding that was linked to their inverse expression patterns of miR-30c and Serpine1 (i.e., miR-30chi Serpine1lo ECs were poorly angiogenic and miR-30clo Serpine1hi ECs were highly angiogenic). Thus, by balancing Serpine1 expression in ECs downstream of TGF-ß, miR-30c functions as a tumor suppressor in the tumor microenvironment through its ability to promote fibrin degradation and inhibit blood vessel formation.

SETD2 Haploinsufficiency for Microtubule Methylation Is an Early Driver of Genomic Instability in Renal Cell Carcinoma.

Loss of the short arm of chromosome 3 (3p) occurs early in >95% of clear cell renal cell carcinoma (ccRCC). Nearly ubiquitous 3p loss in ccRCC suggests haploinsufficiency for 3p tumor suppressors as early drivers of tumorigenesis. We previously reported methyltransferase SETD2, which trimethylates H3 histones on lysine 36 (H3K36me3) and is located in the 3p deletion, to also trimethylate microtubules on lysine 40 (αTubK40me3) during mitosis, with αTubK40me3 required for genomic stability. We now show that monoallelic, Setd2-deficient cells retaining H3K36me3, but not αTubK40me3, exhibit a dramatic increase in mitotic defects and micronuclei count, with increased viability compared with biallelic loss. In SETD2-inactivated human kidney cells, rescue with a pathogenic SETD2 mutant deficient for microtubule (αTubK40me3), but not histone (H3K36me3) methylation, replicated this phenotype. Genomic instability (micronuclei) was also a hallmark of patient-derived cells from ccRCC. These data show that the SETD2 tumor suppressor displays a haploinsufficiency phenotype disproportionately impacting microtubule methylation and serves as an early driver of genomic instability.Significance: Loss of a single allele of a chromatin modifier plays a role in promoting oncogenesis, underscoring the growing relevance of tumor suppressor haploinsufficiency in tumorigenesis. Cancer Res; 78(12); 3135-46. ©2018 AACR.

Transcriptome-wide identification and study of cancer-specific splicing events across multiple tumors.

Dysregulation of alternative splicing (AS) is one of the molecular hallmarks of cancer, with splicing alteration of numerous genes in cancer patients. However, studying splicing mis-regulation in cancer is complicated by the large noise generated from tissue-specific splicing. To obtain a global picture of cancer-specific splicing, we analyzed transcriptome sequencing data from 1149 patients in The Cancer Genome Atlas project, producing a core set of AS events significantly altered across multiple cancer types. These cancer-specific AS events are highly conserved, are more likely to maintain protein reading frame, and mainly function in cell cycle, cell adhesion/migration, and insulin signaling pathways. Furthermore, these events can serve as new molecular biomarkers to distinguish cancer from normal tissues, to separate cancer subtypes, and to predict patient survival. We also found that most genes whose expression is closely associated with cancer-specific splicing are key regulators of the cell cycle. This study uncovers a common set of cancer-specific AS events altered across multiple cancers, providing mechanistic insight into how splicing is mis-regulated in cancers.

The splicing activator DAZAP1 integrates splicing control into MEK/Erk-regulated cell proliferation and migration.

Alternative splicing of pre-messenger RNA (mRNA) is a critical stage of gene regulation in response to environmental stimuli. Here we show that DAZAP1, an RNA-binding protein involved in mammalian development and spermatogenesis, promotes inclusion of weak exons through specific recognition of diverse cis-elements. The carboxy-terminal proline-rich domain of DAZAP1 interacts with and neutralizes general splicing inhibitors, and is sufficient to activate splicing when recruited to pre-mRNA. This domain is phosphorylated by the MEK/Erk (extracellular signal-regulated protein kinase) pathway and this modification is essential for the splicing regulatory activity and the nuclear/cytoplasmic translocation of DAZAP1. Using mRNA-seq, we identify endogenous splicing events regulated by DAZAP1, many of which are involved in maintaining cell growth. Knockdown or over-expression of DAZAP1 causes a cell proliferation defect. Taken together, these studies reveal a molecular mechanism that integrates splicing control into MEK/Erk-regulated cell proliferation.

The splicing factor RBM4 controls apoptosis, proliferation, and migration to suppress tumor progression.

Splicing dysregulation is one of the molecular hallmarks of cancer. However, the underlying molecular mechanisms remain poorly defined. Here we report that the splicing factor RBM4 suppresses proliferation and migration of various cancer cells by specifically controlling cancer-related splicing. Particularly, RBM4 regulates Bcl-x splicing to induce apoptosis, and coexpression of Bcl-xL partially reverses the RBM4-mediated tumor suppression. Moreover, RBM4 antagonizes an oncogenic splicing factor, SRSF1, to inhibit mTOR activation. Strikingly, RBM4 expression is decreased dramatically in cancer patients, and the RBM4 level correlates positively with improved survival. In addition to providing mechanistic insights of cancer-related splicing dysregulation, this study establishes RBM4 as a tumor suppressor with therapeutic potential and clinical values as a prognostic factor.

Engineering RNA endonucleases with customized sequence specificities.

Specific cleavage of RNAs is critical for in vitro manipulation of RNA and for in vivo gene silencing. Here we engineer artificial site-specific RNA endonucleases to function analogously to DNA restriction enzymes. We combine a general RNA cleavage domain with a series of Pumilio/fem-3-binding factor domains that specifically recognize different 8-nucleotide RNA sequences. The resulting artificial site-specific RNA endonucleases specifically recognize RNA substrates and efficiently cleave near their binding sites. The artificial site-specific RNA endonucleases can be devised to recognize and cleave various RNA target sequences, providing a useful tool to manipulate RNAs in vitro. In addition, we generate designer artificial site-specific RNA endonucleases to specifically silence an endogenous gene in Escherichia coli, as well as a mitochondrial-encoded gene in human cells, suggesting that artificial site-specific RNA endonucleases can serve as a gene-silencing tool with designed specificity.

Precursor ion independent algorithm for top-down shotgun proteomics.

We present a precursor ion independent top-down algorithm (PIITA) for use in automated assignment of protein identifications from tandem mass spectra of whole proteins. To acquire the data, we utilize data-dependent acquisition to select protein precursor ions eluting from a C4-based HPLC column for collision induced dissociation in the linear ion trap of an LTQ-Orbitrap mass spectrometer. Gas-phase fractionation is used to increase the number of acquired tandem mass spectra, all of which are recorded in the Orbitrap mass analyzer. To identify proteins, the PIITA algorithm compares deconvoluted, deisotoped, observed tandem mass spectra to all possible theoretical tandem mass spectra for each protein in a genomic sequence database without regard for measured parent ion mass. Only after a protein is identified, is any difference in measured and theoretical precursor mass used to identify and locate post-translation modifications. We demonstrate the application of PIITA to data generated via our wet-lab approach on a Salmonella typhimurium outer membrane extract and compare these results to bottom-up analysis. From these data, we identify 154 proteins by top-down analysis, 73 of which were not identified in a parallel bottom-up analysis. We also identify 201 unique isoforms of these 154 proteins at a false discovery rate (FDR) of <1%.