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.

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.

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.

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.

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.

Transcriptome sequencing across a prostate cancer cohort identifies PCAT-1, an unannotated lincRNA implicated in disease progression.

Noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease and to reveal uncharacterized aspects of tumor biology. Here we discover 121 unannotated prostate cancer-associated ncRNA transcripts (PCATs) by ab initio assembly of high-throughput sequencing of polyA(+) RNA (RNA-Seq) from a cohort of 102 prostate tissues and cells lines. We characterized one ncRNA, PCAT-1, as a prostate-specific regulator of cell proliferation and show that it is a target of the Polycomb Repressive Complex 2 (PRC2). We further found that patterns of PCAT-1 and PRC2 expression stratified patient tissues into molecular subtypes distinguished by expression signatures of PCAT-1-repressed target genes. Taken together, our findings suggest that PCAT-1 is a transcriptional repressor implicated in a subset of prostate cancer patients. These findings establish the utility of RNA-Seq to identify disease-associated ncRNAs that may improve the stratification of cancer subtypes.

Therapeutic targeting of SPINK1-positive prostate cancer.

Gene fusions involving ETS (erythroblastosis virus E26 transformation-specific) family transcription factors are found in ~50% of prostate cancers and as such can be used as a basis for the molecular subclassification of prostate cancer. Previously, we showed that marked overexpression of SPINK1 (serine peptidase inhibitor, Kazal type 1), which encodes a secreted serine protease inhibitor, defines an aggressive molecular subtype of ETS fusion-negative prostate cancers (SPINK1+/ETS⁻, ~10% of all prostate cancers). Here, we examined the potential of SPINK1 as an extracellular therapeutic target in prostate cancer. Recombinant SPINK1 protein (rSPINK1) stimulated cell proliferation in benign RWPE as well as cancerous prostate cells. Indeed, RWPE cells treated with either rSPINK1 or conditioned medium from 22RV1 prostate cancer cells (SPINK1+/ETS⁻) significantly increased cell invasion and intravasation when compared with untreated cells. In contrast, knockdown of SPINK1 in 22RV1 cells inhibited cell proliferation, cell invasion, and tumor growth in xenograft assays. 22RV1 cell proliferation, invasion, and intravasation were attenuated by a monoclonal antibody (mAb) to SPINK1 as well. We also demonstrated that SPINK1 partially mediated its neoplastic effects through interaction with the epidermal growth factor receptor (EGFR). Administration of antibodies to SPINK1 or EGFR (cetuximab) in mice bearing 22RV1 xenografts attenuated tumor growth by more than 60 and 40%, respectively, or ~75% when combined, without affecting PC3 xenograft (SPINK1⁻/ETS⁻) growth. Thus, this study suggests that SPINK1 may be a therapeutic target in a subset of patients with SPINK1+/ETS⁻ prostate cancer. Our results provide a rationale for both the development of humanized mAbs to SPINK1 and evaluation of EGFR inhibition in SPINK1+/ETS⁻ prostate cancers.

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.

An integrated network of androgen receptor, polycomb, and TMPRSS2-ERG gene fusions in prostate cancer progression.

Chromosomal rearrangements fusing the androgen-regulated gene TMPRSS2 to the oncogenic ETS transcription factor ERG occur in approximately 50% of prostate cancers, but how the fusion products regulate prostate cancer remains unclear. Using chromatin immunoprecipitation coupled with massively parallel sequencing, we found that ERG disrupts androgen receptor (AR) signaling by inhibiting AR expression, binding to and inhibiting AR activity at gene-specific loci, and inducing repressive epigenetic programs via direct activation of the H3K27 methyltransferase EZH2, a Polycomb group protein. These findings provide a working model in which TMPRSS2-ERG plays a critical role in cancer progression by disrupting lineage-specific differentiation of the prostate and potentiating the EZH2-mediated dedifferentiation program.

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.

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.

Genomic loss of microRNA-101 leads to overexpression of histone methyltransferase EZH2 in cancer.

Enhancer of zeste homolog 2 (EZH2) is a mammalian histone methyltransferase that contributes to the epigenetic silencing of target genes and regulates the survival and metastasis of cancer cells. EZH2 is overexpressed in aggressive solid tumors by mechanisms that remain unclear. Here we show that the expression and function of EZH2 in cancer cell lines are inhibited by microRNA-101 (miR-101). Analysis of human prostate tumors revealed that miR-101 expression decreases during cancer progression, paralleling an increase in EZH2 expression. One or both of the two genomic loci encoding miR-101 were somatically lost in 37.5% of clinically localized prostate cancer cells (6 of 16) and 66.7% of metastatic disease cells (22 of 33). We propose that the genomic loss of miR-101 in cancer leads to overexpression of EZH2 and concomitant dysregulation of epigenetic pathways, resulting in cancer progression.

Golgi protein GOLM1 is a tissue and urine biomarker of prostate cancer.

Prostate cancer is the most common type of tumor found in American men and is the second leading cause of cancer death in males. To identify biomarkers that distinguish prostate cancer from normal, we compared multiple gene expression profiling studies. Through meta-analysis of expression array data from multiple prostate cancer studies, we identified GOLM1 (Golgi membrane protein 1, Golm 1) as consistently up-regulated in clinically localized prostate cancer. This observation was confirmed by reverse transcription-polymerase chain reaction (RT-PCR) and validated at the protein level by immunoblot assay and immunohistochemistry. Prostate epithelial cells were identified as the cellular source of GOLM1 expression using laser capture microdissection. Immunohistochemical staining localized the GOLM1 signal to the subapical cytoplasmic region, typical of a Golgi distribution. Surprisingly, GOLM1 immunoreactivity was detected in the supernatants of prostate cell lines and in the urine of patients with prostate cancer. The mechanism by which intact GOLM1 might be released from cells has not yet been elucidated. GOLM1 transcript levels were measured in urine sediments using quantitative PCR on a cohort of patients presenting for biopsy or radical prostatectomy. We found that urinary GOLM1 mRNA levels were a significant predictor of prostate cancer. Further, GOLM1 outperformed serum prostate-specific antigen (PSA) in detecting prostate cancer. The area under the receiver-operating characteristic curve was 0.622 for GOLM1 (P = .0009) versus 0.495 for serum PSA (P = .902). Our data indicating the up-regulation of GOLM1 expression and its appearance in patients' urine suggest GOLM1 as a potential novel biomarker for clinically localized prostate cancer.

A fluorescence in situ hybridization screen for E26 transformation-specific aberrations: identification of DDX5-ETV4 fusion protein in prostate cancer.

Recurrent gene fusions involving E26 transformation-specific (ETS) transcription factors ERG, ETV1, ETV4, or ETV5 have been identified in 40% to 70% of prostate cancers. Here, we used a comprehensive fluorescence in situ hybridization (FISH) split probe strategy interrogating all 27 ETS family members and their five known 5' fusion partners in a cohort of 110 clinically localized prostate cancer patients. Gene rearrangements were only identified in ETS genes that were previously implicated in prostate cancer gene fusions including ERG, ETV1, and ETV4 (43%, 5%, and 5%, respectively), suggesting that a substantial fraction of prostate cancers (estimated at 30-60%) cannot be attributed to an ETS gene fusion. Among the known 5' gene fusion partners, TMPRSS2 was rearranged in 47% of cases followed by SLC45A3, HNRPA2B1, and C15ORF21 in 2%, 1%, and 1% of cases, respectively. Based on this comprehensive FISH screen, we have made four noteworthy observations. First, by screening the entire ETS transcription factor family for rearrangements, we found that a large fraction of prostate cancers (44%) cannot be ascribed to an ETS gene fusion, an observation which will stimulate research into identifying recurrent non-ETS aberrations in prostate cancers. Second, we identified SLC45A3 as a novel 5' fusion partner of ERG; previously, TMPRSS2 was the only described 5' partner of ERG. Third, we identified two prostate-specific, androgen-induced genes, FLJ35294 and CANT1, as 5' partners to ETV1 and ETV4. Fourth, we identified a ubiquitously expressed, androgen-insensitive gene, DDX5, fused in frame with ETV4, leading to the expression of a DDX5-ETV4 fusion protein.

The role of SPINK1 in ETS rearrangement-negative prostate cancers.

ETS gene fusions have been characterized in a majority of prostate cancers; however, the key molecular alterations in ETS-negative cancers are unclear. Here we used an outlier meta-analysis (meta-COPA) to identify SPINK1 outlier expression exclusively in a subset of ETS rearrangement-negative cancers ( approximately 10% of total cases). We validated the mutual exclusivity of SPINK1 expression and ETS fusion status, demonstrated that SPINK1 outlier expression can be detected noninvasively in urine, and observed that SPINK1 outlier expression is an independent predictor of biochemical recurrence after resection. We identified the aggressive 22RV1 cell line as a SPINK1 outlier expression model and demonstrate that SPINK1 knockdown in 22RV1 attenuates invasion, suggesting a functional role in ETS rearrangement-negative prostate cancers.

A first-generation multiplex biomarker analysis of urine for the early detection of prostate cancer.

Although prostate-specific antigen (PSA) serum level is currently the standard of care for prostate cancer screening in the United States, it lacks ideal specificity and additional biomarkers are needed to supplement or potentially replace serum PSA testing. Emerging evidence suggests that monitoring the noncoding RNA transcript PCA3 in urine may be useful in detecting prostate cancer in patients with elevated PSA levels. Here, we show that a multiplex panel of urine transcripts outperforms PCA3 transcript alone for the detection of prostate cancer. We measured the expression of seven putative prostate cancer biomarkers, including PCA3, in sedimented urine using quantitative PCR on a cohort of 234 patients presenting for biopsy or radical prostatectomy. By univariate analysis, we found that increased GOLPH2, SPINK1, and PCA3 transcript expression and TMPRSS2:ERG fusion status were significant predictors of prostate cancer. Multivariate regression analysis showed that a multiplexed model, including these biomarkers, outperformed serum PSA or PCA3 alone in detecting prostate cancer. The area under the receiver-operating characteristic curve was 0.758 for the multiplexed model versus 0.662 for PCA3 alone (P = 0.003). The sensitivity and specificity for the multiplexed model were 65.9% and 76.0%, respectively, and the positive and negative predictive values were 79.8% and 60.8%, respectively. Taken together, these results provide the framework for the development of highly optimized, multiplex urine biomarker tests for more accurate detection of prostate cancer.

Role of the TMPRSS2-ERG gene fusion in prostate cancer.

TMPRSS2-ERG gene fusions are the predominant molecular subtype of prostate cancer. Here, we explored the role of TMPRSS2-ERG gene fusion product using in vitro and in vivo model systems. Transgenic mice expressing the ERG gene fusion product under androgen-regulation develop mouse prostatic intraepithelial neoplasia (PIN), a precursor lesion of prostate cancer. Introduction of the ERG gene fusion product into primary or immortalized benign prostate epithelial cells induced an invasion-associated transcriptional program but did not increase cellular proliferation or anchorage-independent growth. These results suggest that TMPRSS2-ERG may not be sufficient for transformation in the absence of secondary molecular lesions. Transcriptional profiling of ERG knockdown in the TMPPRSS2-ERG-positive prostate cancer cell line VCaP revealed decreased expression of genes over-expressed in prostate cancer versus PIN and genes overexpressed in ETS-positive versus -negative prostate cancers in addition to inhibiting invasion. ERG knockdown in VCaP cells also induced a transcriptional program consistent with prostate differentiation. Importantly, VCaP cells and benign prostate cells overexpressing ERG directly engage components of the plasminogen activation pathway to mediate cellular invasion, potentially representing a downstream ETS target susceptible to therapeutic intervention. Our results support previous work suggesting that TMPRSS2-ERG fusions mediate invasion, consistent with the defining histologic distinction between PIN and prostate cancer.