Expressed pseudogenes in the transcriptional landscape of human cancers.

Pseudogene transcripts can provide a novel tier of gene regulation through generation of endogenous siRNAs or miRNA-binding sites. Characterization of pseudogene expression, however, has remained confined to anecdotal observations due to analytical challenges posed by the extremely close sequence similarity with their counterpart coding genes. Here, we describe a systematic analysis of pseudogene "transcription" from an RNA-Seq resource of 293 samples, representing 13 cancer and normal tissue types, and observe a surprisingly prevalent, genome-wide expression of pseudogenes that could be categorized as ubiquitously expressed or lineage and/or cancer specific. Further, we explore disease subtype specificity and functions of selected expressed pseudogenes. Taken together, we provide evidence that transcribed pseudogenes are a significant contributor to the transcriptional landscape of cells and are positioned to play significant roles in cellular differentiation and cancer progression, especially in light of the recently described ceRNA networks. Our work provides a transcriptome resource that enables high-throughput analyses of pseudogene expression.

Pan AMPK Activation Protects Tubules in Rat Ischemic Acute Kidney Injury.

Acute Kidney Injury (AKI) is characterized by an abrupt decline in kidney function and has been associated with excess risks of death, kidney disease progression, and cardiovascular events. The kidney has a high energetic demand with mitochondrial health being essential to renal function and damaged mitochondria has been reported across AKI subtypes. 5' adenosine monophosphate-activated protein kinase (AMPK) activation preserves cellular energetics through improvement of mitochondrial function and biogenesis when ATP levels are low such as under ischemia-induced AKI. We developed a selective potent small molecule pan AMPK activator, compound 1, and tested its ability to increase AMPK activity and preserve kidney function during ischemia/reperfusion injury in rats. A single administration of 1 caused sustained activation of AMPK for at least 24 hours, protected against acute tubular necrosis, and reduced clinical markers of tubular injury such as NephroCheck and Fractional Excretion of Sodium (FENa). Reduction in plasma creatinine and increased Glomerular Filtration Rate (GFR) indicated preservation of kidney function. Surprisingly, we observed a strong diuretic effect of AMPK activation associated with natriuresis both with and without AKI. Our findings demonstrate that activation of AMPK leads to protection of tubular function under hypoxic/ischemic conditions which holds promise as a potential novel therapeutic approach for AKI. Significance Statement No approved pharmacological therapies currently exist for acute kidney injury. We developed Compound 1 which dose-dependently activated AMPK in the kidney and protected kidney function and tubules after ischemic renal injury in the rat. This was accompanied by natriuresis in injured as well as uninjured rats.

Heat Shock Protein Beta-1 Modifies Anterior to Posterior Purkinje Cell Vulnerability in a Mouse Model of Niemann-Pick Type C Disease.

Selective neuronal vulnerability is characteristic of most degenerative disorders of the CNS, yet mechanisms underlying this phenomenon remain poorly characterized. Many forms of cerebellar degeneration exhibit an anterior-to-posterior gradient of Purkinje cell loss including Niemann-Pick type C1 (NPC) disease, a lysosomal storage disorder characterized by progressive neurological deficits that often begin in childhood. Here, we sought to identify candidate genes underlying vulnerability of Purkinje cells in anterior cerebellar lobules using data freely available in the Allen Brain Atlas. This approach led to the identification of 16 candidate neuroprotective or susceptibility genes. We demonstrate that one candidate gene, heat shock protein beta-1 (HSPB1), promoted neuronal survival in cellular models of NPC disease through a mechanism that involved inhibition of apoptosis. Additionally, we show that over-expression of wild type HSPB1 or a phosphomimetic mutant in NPC mice slowed the progression of motor impairment and diminished cerebellar Purkinje cell loss. We confirmed the modulatory effect of Hspb1 on Purkinje cell degeneration in vivo, as knockdown by Hspb1 shRNA significantly enhanced neuron loss. These results suggest that strategies to promote HSPB1 activity may slow the rate of cerebellar degeneration in NPC disease and highlight the use of bioinformatics tools to uncover pathways leading to neuronal protection in neurodegenerative disorders.

Renal cell tumors with clear cell histology and intact VHL and chromosome 3p: a histological review of tumors from the Cancer Genome Atlas database.

Clear cell renal cell carcinoma is by far the most common form of kidney cancer; however, a number of histologically similar tumors are now recognized and considered distinct entities. The Cancer Genome Atlas published data set was queried (http://cbioportal.org) for clear cell renal cell carcinoma tumors lacking VHL gene mutation and chromosome 3p loss, for which whole-slide images were reviewed. Of the 418 tumors in the published Cancer Genome Atlas clear cell renal cell carcinoma database, 387 had VHL mutation, copy number loss for chromosome 3p, or both (93%). Of the remaining, 27/31 had whole-slide images for review. One had 3p loss based on karyotype but not sequencing, and three demonstrated VHL promoter hypermethylation. Nine could be reclassified as distinct or emerging entities: translocation renal cell carcinoma (n=3), TCEB1 mutant renal cell carcinoma (n=3), papillary renal cell carcinoma (n=2), and clear cell papillary renal cell carcinoma (n=1). Of the remaining, 6 had other clear cell renal cell carcinoma-associated gene alterations (PBRM1, SMARCA4, BAP1, SETD2), leaving 11 specimens, including 2 high-grade or sarcomatoid renal cell carcinomas and 2 with prominent fibromuscular stroma (not TCEB1 mutant). One of the remaining tumors exhibited gain of chromosome 7 but lacked histological features of papillary renal cell carcinoma. Two tumors previously reported to harbor TFE3 gene fusions also exhibited VHL mutation, chromosome 3p loss, and morphology indistinguishable from clear cell renal cell carcinoma, the significance of which is uncertain. In summary, almost all clear cell renal cell carcinomas harbor VHL mutation, 3p copy number loss, or both. Of tumors with clear cell histology that lack these alterations, a subset can now be reclassified as other entities. Further study will determine whether additional entities exist, based on distinct genetic pathways that may have implications for treatment.

Recurrent reciprocal RNA chimera involving YPEL5 and PPP1CB in chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is the most common form of leukemia in adults in the Western hemisphere. Tumor-specific chromosomal translocations, characteristic findings in several human malignancies that directly lead to malignant transformation, have not been identified in CLL. Using paired-end transcriptome sequencing, we identified recurrent and reciprocal RNA chimeras involving yippee like 5 (YPEL5) and serine/threonine-protein phosphatase PP1-beta-catalytic subunit (PPP1CB) in CLL. Two of seven index cases (28%) harbored the reciprocal RNA chimeras in our initial screening. Using quantitative real-time PCR (q real-time PCR), YPEL5/PPP1CB and PPP1CB/YPEL5 fusion transcripts were detected in 97 of 103 CLL samples (95%) but not in paired normal samples, benign lymphocytes, or various unrelated cancers. Whole-genome sequencing and Southern blotting demonstrated no evidence for a genomic fusion between YPEL5 and PPP1CB. YPEL5/PPP1CB chimera, when introduced into mammalian cells, expressed a truncated PPP1CB protein that demonstrated diminished phosphatase activity. PPP1CB silencing resulted in enhanced proliferation and colony formation of MEC1 and JVM3 cells, implying a role in the pathogenesis of mature B-cell leukemia. These studies uncover a potential role for recurrent RNA chimeras involving phosphatases in the pathogenesis of a common form of leukemia.

Integrative molecular concept modeling of prostate cancer progression.

Despite efforts to profile prostate cancer, the genetic alterations and biological processes that correlate with the observed histological progression are unclear. Using laser-capture microdissection to isolate 101 cell populations, we have profiled prostate cancer progression from benign epithelium to metastatic disease. By analyzing expression signatures in the context of over 14,000 'molecular concepts', or sets of biologically connected genes, we generated an integrative model of progression. Molecular concepts that demarcate critical transitions in progression include protein biosynthesis, E26 transformation-specific (ETS) family transcriptional targets, androgen signaling and cell proliferation. Of note, relative to low-grade prostate cancer (Gleason pattern 3), high-grade cancer (Gleason pattern 4) shows an attenuated androgen signaling signature, similar to metastatic prostate cancer, which may reflect dedifferentiation and explain the clinical association of grade with prognosis. Taken together, these data show that analyzing gene expression signatures in the context of a compendium of molecular concepts is useful in understanding cancer biology.

Primary urethral clear-cell adenocarcinoma: comprehensive analysis by surgical pathology, cytopathology, and next-generation sequencing.

Primary clear-cell adenocarcinoma of the urethra, a rare tumor that histomorphologically resembles clear-cell carcinoma of the female genital tract, occurs predominantly in women and is associated with a relatively poor prognosis. The histogenesis of this rare urethral neoplasm has not been completely resolved, but it is thought to arise from either müllerian rests or metaplastic urothelium. Herein, we present comprehensive surgical pathological and cytopathological findings from a patient with primary urethral clear-cell adenocarcinoma and describe next-generation sequencing results for this patient's unique tumor-the first such reported characterization of molecular aberrations in urethral clear-cell adenocarcinoma at the transcriptomic and genomic levels. Transcriptome analysis revealed novel gene fusion candidates, including ANKRD28-FNDC3B. Whole-exome analysis demonstrated focal copy number loss at the SMAD4 and ARID2 loci and 38 somatic mutations, including a truncating mutation in ATM and a novel nonsynonymous mutation in ALK.

Transcriptome sequencing to detect gene fusions in cancer.

Recurrent gene fusions, typically associated with haematological malignancies and rare bone and soft-tissue tumours, have recently been described in common solid tumours. Here we use an integrative analysis of high-throughput long- and short-read transcriptome sequencing of cancer cells to discover novel gene fusions. As a proof of concept, we successfully used integrative transcriptome sequencing to 're-discover' the BCR-ABL1 (ref. 10) gene fusion in a chronic myelogenous leukaemia cell line and the TMPRSS2-ERG gene fusion in a prostate cancer cell line and tissues. Additionally, we nominated, and experimentally validated, novel gene fusions resulting in chimaeric transcripts in cancer cell lines and tumours. Taken together, this study establishes a robust pipeline for the discovery of novel gene chimaeras using high-throughput sequencing, opening up an important class of cancer-related mutations for comprehensive characterization.

Metabolomic profiles delineate potential role for sarcosine in prostate cancer progression.

Multiple, complex molecular events characterize cancer development and progression. Deciphering the molecular networks that distinguish organ-confined disease from metastatic disease may lead to the identification of critical biomarkers for cancer invasion and disease aggressiveness. Although gene and protein expression have been extensively profiled in human tumours, little is known about the global metabolomic alterations that characterize neoplastic progression. Using a combination of high-throughput liquid-and-gas-chromatography-based mass spectrometry, we profiled more than 1,126 metabolites across 262 clinical samples related to prostate cancer (42 tissues and 110 each of urine and plasma). These unbiased metabolomic profiles were able to distinguish benign prostate, clinically localized prostate cancer and metastatic disease. Sarcosine, an N-methyl derivative of the amino acid glycine, was identified as a differential metabolite that was highly increased during prostate cancer progression to metastasis and can be detected non-invasively in urine. Sarcosine levels were also increased in invasive prostate cancer cell lines relative to benign prostate epithelial cells. Knockdown of glycine-N-methyl transferase, the enzyme that generates sarcosine from glycine, attenuated prostate cancer invasion. Addition of exogenous sarcosine or knockdown of the enzyme that leads to sarcosine degradation, sarcosine dehydrogenase, induced an invasive phenotype in benign prostate epithelial cells. Androgen receptor and the ERG gene fusion product coordinately regulate components of the sarcosine pathway. Here, by profiling the metabolomic alterations of prostate cancer progression, we reveal sarcosine as a potentially important metabolic intermediary of cancer cell invasion and aggressivity.

Deep sequencing reveals distinct patterns of DNA methylation in prostate cancer.

Beginning with precursor lesions, aberrant DNA methylation marks the entire spectrum of prostate cancer progression. We mapped the global DNA methylation patterns in select prostate tissues and cell lines using MethylPlex-next-generation sequencing (M-NGS). Hidden Markov model-based next-generation sequence analysis identified ∼68,000 methylated regions per sample. While global CpG island (CGI) methylation was not differential between benign adjacent and cancer samples, overall promoter CGI methylation significantly increased from ~12.6% in benign samples to 19.3% and 21.8% in localized and metastatic cancer tissues, respectively (P-value < 2 × 10(-16)). We found distinct patterns of promoter methylation around transcription start sites, where methylation occurred not only on the CGIs, but also on flanking regions and CGI sparse promoters. Among the 6691 methylated promoters in prostate tissues, 2481 differentially methylated regions (DMRs) are cancer-specific, including numerous novel DMRs. A novel cancer-specific DMR in the WFDC2 promoter showed frequent methylation in cancer (17/22 tissues, 6/6 cell lines), but not in the benign tissues (0/10) and normal PrEC cells. Integration of LNCaP DNA methylation and H3K4me3 data suggested an epigenetic mechanism for alternate transcription start site utilization, and these modifications segregated into distinct regions when present on the same promoter. Finally, we observed differences in repeat element methylation, particularly LINE-1, between ERG gene fusion-positive and -negative cancers, and we confirmed this observation using pyrosequencing on a tissue panel. This comprehensive methylome map will further our understanding of epigenetic regulation in prostate cancer progression.

Mining for regulatory programs in the cancer transcriptome.

DNA microarrays have been widely applied to cancer transcriptome analysis. The Oncomine database contains a large collection of such data, as well as hundreds of derived gene-expression signatures. We studied the regulatory mechanisms responsible for gene deregulation in these cancer signatures by searching for the coordinate regulation of genes with common transcription factor binding sites. We found that genes with binding sites for the archetypal cancer transcription factor, E2F, were disproportionately overexpressed in a wide variety of cancers, whereas genes with binding sites for other transcription factors, such as Myc-Max, c-Rel and ATF, were disproportionately overexpressed in specific cancer types. These results suggest that alterations in pathways activating these transcription factors may be responsible for the observed gene deregulation and cancer pathogenesis.

Prognostic Significance of Proteomics-Discovered Circulating Inflammatory Biomarkers in Patients With Pulmonary Arterial Hypertension.

BACKGROUND: Pulmonary arterial hypertension (PAH) ultimately leads to right ventricular failure and premature death. The identification of circulating biomarkers with prognostic utility is considered a priority. As chronic inflammation is recognized as key pathogenic driver, we sought to identify inflammation-related circulating proteins that add incremental value to current risk stratification models for long-term survival in patients with PAH. METHODS AND RESULTS: Plasma levels of 384 inflammatory proteins were measured with the proximity extension assay technology in patients with PAH (n=60) and controls with normal hemodynamics (n=28). Among these, 51 analytes were significantly overexpressed in the plasma of patients with PAH compared with controls. Cox proportional hazard analyses and C-statistics were performed to assess the prognostic value and the incremental prognostic value of differentially expressed proteins. A panel of 6 proteins (CRIM1 [cysteine rich transmembrane bone morphogenetic protein regulator 1], HGF [hepatocyte growth factor], FSTL3 [follistatin-like 3], PLAUR [plasminogen activator, urokinase receptor], CLSTN2 [calsyntenin 2], SPON1 [spondin 1]) were independently associated with death/lung transplantation at the time of PAH diagnosis after adjustment for the 2015 European Society of Cardiology/European Respiratory Society guidelines, the REVEAL (Registry to Evaluate Early and Long-Term PAH Disease Management) 2.0 risk scores, and the refined 4-strata risk assessment. CRIM1, PLAUR, FSTL3, and SPON1 showed incremental prognostic value on top of the predictive models. As determined by Western blot, FSTL3 and SPON1 were significantly upregulated in the right ventricle of patients with PAH and animal models (monocrotaline-injected and pulmonary artery banding-subjected rats). CONCLUSIONS: In addition to revealing new actors likely involved in cardiopulmonary remodeling in PAH, our screening identified promising circulating biomarkers to improve risk prediction in PAH, which should be externally confirmed.

miRConnect 2.0: identification of oncogenic, antagonistic miRNA families in three human cancers.

BACKGROUND: Based on their function in cancer micro(mi)RNAs are often grouped as either tumor suppressors or oncogenes. However, miRNAs regulate multiple tumor relevant signaling pathways raising the question whether two oncogenic miRNAs could be functional antagonists by promoting different steps in tumor progression. We recently developed a method to connect miRNAs to biological function by comparing miRNA and gene array expression data from the NCI60 cell lines without using miRNA target predictions (miRConnect). RESULTS: We have now extended this analysis to three primary human cancers (ovarian cancer, glioblastoma multiforme, and kidney renal clear cell carcinoma) available at the Cancer Genome Atlas (TCGA), and have correlated the expression of the clustered miRNAs with 158 oncogenic signatures (miRConnect 2.0). We have identified functionally antagonistic groups of miRNAs. One group (the agonists), which contains many of the members of the miR-17 family, correlated with c-Myc induced genes and E2F gene signatures. A group that was directly antagonistic to the agonists in all three primary cancers contains miR-221 and miR-222. Since both miR-17 ~ 92 and miR-221/222 are considered to be oncogenic this points to a functional antagonism of different oncogenic miRNAs. Analysis of patient data revealed that in certain patients agonistic miRNAs predominated, whereas in other patients antagonists predominated. In glioblastoma a high ratio of miR-17 to miR-221/222 was predictive of better overall survival suggesting that high miR-221/222 expression is more adverse for patients than high miR-17 expression. CONCLUSION: miRConnect 2.0 is useful for identifying activities of miRNAs that are relevant to primary cancers. The new correlation data on miRNAs and mRNAs deregulated in three primary cancers are available at miRConnect.org.

Gene Fusion Markup Language: a prototype for exchanging gene fusion data.

BACKGROUND: An avalanche of next generation sequencing (NGS) studies has generated an unprecedented amount of genomic structural variation data. These studies have also identified many novel gene fusion candidates with more detailed resolution than previously achieved. However, in the excitement and necessity of publishing the observations from this recently developed cutting-edge technology, no community standardization approach has arisen to organize and represent the data with the essential attributes in an interchangeable manner. As transcriptome studies have been widely used for gene fusion discoveries, the current non-standard mode of data representation could potentially impede data accessibility, critical analyses, and further discoveries in the near future. RESULTS: Here we propose a prototype, Gene Fusion Markup Language (GFML) as an initiative to provide a standard format for organizing and representing the significant features of gene fusion data. GFML will offer the advantage of representing the data in a machine-readable format to enable data exchange, automated analysis interpretation, and independent verification. As this database-independent exchange initiative evolves it will further facilitate the formation of related databases, repositories, and analysis tools. The GFML prototype is made available at http://code.google.com/p/gfml-prototype/. CONCLUSION: The Gene Fusion Markup Language (GFML) presented here could facilitate the development of a standard format for organizing, integrating and representing the significant features of gene fusion data in an inter-operable and query-able fashion that will enable biologically intuitive access to gene fusion findings and expedite functional characterization. A similar model is envisaged for other NGS data analyses.

Quantitative proteomic profiling of prostate cancer reveals a role for miR-128 in prostate cancer.

Multiple, complex molecular events characterize cancer development and progression. Deciphering the molecular networks that distinguish organ-confined disease from metastatic disease may lead to the identification of biomarkers of cancer invasion and disease aggressiveness. Although alterations in gene expression have been extensively quantified during neoplastic progression, complementary analyses of proteomic changes have been limited. Here we interrogate the proteomic alterations in a cohort of 15 prostate-derived tissues that included five each from adjacent benign prostate, clinically localized prostate cancer, and metastatic disease from distant sites. The experimental strategy couples isobaric tags for relative and absolute quantitation with multidimensional liquid phase peptide fractionation followed by tandem mass spectrometry. Over 1000 proteins were quantified across the specimens and delineated into clinically localized and metastatic prostate cancer-specific signatures. Included in these class-specific profiles were both proteins that were known to be dysregulated during prostate cancer progression and new ones defined by this study. Enrichment analysis of the prostate cancer-specific proteomic signature, to gain insight into the functional consequences of these alterations, revealed involvement of miR-128-a/b regulation during prostate cancer progression. This finding was validated using real time PCR analysis for microRNA transcript levels in an independent set of 15 clinical specimens. miR-128 levels were elevated in benign prostate epithelial cell lines compared with invasive prostate cancer cells. Knockdown of miR-128 induced invasion in benign prostate epithelial cells, whereas its overexpression attenuated invasion in prostate cancer cells. Taken together, our profiles of the proteomic alterations of prostate cancer progression revealed miR-128 as a potentially important negative regulator of prostate cancer cell invasion.

Chimeric transcript discovery by paired-end transcriptome sequencing.

Recurrent gene fusions are a prevalent class of mutations arising from the juxtaposition of 2 distinct regions, which can generate novel functional transcripts that could serve as valuable therapeutic targets in cancer. Therefore, we aim to establish a sensitive, high-throughput methodology to comprehensively catalog functional gene fusions in cancer by evaluating a paired-end transcriptome sequencing strategy. Not only did a paired-end approach provide a greater dynamic range in comparison with single read based approaches, but it clearly distinguished the high-level "driving" gene fusions, such as BCR-ABL1 and TMPRSS2-ERG, from potential lower level "passenger" gene fusions. Also, the comprehensiveness of a paired-end approach enabled the discovery of 12 previously undescribed gene fusions in 4 commonly used cell lines that eluded previous approaches. Using the paired-end transcriptome sequencing approach, we observed read-through mRNA chimeras, tissue-type restricted chimeras, converging transcripts, diverging transcripts, and overlapping mRNA transcripts. Last, we successfully used paired-end transcriptome sequencing to detect previously undescribed ETS gene fusions in prostate tumors. Together, this study establishes a highly specific and sensitive approach for accurately and comprehensively cataloguing chimeras within a sample using paired-end transcriptome sequencing.

AGTR1 overexpression defines a subset of breast cancer and confers sensitivity to losartan, an AGTR1 antagonist.

Breast cancer patients have benefited from the use of targeted therapies directed at specific molecular alterations. To identify additional opportunities for targeted therapy, we searched for genes with marked overexpression in subsets of tumors across a panel of breast cancer profiling studies comprising 3,200 microarray experiments. In addition to prioritizing ERBB2, we found AGTR1, the angiotensin II receptor type I, to be markedly overexpressed in 10-20% of breast cancer cases across multiple independent patient cohorts. Validation experiments confirmed that AGTR1 is highly overexpressed, in several cases more than 100-fold. AGTR1 overexpression was restricted to estrogen receptor-positive tumors and was mutually exclusive with ERBB2 overexpression across all samples. Ectopic overexpression of AGTR1 in primary mammary epithelial cells, combined with angiotensin II stimulation, led to a highly invasive phenotype that was attenuated by the AGTR1 antagonist losartan. Similarly, losartan reduced tumor growth by 30% in AGTR1-positive breast cancer xenografts. Taken together, these observations indicate that marked AGTR1 overexpression defines a subpopulation of ER-positive, ERBB2-negative breast cancer that may benefit from targeted therapy with AGTR1 antagonists, such as losartan.

HPeak: an HMM-based algorithm for defining read-enriched regions in ChIP-Seq data.

BACKGROUND: Protein-DNA interaction constitutes a basic mechanism for the genetic regulation of target gene expression. Deciphering this mechanism has been a daunting task due to the difficulty in characterizing protein-bound DNA on a large scale. A powerful technique has recently emerged that couples chromatin immunoprecipitation (ChIP) with next-generation sequencing, (ChIP-Seq). This technique provides a direct survey of the cistrom of transcription factors and other chromatin-associated proteins. In order to realize the full potential of this technique, increasingly sophisticated statistical algorithms have been developed to analyze the massive amount of data generated by this method. RESULTS: Here we introduce HPeak, a Hidden Markov model (HMM)-based Peak-finding algorithm for analyzing ChIP-Seq data to identify protein-interacting genomic regions. In contrast to the majority of available ChIP-Seq analysis software packages, HPeak is a model-based approach allowing for rigorous statistical inference. This approach enables HPeak to accurately infer genomic regions enriched with sequence reads by assuming realistic probability distributions, in conjunction with a novel weighting scheme on the sequencing read coverage. CONCLUSIONS: Using biologically relevant data collections, we found that HPeak showed a higher prevalence of the expected transcription factor binding motifs in ChIP-enriched sequences relative to the control sequences when compared to other currently available ChIP-Seq analysis approaches. Additionally, in comparison to the ChIP-chip assay, ChIP-Seq provides higher resolution along with improved sensitivity and specificity of binding site detection. Additional file and the HPeak program are freely available at http://www.sph.umich.edu/csg/qin/HPeak.

Humoral response profiling reveals pathways to prostate cancer progression.

There is considerable evidence for an association between prostate cancer development and inflammation, which results in autoantibody generation against tumor proteins. This immune system-driven amplification of the autoantibody response to intracellular antigens can serve as a sensitive tool to detect low abundance serum proteomic tumor markers for prostate cancer as well as provide insight into biological processes perturbed during cancer development. Here we examine serum humoral responses in a cohort of 34 patients with either benign prostatic hyperplasia or clinically localized prostate cancer (PCa). The experimental strategy couples multidimensional liquid-phase protein fractionation of localized and metastatic prostate cancer tissue lysates to protein microarrays and subsequent mass spectrometry. A supervised learning analysis of the humoral response arrays generated a parsimonious predictor having 78% sensitivity and 75% specificity in distinguishing PCa from benign prostatic hyperplasia in a cohort of American males with elevated prostate-specific antigen. Enrichment analysis of the PCa-specific humoral signature revealed large scale immune reprogramming mediated by STAT transcription factors and the generation of autoantibodies to enzymes involved in nitrogen metabolism. Meta-analysis of independent prostate cancer gene expression data validated the presence of STAT-induced immunomodulation. Concomitant validation of elevated levels of the nitrogen metabolism pathway was obtained by direct measurement of metabolic levels of glutamate and aspartate in prostate cancer tissues. Thus, in addition to functioning as markers in prostate cancer detection, humoral response profiles can serve as powerful tools revealing pathway dysregulation that might otherwise be suppressed by the complexity of the cancer proteome.