Search

The quantitation of proteins using shotgun proteomics has gained popularity in the last decades, simplifying sample handling procedures, removing extensive protein separation steps and achieving a relatively high throughput readout. The process starts with the digestion of the protein mixture into peptides, which are then separated by liquid chromatography and sequenced by tandem mass spectrometry (MS/MS). At the end of the workflow, recovering the identity of the proteins originally present in the sample is often a difficult and ambiguous process, because more than one protein identifier may match a set of peptides identified from the MS/MS spectra. To address this identification problem, many MS/MS data processing software tools combine all plausible protein identifiers matching a common set of peptides into a protein group. However, this solution introduces new challenges in studies with multiple experimental runs, which can be characterized by three main factors: i) protein groups’ identifiers are local, i.e., they vary run to run, ii) the composition of each group may change across runs, and iii) the supporting evidence of proteins within each group may also change across runs. Since in general there is no conclusive evidence about the absence of proteins in the groups, protein groups need to be linked across different runs in subsequent statistical analyses. We propose an algorithm, called Protein Group Code Algorithm (PGCA), to link groups from multiple experimental runs by forming global protein groups from connected local groups. The algorithm is computationally inexpensive and enables the connection and analysis of lists of protein groups across runs needed in biomarkers studies. We illustrate the identification problem and the stability of the PGCA mapping using 65 iTRAQ experimental runs. Further, we use two biomarker studies to show how PGCA enables the discovery of relevant candidate protein group markers with similar but non-identical compositions in different runs. [ABSTRACT FROM AUTHOR]

Localized hypoxia in solid tumors activates transcriptional programs that promote the metastatic transformation of cells. Like hypoxia-inducible hyper-vascularization, loss of the retinoblastoma protein (Rb) is a trait common to advanced stages of tumor progression in many metastatic cancers. However, no link between the role of Rb and hypoxia-driven metastatic processes has been established. We demonstrated that Rb is a key mediator of the hypoxic response mediated by HIF1α/β, the master regulator of the hypoxia response, and its essential co-activator, the thyroid hormone receptor/retinoblastoma-interacting protein (TRIP230). Furthermore, loss of Rb unmasks the full co-activation potential of TRIP230. Using small inhibitory RNA approaches in vivo, we established that Rb attenuates the normal physiological response to hypoxia by HIF1α. Notably, loss of Rb results in hypoxia-dependent biochemical changes that promote acquisition of an invasive phenotype in MCF7 breast cancer cells. In addition, Rb is present in HIF1α-ARNT/HIF1β transcriptional complexes associated with TRIP230 as determined by co-immuno-precipitation, GST-pull-down and ChIP assays. These results demonstrate that Rb is a negative modulator of hypoxia-regulated transcription by virtue of its direct effects on the HIF1 complex. This work represents the first link between the functional ablation of Rb in tumor cells and HIF1α-dependent transcriptional activation and invasion. [ABSTRACT FROM AUTHOR]

The activated AHR/ARNT complex (AHRC) regulates the expression of target genes upon exposure to environmental contaminants such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Importantly, evidence has shown that TCDD represses estrogen receptor (ER) target gene activation through the AHRC. Our data indicates that AHR and ARNT act independently from each other at non-dioxin response element sites. Therefore, we sought to determine the specific functions of AHR and ARNT in estrogen-dependent signaling in human MCF7 breast cancer and human ECC-1 endometrial carcinoma cells. Knockdown of AHR with siRNA abrogates dioxin-inducible repression of estrogen-dependent gene transcription. Intriguingly, knockdown of ARNT does not effect TCDD-mediated repression of estrogen-regulated transcription, suggesting that AHR represses ER function independently of ARNT. This theory is supported by the ability of the selective AHR modulator 39,49-dimethoxy-a-naphthoflavone (DiMNF) to repress estrogen-inducible transcription. Furthermore, basal and estrogen-activated transcription of the genes encoding cathepsin-D and pS2 are down-regulated in MCF7 cells but upregulated in ECC-1 cells in response to loss of ARNT. These responses are mirrored at the protein level with cathepsin-D. Furthermore, knock-down of ARNT led to opposite but corresponding changes in estrogen-stimulated proliferation in both MCF7 and ECC-1 cells. We have obtained experimental evidence demonstrating a dioxin-dependent repressor function for AHR and a dioxin-independent co-activator/co-repressor function for ARNT in estrogen signalling. These results provide us with further insight into the mechanisms of transcription factor crosstalk and putative therapeutic targets in estrogenpositive cancers. [ABSTRACT FROM AUTHOR]

Purpose: To validate a previously characterized 10-gene signature in prostate cancer with implication to distinguish aggressive and indolent disease within low and intermediate patients’ risk groups.Methods: A case-control study design used to select 545 patients from the Mayo clinic tumor registry who underwent radical prostatectomy. A training set from this cohort (n = 359) was used to build a 10-gene model, based on high-dimensional discriminant analysis (HDDA10) to predict several endpoints of clinical patients’ outcome. An independent set (n = 219) from the same institution was used as validation set.Results: HDDA10 showed significant performance for predicting metastasis (Mets) (AUC 0.68, p = 6.4E - 6) and biochemical recurrence (BCR) (AUC = 0.65, p = 0.003) in the validation set outperforming Gleason grade grouping (GG) for BCR (AUC 0.57, p = 0.03) and with comparable performance for Mets endpoint (GG AUC 0.66, p = 8.1E - 5). HDDA10 prognostic significance was superior to any clinical-pathological parameter within GG2 + 3 (GS7) patients achieving an AUC of 0.74 (p = 0.0037) for BCR compared to Gleason pattern 4 (AUC 0.64) (p = 0.015) and AUC for Mets of 0.68 versus AUC of 0.65 for Gleason pattern 4 (p = 0.01). HDDA10 remained significant for both BCR and Mets in multivariate analysis, suggesting that it can be used to increase accuracy in stratifying patients eligible for active surveillance.Conclusion: HDDA10 is of added value to GG and other clinical-pathological parameters in predicting BCR and Mets endpoint, especially in the low to intermediate patients’ risk groups. HDDA10 prognostic value should be further validated prospectively in stratifying patients specifically in low to intermediate GS (GG1-2), such as active surveillance programs. [ABSTRACT FROM AUTHOR]

Background: Gene network inference (GNI) algorithms can be used to identify sets of coordinately expressed genes, termed network modules from whole transcriptome gene expression data. The identification of such modules has become a popular approach to systems biology, with important applications in translational research. Although diverse computational and statistical approaches have been devised to identify such modules, their performance behavior is still not fully understood, particularly in complex human tissues. Given human heterogeneity, one important question is how the outputs of these computational methods are sensitive to the input sample set, or stability. A related question is how this sensitivity depends on the size of the sample set. We describe here the SABRE (Similarity Across Bootstrap RE-sampling) procedure for assessing the stability of gene network modules using a re-sampling strategy, introduce a novel criterion for identifying stable modules, and demonstrate the utility of this approach in a clinically-relevant cohort, using two different gene network module discovery algorithms. Results: The stability of modules increased as sample size increased and stable modules were more likely to be replicated in larger sets of samples. Random modules derived from permutated gene expression data were consistently unstable, as assessed by SABRE, and provide a useful baseline value for our proposed stability criterion. Gene module sets identified by different algorithms varied with respect to their stability, as assessed by SABRE. Finally, stable modules were more readily annotated in various curated gene set databases. Conclusions: The SABRE procedure and proposed stability criterion may provide guidance when designing systems biology studies in complex human disease and tissues. [ABSTRACT FROM AUTHOR]

PURPOSE: Multiparametric magnetic resonance imaging (mpMRI) is used widely for prostate cancer (PCa) evaluation. Approximately 35% of aggressive tumors, however, are not visible on mpMRI. We sought to identify the molecular alterations associated with mpMRI-invisible tumors and determine whether mpMRI visibility is associated with PCa prognosis. METHODS: Discovery and validation cohorts included patients who underwent mpMRI before radical prostatectomy and were found to harbor both mpMRI-visible (Prostate Imaging and Reporting Data System 3 to 5) and -invisible (Prostate Imaging and Reporting Data System 1 or 2) foci on surgical pathology. Next-generation sequencing was performed to determine differential gene expression between mpMRI-visible and -invisible foci. A genetic signature for tumor mpMRI visibility was derived in the discovery cohort and assessed in an independent validation cohort. Its association with long-term oncologic outcomes was evaluated in a separate testing cohort. RESULTS: The discovery cohort included 10 patients with 26 distinct PCa foci on surgical pathology, of which 12 (46%) were visible and 14 (54%) were invisible on preoperative mpMRI. Next-generation sequencing detected prioritized genetic mutations in 14 (54%) tumor foci (n = 8 mpMRI visible, n = 6 mpMRI invisible). A nine-gene signature (composed largely of cell organization/structure genes) associated with mpMRI visibility was derived (area under the curve = 0.89), and the signature predicted MRI visibility with 75% sensitivity and 100% specificity (area under the curve = 0.88) in the validation cohort. In the testing cohort (n = 375, median follow-up 8 years) there was no significant difference in biochemical recurrence, distant metastasis, or cancer-specific mortality in patients with predicted mpMRI-visible versus -invisible tumors (all P >.05). CONCLUSION: Compared with mpMRI-invisible disease, mpMRI-visible tumors are associated with underexpression of cellular organization genes. mpMRI visibility does not seem to be predictive of long-term cancer outcomes, highlighting the need for biopsy strategies that detect mpMRI-invisible tumors. [ABSTRACT FROM AUTHOR]

Objectives To determine the prognostic potential of a 24-gene signature, Sig24, for identifying patients with prostate cancer who are at risk of developing metastases or of prostate cancer-specific mortality (PCSM) after radical prostatectomy (RP). Patients and Methods Sig24 scores were calculated from previously collected gene expression microarray data from the Cleveland Clinic and Mayo Clinic (I and II). The performance of Sig24 was determined using time-dependent c-index analysis, Cox proportional hazards regression and Kaplan-Meier survival analysis. Results Higher Sig24 scores were significantly associated with higher pathological Gleason scores in all three cohorts. Analysis of the Mayo Clinic II cohort, which included time-to-event information, indicated that patients with high Sig24 scores also had a higher risk of developing metastasis (hazard ratio [HR] 3.78, 95% confidence interval [CI]: 1.96-7.29; P < 0.001) or of PCSM (HR 6.54, 95% CI: 2.16-19.83; P < 0.001). Conclusions The findings of the present study show the applicability of Sig24 for the prognosis of metastasis or PCSM after RP. Future studies investigating the combination of Sig24 with available prognostic tests may provide new approaches to improve risk stratification for patients with prostate cancer. [ABSTRACT FROM AUTHOR]

Recent technical advances in the field of quantitative proteomics have stimulated a large number of biomarker discovery studies of various diseases, providing avenues for new treatments and diagnostics. However, inherent challenges have limited the successful translation of candidate biomarkers into clinical use, thus highlighting the need for a robust analytical methodology to transition from biomarker discovery to clinical implementation. We have developed an end-to-end computational proteomic pipeline for biomarkers studies. At the discovery stage, the pipeline emphasizes different aspects of experimental design, appropriate statistical methodologies, and quality assessment of results. At the validation stage, the pipeline focuses on the migration of the results to a platform appropriate for external validation, and the development of a classifier score based on corroborated protein biomarkers. At the last stage towards clinical implementation, the main aims are to develop and validate an assay suitable for clinical deployment, and to calibrate the biomarker classifier using the developed assay. The proposed pipeline was applied to a biomarker study in cardiac transplantation aimed at developing a minimally invasive clinical test to monitor acute rejection. Starting with an untargeted screening of the human plasma proteome, five candidate biomarker proteins were identified. Rejection-regulated proteins reflect cellular and humoral immune responses, acute phase inflammatory pathways, and lipid metabolism biological processes. A multiplex multiple reaction monitoring mass-spectrometry (MRM-MS) assay was developed for the five candidate biomarkers and validated by enzyme-linked immune-sorbent (ELISA) and immunonephelometric assays (INA). A classifier score based on corroborated proteins demonstrated that the developed MRM-MS assay provides an appropriate methodology for an external validation, which is still in progress. Plasma proteomic biomarkers of acute cardiac rejection may offer a relevant post-transplant monitoring tool to effectively guide clinical care. The proposed computational pipeline is highly applicable to a wide range of biomarker proteomic studies. [ABSTRACT FROM AUTHOR]

Background: Biomarker panels derived separately from genomic and proteomic data and with a variety of computational methods have demonstrated promising classification performance in various diseases. An open question is how to create effective proteo-genomic panels. The framework of ensemble classifiers has been applied successfully in various analytical domains to combine classifiers so that the performance of the ensemble exceeds the performance of individual classifiers. Using blood-based diagnosis of acute renal allograft rejection as a case study, we address the following question in this paper: Can acute rejection classification performance be improved by combining individual genomic and proteomic classifiers in an ensemble? Results: The first part of the paper presents a computational biomarker development pipeline for genomic and proteomic data. The pipeline begins with data acquisition (e.g., from bio-samples to microarray data), quality control, statistical analysis and mining of the data, and finally various forms of validation. The pipeline ensures that the various classifiers to be combined later in an ensemble are diverse and adequate for clinical use. Five mRNA genomic and five proteomic classifiers were developed independently using single time-point blood samples from 11 acute-rejection and 22 non-rejection renal transplant patients. The second part of the paper examines five ensembles ranging in size from two to 10 individual classifiers. Performance of ensembles is characterized by area under the curve (AUC), sensitivity, and specificity, as derived from the probability of acute rejection for individual classifiers in the ensemble in combination with one of two aggregation methods: (1) Average Probability or (2) Vote Threshold. One ensemble demonstrated superior performance and was able to improve sensitivity and AUC beyond the best values observed for any of the individual classifiers in the ensemble, while staying within the range of observed specificity. The Vote Threshold aggregation method achieved improved sensitivity for all 5 ensembles, but typically at the cost of decreased specificity. Conclusion: Proteo-genomic biomarker ensemble classifiers show promise in the diagnosis of acute renal allograft rejection and can improve classification performance beyond that of individual genomic or proteomic classifiers alone. Validation of our results in an international multicenter study is currently underway. [ABSTRACT FROM AUTHOR]

Background Treatment-induced neuroendocrine prostate cancer (tNEPC) is an aggressive variant of late-stage metastatic castrate-resistant prostate cancer that commonly arises through neuroendocrine transdifferentiation (NEtD). Treatment options are limited, ineffective, and, for most patients, result in death in less than a year. We previously developed a first-in-field patient-derived xenograft (PDX) model of NEtD. Longitudinal deep transcriptome profiling of this model enabled monitoring of dynamic transcriptional changes during NEtD and in the context of androgen deprivation. Long non-coding RNA (lncRNA) are implicated in cancer where they can control gene regulation. Until now, the expression of lncRNAs during NEtD and their clinical associations were unexplored. Results We implemented a next-generation sequence analysis pipeline that can detect transcripts at low expression levels and built a genome-wide catalogue (n = 37,749) of lncRNAs. We applied this pipeline to 927 clinical samples and our high-fidelity NEtD model LTL331 and identified 821 lncRNAs in NEPC. Among these are 122 lncRNAs that robustly distinguish NEPC from prostate adenocarcinoma (AD) patient tumours. The highest expressed lncRNAs within this signature are H19, LINC00617, and SSTR5-AS1. Another 742 are associated with the NEtD process and fall into four distinct patterns of expression (NEtD lncRNA Class I, II, III, and IV) in our PDX model and clinical samples. Each class has significant (z-scores >2) and unique enrichment for transcription factor binding site (TFBS) motifs in their sequences. Enriched TFBS include (1) TP53 and BRN1 in Class I, (2) ELF5, SPIC, and HOXD1 in Class II, (3) SPDEF in Class III, (4) HSF1 and FOXA1 in Class IV, and (5) TWIST1 when merging Class III with IV. Common TFBS in all NEtD lncRNA were also identified and include E2F, REST, PAX5, PAX9, and STAF. Interrogation of the top deregulated candidates (n = 100) in radical prostatectomy adenocarcinoma samples with long-term follow-up (median 18 years) revealed significant clinicopathological associations. Specifically, we identified 25 that are associated with rapid metastasis following androgen deprivation therapy (ADT). Two of these lncRNAs (SSTR5-AS1 and LINC00514) stratified patients undergoing ADT based on patient outcome. Discussion To date, a comprehensive characterization of the dynamic landscape of lncRNAs during the NEtD process has not been performed. A temporal analysis of the PDX-based NEtD model has for the first time provided this dynamic landscape. TFBS analysis identified NEPC-related TF motifs present within the NEtD lncRNA sequences, suggesting functional roles for these lncRNAs in NEPC pathogenesis. Furthermore, select NEtD lncRNAs appear to be associated with metastasis and patients receiving ADT. Treatment-related metastasis is a clinical consequence of NEPC tumours. Top candidate lncRNAs FENDRR, H19, LINC00514, LINC00617, and SSTR5-AS1 identified in this study are implicated in the development of NEPC. We present here for the first time a genome-wide catalogue of NEtD lncRNAs that characterize the transdifferentiation process and a robust NEPC lncRNA patient expression signature. To accomplish this, we carried out the largest integrative study that applied a PDX NEtD model to clinical samples. These NEtD and NEPC lncRNAs are strong candidates for clinical biomarkers and therapeutic targets and warrant further investigation. [ABSTRACT FROM AUTHOR]