Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

Abstract 1905: Defining a therapeutic classification of breast cancer by actionable targets

Breast cancer remains a frequent cause of female cancer death despite the great strides in elucidation of biological subtypes and their reported clinical and prognostic significance. We have defined a general cohort of breast cancers in terms of putative actionable targets, involving growth and proliferative factors, the cell cycle, and apoptotic pathways, both as single biomarkers across a general cohort and within intrinsic molecular subtypes.

We identified 293 patients treated with adjuvant chemotherapy. Additional hormonal therapy and trastuzumab was administered depending on hormonal and HER2 status respectively. We performed immunohistochemistry for ER, PR, HER2, MM1, CK5/6, p53, TOP2A, EGFR, IGF1R, PTEN, p-mTOR and e-cadherin. The cohort was classified into luminal (62%) and non-luminal (38%) tumors as well as luminal A (27%), luminal B HER2 negative (22%) and positive (12%), HER2 enriched (14%) and triple negative (25%). Patients with luminal tumors and co-overexpression of TOP2A or IGF1R loss displayed worse overall survival (p=0.0251 and p=0.0008 respectively). Non-luminal tumors had much greater heterogeneous expression profiles with no individual markers of prognostic significance. Non-luminal tumors were characterised by EGFR and TOP2A overexpression, IGF1R, PTEN and p-mTOR negativity and extreme p53 expression.

Our results indicate that only a minority of intrinsic subtype tumors purely express single novel actionable targets. This lack of pure biomarker expression is particular prevalent in the triple negative subgroup and may allude to the mechanism of targeted therapy inaction and myriad disappointing trial results. Utilising a combinatorial biomarker approach may enhance studies of targeted therapies providing additional information during design and patient selection while also helping decipher negative trial results.

The prognostic value of the stem-like group in colorectal cancer using a panel of immunohistochemistry markers

Colorectal cancer (CRC) is the second leading cause of cancer-related deaths in the Western world. It is becoming increasingly clear that CRC is a diverse disease, as exemplified by the identification of subgroups of CRC tumours that are driven by distinct biology. Recently, a number of studies have begun to define panels of diagnostically relevant markers to align patients into individual subgroups in an attempt to give information on prognosis and treatment response. We examined the immunohistochemical expression profile of 18 markers, each representing a putative role in cancer development, in 493 primary colorectal carcinomas using tissue microarrays. Through unsupervised clustering in stage II cancers, we identified two cluster groups that are broadly defined by inflammatory or immune-related factors (CD3, CD8, COX-2 and FOXP3) and stem-like factors (CD44, LGR5, SOX2, OCT4). The expression of the stem-like group markers was associated with a significantly worse prognosis compared to cases with lower expression. In addition, patients classified in the stem-like subgroup displayed a trend towards a benefit from adjuvant treatment. The biologically relevant and poor prognostic stem-like group could also be identified in early stage I cancers, suggesting a potential opportunity for the identification of aggressive tumors at a very early stage of the disease.

Protein expression differences between lung adenocarcinoma and squamous cell carcinoma with brain metastasis

AIM: We investigated tissue biomarkers in non-small cell lung cancer (NSCLC) to find indicators of brain metastasis and peritumoral brain edema.

PATIENTS AND METHODS: Fifty-two cases were studied out of which 26 had corresponding brain metastatic tissue. Clinicopathological characteristics of tumors were correlated with biomarkers of cell adhesion, cell growth, cell cycle and apoptosis regulation that were previously immunohistochemically studied but never analyzed separately according to histological subgroups, gender and smoking history.

RESULTS: Increased collagen XVII in adenocarcinoma (ADC) and increased caspase-9, CD44v6, and decreased cellular apoptosis susceptibility protein (CAS) and Ki-67 in squamous cell carcinoma (SCC) correlated significantly with brain metastasis. Increased β-catenin, E-cadherin and decreased caspase-9 expression in primary SCC, and decreased CD44v6 expression in brain metastatic SCC tissues showed a significant correlation with the extent of peritumoral brain edema. Positive correlation between smoking and biomarker expression could be observed in metastatic ADCs with p16 and caspase-8, while-negative correlation was found in SCC without brain metastasis with caspase-3, and in SCC with brain metastasis with p27.

CONCLUSION: Our results highlight the importance of separate analysis of biomarker expression in histological subtypes of NSCLC.

Tumorigenesis and peritoneal colonization from fallopian tube epithelium

Ovarian cancer is the most lethal gynecological malignancy, primarily because its origin and initiation factors are unknown. A secretory murine oviductal epithelial (MOE) model was generated to address the hypothesis that the fallopian tube is an origin for high-grade serous cancer. MOE cells were stably altered to express mutation in p53, silence PTEN, activate AKT, and amplify KRAS alone and in combination, to define if this cell type gives rise to tumors and what genetic alterations are required to drive malignancy. Cell lines were characterized in vitro and allografted into mice. Silencing PTEN formed high-grade carcinoma with wide spread tumor explants including metastasis into the ovary. Addition of p53 mutation to PTEN silencing did not enhance this phenotype, whereas addition of KRAS mutation reduced survival. Interestingly, PTEN silencing and KRAS mutation originating from ovarian surface epithelium generated endometrioid carcinoma, suggesting that different cellular origins with identical genetic manipulations can give rise to distinct cancer histotypes. Defining the roles of specific signaling modifications in tumorigenesis from the fallopian tube/oviduct is essential for early detection and development of targeted therapeutics. Further, syngeneic MOE allografts provide an ideal model for pre-clinical testing in an in vivo environment with an intact immune system.

Inter- and Intrafraction Target Motion in Highly Focused Single Vocal Cord Irradiation of T1a Larynx Cancer Patients

PURPOSE: The purpose of this study was to verify clinical target volume-planning target volume (CTV-PTV) margins in single vocal cord irradiation (SVCI) of T1a larynx tumors and characterize inter- and intrafraction target motion.

METHODS AND MATERIALS: For 42 patients, a single vocal cord was irradiated using intensity modulated radiation therapy at a total dose of 58.1 Gy (16 fractions × 3.63 Gy). A daily cone beam computed tomography (CBCT) scan was performed to online correct the setup of the thyroid cartilage after patient positioning with in-room lasers (interfraction motion correction). To monitor intrafraction motion, CBCT scans were also acquired just after patient repositioning and after dose delivery. A mixed online-offline setup correction protocol ("O2 protocol") was designed to compensate for both inter- and intrafraction motion.

RESULTS: Observed interfraction, systematic (Σ), and random (σ) setup errors in left-right (LR), craniocaudal (CC), and anteroposterior (AP) directions were 0.9, 2.0, and 1.1 mm and 1.0, 1.6, and 1.0 mm, respectively. After correction of these errors, the following intrafraction movements derived from the CBCT acquired after dose delivery were: Σ = 0.4, 1.3, and 0.7 mm, and σ = 0.8, 1.4, and 0.8 mm. More than half of the patients showed a systematic non-zero intrafraction shift in target position, (ie, the mean intrafraction displacement over the treatment fractions was statistically significantly different from zero; P<.05). With the applied CTV-PTV margins (for most patients 3, 5, and 3 mm in LR, CC, and AP directions, respectively), the minimum CTV dose, estimated from the target displacements observed in the last CBCT, was at least 94% of the prescribed dose for all patients and more than 98% for most patients (37 of 42). The proposed O2 protocol could effectively reduce the systematic intrafraction errors observed after dose delivery to almost zero (Σ = 0.1, 0.2, 0.2 mm).

CONCLUSIONS: With adequate image guidance and CTV-PTV margins in LR, CC, and AP directions of 3, 5, and 3 mm, respectively, excellent target coverage in SVCI could be ensured.

Infrequent mutations of Archipelago (hAGO, hCDC4, Fbw7) in primary ovarian cancer

Objective: Archipelago (AGO, also known as hCdc4, Fbw7, or Sel-10) is an F-box containing component of the SCF complex implicated in the ubiquitination and proteolysis of cyclin E and c-Myc, and found to be mutated in 16% of endometrial carcinomas. We have previously reported somatic mutations in AGO in 3/10 ovarian cancer cell lines, but the frequency of such mutations in primary ovarian cancer is unknown.

Methods: The coding sequence of AGO was analyzed in 95 primary sporadic ovarian tumors and 16 cases of familial ovarian cancer, and correlated with levels of cyclin E and c-Myc protein expression. Constructs encoding mutations in AGO were transfected into an AGO-null cell line to directly test their ability to regulate cyclin E and c-Myc levels.

Results: Mutations were present in only 2 of 95 sporadic cases: a premature stop within the WD domain (471 Ter) and a missense change near the F-box (S245T). Both primary tumor specimens containing these mutations showed high levels of cyclin E and c-Myc, but reconstitution of an AGO-null cell line with constructs encoding these mutations showed 471 Ter to be inactive in regulating endogenous cyclin E and c-Myc levels, while the S245T mutant was indistinguishable from wild-type. No germ-line mutations were found in familial cases of ovarian cancer.

Conclusion: Somatic AGO mutations are infrequent in primary ovarian cancers and are unlikely to contribute to familial ovarian cancer. Reconstitution experiments, rather than measuring tumor levels of cyclin E and c-Myc, provide an effective approach to determine the functional significance of AGO mutations identified in human cancers.

PTEN mRNA detection by chromogenic, RNA in situ technologies: a reliable alternative to PTEN immunohistochemistry

Immunohistochemical staining for phosphatase and tensin homolog (PTEN) does not have either an acceptable standard protocol or concordance of scoring between pathologists. Evaluation of PTEN mRNA with a unique and verified sequence probe may offer a realistic alternative providing a robust and reproducible protocol. In this study, we have evaluated an in situ hybridization (ISH) protocol for PTEN mRNA using RNAScope technology and compared it with a standard protocol for PTEN immunohistochemistry (IHC). PTEN mRNA expression by ISH was consistently more sensitive than PTEN IHC, with 56% of samples on a mixed-tumor tissue microarray (TMA) showing high expression by ISH compared with 42% by IHC. On a prostate TMA, 49% of cases showed high expression by ISH compared with 43% by IHC. Variations in PTEN mRNA expression within malignant epithelium were quantifiable using image analysis on the prostate TMAs. Within tumors, clear overexpression of PTEN mRNA on malignant epithelium compared with benign epithelium was frequently observed and quantified. The use of SpotStudio software in the mixed-tumor TMA allowed for clear demonstration of varying levels of PTEN mRNA between tumor samples by the mRNA methodology. This was evident by the quantifiable differences between distinct oropharyngeal tumors (up to 3-fold increase in average number of spots per cell between 2 cases). mRNA detection of PTEN or other biomarkers, for which optimal or standardized immunohistochemical techniques are not available, represents a means by which heterogeneity of expression within focal regions of tumor can be explored with more confidence.

Gene expression and epigenetic discovery screen reveal methylation of SFRP2 in prostate cancer

Aberrant activation of Wnts is common in human cancers, including prostate. Hypermethylation associated transcriptional silencing of Wnt antagonist genes SFRPs (Secreted Frizzled-Related Proteins) is a frequent oncogenic event. The significance of this is not known in prostate cancer. The objectives of our study were to (i) profile Wnt signaling related gene expression and (ii) investigate methylation of Wnt antagonist genes in prostate cancer. Using TaqMan Low Density Arrays, we identified 15 Wnt signaling related genes with significantly altered expression in prostate cancer; the majority of which were upregulated in tumors. Notably, histologically benign tissue from men with prostate cancer appeared more similar to tumor (r = 0.76) than to benign prostatic hyperplasia (BPH; r = 0.57, p < 0.001). Overall, the expression profile was highly similar between tumors of high (≥ 7) and low (≤ 6) Gleason scores. Pharmacological demethylation of PC-3 cells with 5-Aza-CdR reactivated 39 genes (≥ 2-fold); 40% of which inhibit Wnt signaling. Methylation frequencies in prostate cancer were 10% (2/20) (SFRP1), 64.86% (48/74) (SFRP2), 0% (0/20) (SFRP4) and 60% (12/20) (SFRP5). SFRP2 methylation was detected at significantly lower frequencies in high-grade prostatic intraepithelial neoplasia (HGPIN; 30%, (6/20), p = 0.0096), tumor adjacent benign areas (8.82%, (7/69), p < 0.0001) and BPH (11.43% (4/35), p < 0.0001). The quantitative level of SFRP2 methylation (normalized index of methylation) was also significantly higher in tumors (116) than in the other samples (HGPIN = 7.45, HB = 0.47, and BPH = 0.12). We show that SFRP2 hypermethylation is a common event in prostate cancer. SFRP2 methylation in combination with other epigenetic markers may be a useful biomarker of prostate cancer.

In silico analysis and DHPLC screening strategy identifies novel apoptotic gene targets of aberrant promoter hypermethylation in prostate cancer

BACKGROUND: Aberrant DNA methylation has been implicated as a key survival mechanism in cancer, whereby promoter hypermethylation silences genes essential for many cellular processes including apoptosis. Limited data is available on the methylation profile of apoptotic genes in prostate cancer (CaP). The aim of this study was to profile methylation of apoptotic-related genes in CaP using denaturing high performance liquid chromatography (DHPLC).

METHODS: Based on an in silico selection process, 13 genes were screened for methylation in CaP cell lines using DHPLC. Quantitative methylation specific PCR was employed to determine methylation levels in prostate tissue specimens (n = 135), representing tumor, histologically benign prostate, high-grade prostatic intraepithelial neoplasia and benign prostatic hyperplasia. Gene expression was measured by QRT-PCR in cell lines and tissue specimens.

RESULTS: The promoters of BIK, BNIP3, cFLIP, TMS1, DCR1, DCR2, and CDKN2A appeared fully or partially methylated in a number of malignant cell lines. This is the first report of aberrant methylation of BIK, BNIP3, and cFLIP in CaP. Quantitative methylation analysis in prostate tissues identified 5 genes (BNIP3, CDKN2A, DCR1, DCR2 and TMS1) which were frequently methylated in tumors but were unmethylated in 100% of benign tissues. Furthermore, 69% of tumors were methylated in at least one of the five-gene panel. In the case of all genes, except BNIP3, promoter hypermethylation was associated with concurrent downregulation of gene expression.

CONCLUSION: Future examination of this "CaP apoptotic methylation signature" in a larger cohort of patients is justified to further evaluate its value as a diagnostic and prognostic marker.