Associations between fungal and bacterial microbiota of airways and asthma endotypes.

BACKGROUND: The relationship between asthma, atopy, and underlying type 2 (T2) airway inflammation is complex. Although the bacterial airway microbiota is known to differ in asthmatic patients, the fungal and bacterial markers that discriminate T2-high (eosinophilic) and T2-low (neutrophilic/mixed-inflammation) asthma and atopy are still incompletely identified. OBJECTIVES: The aim of this study was to demonstrate the fungal microbiota structure of airways in asthmatic patients associated with T2 inflammation, atopy, and key clinical parameters. METHODS: We collected endobronchial brush (EB) and bronchoalveolar lavage (BAL) samples from 39 asthmatic patients and 19 healthy subjects followed by 16S gene and internal transcribed spacer-based microbiota sequencing. The microbial sequences were classified into exact sequence variants. The T2 phenotype was defined by using a blood eosinophil count with a threshold of 300 cells/µL. RESULTS: Fungal diversity was significantly lower in EB samples from patients with T2-high compared with T2-low inflammation; key fungal genera enriched in patients with T2-high inflammation included Trichoderma species, whereas Penicillium species was enriched in patients with atopy. In BAL fluid samples the dominant genera were Cladosporium, Fusarium, Aspergillus, and Alternaria. Using generalized linear models, we identified significant associations between specific fungal exact sequence variants and FEV1, fraction of exhaled nitric oxide values, BAL fluid cell counts, and corticosteroid use. Investigation of interkingdom (bacterial-fungal) co-occurrence patterns revealed different topologies between asthmatic patients and healthy control subjects. Random forest models with fungal classifiers predicted asthma status with 75% accuracy for BAL fluid samples and 80% accuracy for EB samples. CONCLUSIONS: We demonstrate clear differences in bacterial and fungal microbiota in asthma-associated phenotypes. Our study provides additional support for considering microbial signatures in delineating asthma phenotypes.

Integrative genomics analysis reveals silencing of beta-adrenergic signaling by polycomb in prostate cancer.

The Polycomb group (PcG) protein EZH2 possesses oncogenic properties for which the underlying mechanism is unclear. We integrated in vitro cell line, in vivo tumor profiling, and genome-wide location data to nominate key targets of EZH2. One of the candidates identified was ADRB2 (Adrenergic Receptor, Beta-2), a critical mediator of beta-adrenergic signaling. EZH2 is recruited to the ADRB2 promoter and represses ADRB2 expression. ADRB2 inhibition confers cell invasion and transforms benign prostate epithelial cells, whereas ADRB2 overexpression counteracts EZH2-mediated oncogenesis. ADRB2 is underexpressed in metastatic prostate cancer, and clinically localized tumors that express lower levels of ADRB2 exhibit a poor prognosis. Taken together, we demonstrate the power of integrating multiple diverse genomic data to decipher targets of disease-related genes.

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.

Maternal asthma and microRNA regulation of soluble HLA-G in the airway.

BACKGROUND: We previously reported an interaction between maternal asthma and the child's HLA-G genotype on the child's subsequent risk for asthma. The implicated single nucleotide polymorphism at +3142 disrupted a target site for the microRNA (miR)-152 family. We hypothesized that the interaction effect might be mediated by these miRs. OBJECTIVE: The objective of this study was to test this hypothesis in adults with asthma who are a subset of the same subjects who participated in our earlier family-based studies. METHODS: We measured soluble HLA-G (sHLA-G) concentrations in bronchoalveolar lavage fluid (n = 36) and plasma (n = 57) from adult asthmatic subjects with and without a mother with asthma, and HLA-G and miR-152 family (miR-148a, miR-148b, and miR-152) transcript levels in airway epithelial cells from the same subjects. RESULTS: miR-148b levels were significantly increased in airway epithelial cells from asthmatic subjects with an asthmatic mother compared with those seen in asthmatic subjects without an asthmatic mother, and +3142 genotypes were associated with sHLA-G concentrations in bronchoalveolar lavage fluid among asthmatic subjects with an asthmatic mother but not among those with a nonasthmatic mother. Neither effect was observed in the plasma (sHLA-G) or white blood cells (miRNA). CONCLUSION: These combined results are consistent with +3142 allele-specific targeting of HLA-G by the miR-152 family and support our hypothesis that miRNA regulation of sHLA-G in the airway is influenced by both the asthma status of the subject's mother and the subject's genotype. Moreover, we demonstrate that the effects of maternal asthma on the gene regulatory landscape in the airways of the mother's children persist into adulthood.

Human leukocyte antigen-G expression in differentiated human airway epithelial cells: lack of modulation by Th2-associated cytokines.

BACKGROUND: Human leukocyte antigen (HLA)-G is a nonclassical class I antigen with immunomodulatory roles including up-regulation of suppressor T regulatory lymphocytes. HLA-G was recently identified as an asthma susceptibility gene, and expression of a soluble isoform, HLA-G5, has been demonstrated in human airway epithelium. Increased presence of HLA-G5 has been demonstrated in bronchoalveolar lavage fluid recovered from patients with mild asthma; this suggests a role for this isoform in modulating airway inflammation though the mechanisms by which this occurs is unclear. Airway inflammation associated with Th2 cytokines such as IL-4 and IL-13 is a principal feature of asthma, but whether these cytokines elicit expression of HLA-G is not known. METHODS: We examined gene and protein expression of both soluble (G5) and membrane-bound (G1) HLA-G isoforms in primary differentiated human airway epithelial cells collected from normal lungs and grown in air-liquid interface culture. Cells were treated with up to 10 ng/ml of either IL-4, IL-5, or IL-13, or 100 ng/ml of the immunomodulatory cytokine IL-10, or 10,000 U/ml of the Th1-associated cytokine interferon-beta, for 24 hr, after which RNA was isolated for evaluation by quantitative PCR and protein was collected for Western blot analysis. RESULTS: HLA-G5 but not G1 was present in dAEC as demonstrated by quantitative PCR, western blot and confocal microscopy. Neither G5 nor G1 expression was increased by the Th2-associated cytokines IL-4, IL-5 or IL-13 over 24 hr, nor after treatment with IL-10, but was increased 4.5 ± 1.4 fold after treatment with 10,000 U/ml interferon-beta. CONCLUSIONS: These data demonstrate the constitutive expression of a T lymphocyte regulatory molecule in differentiated human airway epithelial cells that is not modulated by Th2-associated cytokines.

Integrative genomic and proteomic analysis of prostate cancer reveals signatures of metastatic progression.

Molecular profiling of cancer at the transcript level has become routine. Large-scale analysis of proteomic alterations during cancer progression has been a more daunting task. Here, we employed high-throughput immunoblotting in order to interrogate tissue extracts derived from prostate cancer. We identified 64 proteins that were altered in prostate cancer relative to benign prostate and 156 additional proteins that were altered in metastatic disease. An integrative analysis of this compendium of proteomic alterations and transcriptomic data was performed, revealing only 48%-64% concordance between protein and transcript levels. Importantly, differential proteomic alterations between metastatic and clinically localized prostate cancer that mapped concordantly to gene transcripts served as predictors of clinical outcome in prostate cancer as well as other solid tumors.

Distinct classes of chromosomal rearrangements create oncogenic ETS gene fusions in prostate cancer.

Recently, we identified recurrent gene fusions involving the 5' untranslated region of the androgen-regulated gene TMPRSS2 and the ETS (E26 transformation-specific) family genes ERG, ETV1 or ETV4 in most prostate cancers. Whereas TMPRSS2-ERG fusions are predominant, fewer TMPRSS2-ETV1 cases have been identified than expected on the basis of the frequency of high (outlier) expression of ETV1 (refs 3-13). Here we explore the mechanism of ETV1 outlier expression in human prostate tumours and prostate cancer cell lines. We identified previously unknown 5' fusion partners in prostate tumours with ETV1 outlier expression, including untranslated regions from a prostate-specific androgen-induced gene (SLC45A3) and an endogenous retroviral element (HERV-K_22q11.23), a prostate-specific androgen-repressed gene (C15orf21), and a strongly expressed housekeeping gene (HNRPA2B1). To study aberrant activation of ETV1, we identified two prostate cancer cell lines, LNCaP and MDA-PCa 2B, that had ETV1 outlier expression. Through distinct mechanisms, the entire ETV1 locus (7p21) is rearranged to a 1.5-megabase prostate-specific region at 14q13.3-14q21.1 in both LNCaP cells (cryptic insertion) and MDA-PCa 2B cells (balanced translocation). Because the common factor of these rearrangements is aberrant ETV1 overexpression, we recapitulated this event in vitro and in vivo, demonstrating that ETV1 overexpression in benign prostate cells and in the mouse prostate confers neoplastic phenotypes. Identification of distinct classes of ETS gene rearrangements demonstrates that dormant oncogenes can be activated in prostate cancer by juxtaposition to tissue-specific or ubiquitously active genomic loci. Subversion of active genomic regulatory elements may serve as a more generalized mechanism for carcinoma development. Furthermore, the identification of androgen-repressed and insensitive 5' fusion partners may have implications for the anti-androgen treatment of advanced prostate cancer.

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.

Expression of IL-4/IL-13 receptors in differentiating human airway epithelial cells.

IL-4 and IL-13 elicit several important responses in airway epithelium including chemokine secretion and mucous secretion that may contribute to airway inflammation, cell migration, and differentiation. These cytokines have overlapping but not identical effector profiles likely due to shared subunits in their receptor complexes. These receptors are variably described in epithelial cells, and the relative expression, localization, and function of these receptors in differentiated and repairing epithelial cells are not clear. We examined IL-4/IL-13 receptor expression and localization in primary airway epithelial cells collected from normal human lungs and grown under conditions yielding both undifferentiated and differentiated cells inclusive of basal, goblet, and ciliated cell phenotypes. Gene expression of the IL-4Rα, IL-2Rγc, IL-13Rα1, and IL-13Rα2 receptor subunits increased with differentiation, but different patterns of localization and protein abundance were seen for each subunit based on both differentiation and the cell subtypes present. Increased expression of receptor subunits observed in more differentiated cells was associated with more substantial functional responses to IL-4 stimulation including increased eotaxin-3 expression and accelerated migration after injury. We demonstrate substantial differences in IL-4/IL-13 receptor subunit expression and responsiveness to IL-4 based on the extent of airway epithelial cell differentiation and suggest that these differences may have functional consequences in airway inflammation.

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.

Ultrasensitive digital quantification of cytokines and bacteria predicts septic shock outcomes.

Quantification of pathogen and host biomarkers is essential for the diagnosis, monitoring, and treatment of infectious diseases. Here, we demonstrate sensitive and rapid quantification of bacterial load and cytokines from human biological samples to generate actionable hypotheses. Our digital assay measures IL-6 and TNF-α proteins, gram-negative (GN) and gram-positive (GP) bacterial DNA, and the antibiotic-resistance gene blaTEM with femtomolar sensitivity. We use our method to characterize bronchoalveolar lavage fluid from patients with asthma, and find elevated GN bacteria and IL-6 levels compared to healthy subjects. We then analyze plasma from patients with septic shock and find that increasing levels of IL-6 and blaTEM are associated with mortality, while decreasing IL-6 levels are associated with recovery. Surprisingly, lower GN bacteria levels are associated with higher probability of death. Applying decision-tree analysis to our measurements, we are able to predict mortality and rate of recovery from septic shock with over 90% accuracy.

Altered transcriptional and chromatin responses to rhinovirus in bronchial epithelial cells from adults with asthma.

There is a life-long relationship between rhinovirus (RV) infection and the development and clinical manifestations of asthma. In this study we demonstrate that cultured primary bronchial epithelial cells from adults with asthma (n = 9) show different transcriptional and chromatin responses to RV infection compared to those without asthma (n = 9). Both the number and magnitude of transcriptional and chromatin responses to RV were muted in cells from asthma cases compared to controls. Pathway analysis of the transcriptionally responsive genes revealed enrichments of apoptotic pathways in controls but inflammatory pathways in asthma cases. Using promoter capture Hi-C we tethered regions of RV-responsive chromatin to RV-responsive genes and showed enrichment of these regions and genes at asthma GWAS loci. Taken together, our studies indicate a delayed or prolonged inflammatory state in cells from asthma cases and highlight genes that may contribute to genetic risk for asthma.

Real-time evaluation of p53 oscillatory behavior in vivo using bioluminescent imaging.

p53 is a key mediator of cellular response to stress, and, although its function has been carefully evaluated in vitro, noninvasive evaluation of the transcriptional activity of p53 in live animals has not been reported. To this end, we developed a transgenic mouse model wherein the firefly luciferase gene expression was dependent on the p53-responsive P2 promoter from the murine double minute 2 (MDM2) gene. Bioluminescence activity following ionizing radiation was shown to be dose, time, and p53 dependent. In addition, expression of both p53 and its activated form as well as the expression of p53 target genes (MDM2 and p21) correlated with bioluminescence activity. Temporal evaluation of p53 activity following ionizing radiation showed a distinct oscillatory pattern, which confirmed the oscillations observed previously in cultured cells. In addition, the kinetics of oscillations were altered by pretreatment with radiation-modifying agents. These results show the use of this mouse model in enhancing our understanding of the transcriptional role of p53 in vivo.

Chemokine expression in the early response to injury in human airway epithelial cells.

Basal airway epithelial cells (AEC) constitute stem/progenitor cells within the central airways and respond to mucosal injury in an ordered sequence of spreading, migration, proliferation, and differentiation to needed cell types. However, dynamic gene transcription in the early events after mucosal injury has not been studied in AEC. We examined gene expression using microarrays following mechanical injury (MI) in primary human AEC grown in submersion culture to generate basal cells and in the air-liquid interface to generate differentiated AEC (dAEC) that include goblet and ciliated cells. A select group of ~150 genes was in differential expression (DE) within 2-24 hr after MI, and enrichment analysis of these genes showed over-representation of functional categories related to inflammatory cytokines and chemokines. Network-based gene prioritization and network reconstruction using the PINTA heat kernel diffusion algorithm demonstrated highly connected networks that were richer in differentiated AEC compared to basal cells. Similar experiments done in basal AEC collected from asthmatic donor lungs demonstrated substantial changes in DE genes and functional categories related to inflammation compared to basal AEC from normal donors. In dAEC, similar but more modest differences were observed. We demonstrate that the AEC transcription signature after MI identifies genes and pathways that are important to the initiation and perpetuation of airway mucosal inflammation. Gene expression occurs quickly after injury and is more profound in differentiated AEC, and is altered in AEC from asthmatic airways. Our data suggest that the early response to injury is substantially different in asthmatic airways, particularly in basal airway epithelial cells.

The metabolic footprint of the airway bacterial community in cystic fibrosis.

BACKGROUND: Progressive, chronic bacterial infection of the airways is a leading cause of death in cystic fibrosis (CF). Culture-independent methods based on sequencing of the bacterial 16S rRNA gene describe a distinct microbial community that decreases in richness and diversity with disease progression. Understanding the functional characteristics of the microbial community may aid in identifying potential therapies and may assist in management, but current methods are cumbersome. Here, we demonstrate the use of an oxidative metabolic assay as a complement to sequencing methods to describe the microbiome in the airways of patients with CF. METHODS: Expectorated sputum was collected from 16 CF subjects and 8 control subjects. The Biolog Gen III Microplate was used in a community-level physiological profiling (CLPP)-based assay to examine oxidative metabolic activity. 16S rRNA V4 amplicon sequencing was used to characterize the taxonomy and diversity of the samples. Correlations were then identified among the oxidative activity and taxonomy data. In an additional paired analysis, sputum from seven CF subjects were collected at two separate clinic visits and compared for oxidative activity, taxonomy, and diversity. RESULTS: Significant differences in oxidative metabolic activity, microbial taxonomy, and diversity were found between the CF and control sputum samples. Oxidative activity correlated positively with total genera but not with other measures of diversity or taxonomy, demonstrating that the metabolic assay complements the structural aspects of the microbiome. As expected, Pseudomonas was significantly enriched in CF samples, while Streptococcus and Prevotella were similarly abundant in both CF and control samples. Paired analysis of CF samples at separate clinic visits revealed comparable oxidative activity that correlated with similar stability in taxonomy and diversity. CONCLUSIONS: The CLPP assay used in this study complements existing sequencing methods to delineate the oxidative metabolic footprint of the CF airway bacterial community. This method may be useful to study the CF microbial community over time and with changes in disease state.

Molecular profiling of human prostate tissues: insights into gene expression patterns of prostate development during puberty.

Testosterone production surges during puberty and orchestrates massive growth and reorganization of the prostate gland, and this glandular architecture is maintained thereafter throughout adulthood. Benign prostatic hyperplasia (BPH) and prostate adenocarcinoma (PCA) are common diseases in adulthood that do not develop in the absence of androgens. Our objective was to gain insight into gene expression changes of the prostate gland at puberty, a crucial juncture in prostate development that is androgen dependent. Understanding the role played by androgens in normal prostate development may provide greater insight into androgen involvement in prostatic diseases. Benign prostate tissues obtained from pubertal and adult age group cadaveric organ donors were harvested and profiled using 20,000 element cDNA microarrays. Statistical analysis of the microarray data identified 375 genes that were differentially expressed in pubertal prostates relative to adult prostates including genes such as Nkx3.1, TMEPAI, TGFBR3, FASN, ANKH, TGFBR2, FAAH, S100P, HoxB13, fibronectin, and TSC2 among others. Comparisons of pubertal and BPH expression profiles revealed a subset of genes that shared the expression pattern between the two groups. In addition, we observed that several genes from this list were previously demonstrated to be regulated by androgen and hence could also be potential in vivo targets of androgen action in the pubertal human prostate. Promoter searches revealed the presence of androgen response elements in a cohort of genes including tumor necrosis factor-alpha induced adipose related protein, which was found to be induced by androgen. In summary, this is the first report that provides a comprehensive view of the molecular events that occur during puberty in the human prostate and provides a cohort of genes that could be potential in vivo targets of androgenic action during puberty.

Development of a sensitive microarray immunoassay and comparison with standard enzyme-linked immunoassay for cytokine analysis.

Cytokine and cytokine inhibitors represent important components of the inflammatory response in patients with trauma, shock, and sepsis. Many investigators wish to quantify cytokines and it would be advantageous to measure multiple cytokines in a multiplex manner to obtain an inflammatory profile rather than a single value. Using the well-accepted standard enzyme-linked immunoassay (ELISA) as a basis, a microarray immunoassay (MI) was designed to measure 16 different human cytokines simultaneously. The MI was performed by spotting antibodies on nitrocellulose pads affixed to glass slides. Detection of the mediators was performed with biotin-conjugated antibodies followed by fluorescently labeled streptavidin. All antibodies and other reagents were purchased commercially. The MI achieved a lower limit of detection that was generally similar to traditional ELISAs (approximately 4-12 pg/mL) and also had a similar coefficient of variation. In the multiplexed MI, there was no cross reactivity between mediators. To verify the utility of the MI, cytokines and cytokine inhibitors were measured in endotoxin stimulated human blood by both ELISA and MI. Virtually identical cytokine concentrations were measured by both methods. These results describe the development of a sensitive, specific and cost-effective multiplexed microarray immunoassay that produces values similar to traditional ELISAs.

Combining urinary detection of TMPRSS2:ERG and PCA3 with serum PSA to predict diagnosis of prostate cancer.

OBJECTIVES: We sought to develop a clinical algorithm combining serum PSA with detection of TMPRSS2:ERG fusion and PCA3 in urine collected after digital rectal exam (post-DRE urine) to predict prostate cancer on subsequent biopsy. MATERIALS AND METHODS: Post-DRE urine was collected in 48 consecutive patients before prostate biopsy at 2 centers; quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used to detect PCA3 and TMPRSS2:ERG fusion transcript expression. Serum PSA was measured by clinical assay. The performance of TMPRSS2:ERG fusion, PCA3, and serum PSA as biomarkers predicting prostate cancer at biopsy was measured; a clinically practical algorithm combining serum PSA with TMPRSS2:ERG and PCA3 in post-DRE urine to predict prostate cancer was developed. RESULTS: Post-DRE urine sediment provided informative RNA in 45 patients; prostate cancer was present on subsequent biopsy in 15. TMPRSS2:ERG in post-DRE urine was associated with prostate cancer (OR = 12.02; P < 0.001). PCA3 had the highest sensitivity in predicting prostate cancer diagnosis (93%), whereas TMPRSS2:ERG had the highest specificity (87%). TMPRSS2:ERG had the greatest discriminatory value in predicting prostate cancer (AUC = 0.77 compared with 0.65 for PCA3 and 0.72 for serum PSA alone). Combining serum PSA, PCA3, and TMPRSS2:ERG in a multivariable algorithm optimized for clinical utility improved cancer prediction (AUC = 0.88; specificity = 90% at 80% sensitivity). CONCLUSIONS: A clinical algorithm specifying biopsy for all patients with PSA ≥ 10 ng/ml, while restricting biopsy among those with PSA <10 ng/ml to only those with detectable PCA3 or TMPRSS2:ERG in post-DRE urine, performed better than the individual biomarkers alone in predicting prostate cancer.