Pan-cancer analysis of mRNA stability for decoding tumour post-transcriptional programs.

Measuring mRNA decay in tumours is a prohibitive challenge, limiting our ability to map the post-transcriptional programs of cancer. Here, using a statistical framework to decouple transcriptional and post-transcriptional effects in RNA-seq data, we uncover the mRNA stability changes that accompany tumour development and progression. Analysis of 7760 samples across 18 cancer types suggests that mRNA stability changes are ~30% as frequent as transcriptional events, highlighting their widespread role in shaping the tumour transcriptome. Dysregulation of programs associated with >80 RNA-binding proteins (RBPs) and microRNAs (miRNAs) drive these changes, including multi-cancer inactivation of RBFOX and miR-29 families. Phenotypic activation or inhibition of RBFOX1 highlights its role in calcium signaling dysregulation, while modulation of miR-29 shows its impact on extracellular matrix organization and stemness genes. Overall, our study underlines the integral role of mRNA stability in shaping the cancer transcriptome, and provides a resource for systematic interrogation of cancer-associated stability pathways.

Alkaloid production and response to natural adverse conditions in Peganum harmala: in silico transcriptome analyses.

Peganum harmala is a valuable wild plant that grows and survives under adverse conditions and produces pharmaceutical alkaloid metabolites. Using different assemblers to develop a transcriptome improves the quality of assembled transcriptome. In this study, a concrete and accurate method for detecting stress-responsive transcripts by comparing stress-related gene ontology (GO) terms and public domains was designed. An integrated transcriptome for P. harmala including 42 656 coding sequences was created by merging de novo assembled transcriptomes. Around 35 000 transcripts were annotated with more than 90% resemblance to three closely related species of Citrus, which confirmed the robustness of the assembled transcriptome; 4853 stress-responsive transcripts were identified. CYP82 involved in alkaloid biosynthesis showed a higher number of transcripts in P. harmala than in other plants, indicating its diverse alkaloid biosynthesis attributes. Transcription factors (TFs) and regulatory elements with 3887 transcripts comprised 9% of the transcriptome. Among the TFs of the integrated transcriptome, cystein2/histidine2 (C2H2) and WD40 repeat families were the most abundant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) MAPK (mitogen-activated protein kinase) signaling map and the plant hormone signal transduction map showed the highest assigned genes to these pathways, suggesting their potential stress resistance. The P. harmala whole-transcriptome survey provides important resources and paves the way for functional and comparative genomic studies on this plant to discover stress-tolerance-related markers and response mechanisms in stress physiology, phytochemistry, ecology, biodiversity, and evolution. P. harmala can be a potential model for studying adverse environmental cues and metabolite biosynthesis and a major source for the production of various alkaloids.

Rational Development of Liquid Biopsy Analysis in Renal Cell Carcinoma.

Renal cell carcinoma (RCC) is known for its variable clinical behavior and outcome, including heterogeneity in developing relapse or metastasis. Recent data highlighted the potential of somatic mutations as promising biomarkers for risk stratification in RCC. Likewise, the analysis of circulating tumor DNA (ctDNA) for such informative somatic mutations (liquid biopsy) is considered an important advance for precision oncology in RCC, allowing to monitor molecular disease evolution in real time. However, our knowledge about the utility of ctDNA analysis in RCC is limited, in part due to the lack of RCC-appropriate assays for ctDNA analysis. Here, by interrogating different blood compartments in xenograft models, we identified plasma cell-free (cf) DNA and extracellular vesicles (ev) DNA enriched for RCC-associated ctDNA. Additionally, we developed sensitive targeted sequencing and bioinformatics workflows capable of detecting somatic mutations in RCC-relevant genes with allele frequencies ≥ 0.5%. Applying this assay to patient-matched tumor and liquid biopsies, we captured tumor mutations in cf- and ev-DNA fractions isolated from the blood, highlighting the potentials of both fractions for ctDNA analysis. Overall, our study presents an RCC-appropriate sequencing assay and workflow for ctDNA analysis and provides a proof of principle as to the feasibility of detecting tumor-specific mutations in liquid biopsy in RCC patients.

Targeting lysyl oxidase (LOX) overcomes chemotherapy resistance in triple negative breast cancer.

Chemoresistance is a major obstacle in triple negative breast cancer (TNBC), the most aggressive breast cancer subtype. Here we identify hypoxia-induced ECM re-modeler, lysyl oxidase (LOX) as a key inducer of chemoresistance by developing chemoresistant TNBC tumors in vivo and characterizing their transcriptomes by RNA-sequencing. Inhibiting LOX reduces collagen cross-linking and fibronectin assembly, increases drug penetration, and downregulates ITGA5/FN1 expression, resulting in inhibition of FAK/Src signaling, induction of apoptosis and re-sensitization to chemotherapy. Similarly, inhibiting FAK/Src results in chemosensitization. These effects are observed in 3D-cultured cell lines, tumor organoids, chemoresistant xenografts, syngeneic tumors and PDX models. Re-expressing the hypoxia-repressed miR-142-3p, which targets HIF1A, LOX and ITGA5, causes further suppression of the HIF-1α/LOX/ITGA5/FN1 axis. Notably, higher LOX, ITGA5, or FN1, or lower miR-142-3p levels are associated with shorter survival in chemotherapy-treated TNBC patients. These results provide strong pre-clinical rationale for developing and testing LOX inhibitors to overcome chemoresistance in TNBC patients.

Unfolded protein response-mediated modulation of mesenchymal stem cells.

The endoplasmic reticulum (ER) receives unfolded proteins predestined for the secretory pathway or to be incorporated as transmembrane proteins. The ER has to accommodate the proper folding and glycosylation of these proteins and also to properly incorporate transmembrane proteins. However, under various circumstances, the proteins shuttling through the ER can be misfolded and undergo aggregation, which causes activation of the unfolded protein response (UPR). The UPR is mediated through three primary pathways: activating transcription factor-6, inositol-requiring enzyme-1 (IRE1), and PKR-like endoplasmic reticulum kinase, which up-regulate ER folding chaperones and temporarily suppress protein translation. The UPR can be both cytoprotective and/or cytotoxic depending on the duration of UPR activation and the type of host cell. Proteostasis controls stem cell function, while stress responses affect stem cell identity and differentiation. The present review aimed to explore and discuss the effects of the UPR pathways on mesenchymal stem cells.

GHSR DNA hypermethylation is a new epigenetic biomarker for gastric adenocarcinoma and beyond.

Aberrations of DNA methylation are early events in the development of tumors. In this study, we investigated the DNA methylation status of growth hormone secretagogue receptor (GHSR), a promising pan-cancer biomarker, in gastric cancer (GC). Initially, data sets from DNA methylation and gene expression studies available at Gene Expression Omnibus (GEO) were analyzed. Confirmation was done on primary tumor specimens and adjacent normal stomach tissue samples. Both analyses showed significant hypermethylation of GHSR. For further validation, The Cancer Genome Atlas data on stomach cancer was used. A receiver operating characteristic curve analysis yielded an area under the curve value of 0.85, corroborating its usefulness as a diagnostic marker. A genome-wide comethylation analysis revealed several correlated genes. CREB1 was found to act as an upstream regulator of this gene network. Furthermore, GHSR methylation was found to be a biomarker in several other tumor entities, namely cancers of the bladder, endometrium, esophagus, head and neck, liver, thyroid, kidney, and ovary. Our findings along with previous reports on other types of cancer suggest a high potential of GHSR gene methylation as a pan-cancer biomarker, which could be considered for liquid biopsy applications.

Aberration hubs in protein interaction networks highlight actionable targets in cancer.

Despite efforts for extensive molecular characterization of cancer patients, such as the international cancer genome consortium (ICGC) and the cancer genome atlas (TCGA), the heterogeneous nature of cancer and our limited knowledge of the contextual function of proteins have complicated the identification of targetable genes. Here, we present Aberration Hub Analysis for Cancer (AbHAC) as a novel integrative approach to pinpoint aberration hubs, i.e. individual proteins that interact extensively with genes that show aberrant mutation or expression. Our analysis of the breast cancer data of the TCGA and the renal cancer data from the ICGC shows that aberration hubs are involved in relevant cancer pathways, including factors promoting cell cycle and DNA replication in basal-like breast tumors, and Src kinase and VEGF signaling in renal carcinoma. Moreover, our analysis uncovers novel functionally relevant and actionable targets, among which we have experimentally validated abnormal splicing of spleen tyrosine kinase as a key factor for cell proliferation in renal cancer. Thus, AbHAC provides an effective strategy to uncover novel disease factors that are only identifiable by examining mutational and expression data in the context of biological networks.

A General Framework for Interrogation of mRNA Stability Programs Identifies RNA-Binding Proteins that Govern Cancer Transcriptomes.

Widespread remodeling of the transcriptome is a signature of cancer; however, little is known about the post-transcriptional regulatory factors, including RNA-binding proteins (RBPs) that regulate mRNA stability, and the extent to which RBPs contribute to cancer-associated pathways. Here, by modeling the global change in gene expression based on the effect of sequence-specific RBPs on mRNA stability, we show that RBP-mediated stability programs are recurrently deregulated in cancerous tissues. Particularly, we uncovered several RBPs that contribute to the abnormal transcriptome of renal cell carcinoma (RCC), including PCBP2, ESRP2, and MBNL2. Modulation of these proteins in cancer cell lines alters the expression of pathways that are central to the disease and highlights RBPs as driving master regulators of RCC transcriptome. This study presents a framework for the screening of RBP activities based on computational modeling of mRNA stability programs in cancer and highlights the role of post-transcriptional gene dysregulation in RCC.

Reactivation of cAMP Pathway by PDE4D Inhibition Represents a Novel Druggable Axis for Overcoming Tamoxifen Resistance in ER-positive Breast Cancer.

Purpose: Tamoxifen remains an important hormonal therapy for ER-positive breast cancer; however, development of resistance is a major obstacle in clinics. Here, we aimed to identify novel mechanisms of tamoxifen resistance and provide actionable drug targets overcoming resistance.Experimental Design: Whole-transcriptome sequencing, downstream pathway analysis, and drug repositioning approaches were used to identify novel modulators [here: phosphodiesterase 4D (PDE4D)] of tamoxifen resistance. Clinical data involving tamoxifen-treated patients with ER-positive breast cancer were used to assess the impact of PDE4D in tamoxifen resistance. Tamoxifen sensitization role of PDE4D was tested in vitro and in vivo Cytobiology, biochemistry, and functional genomics tools were used to elucidate the mechanisms of PDE4D-mediated tamoxifen resistance.Results: PDE4D, which hydrolyzes cyclic AMP (cAMP), was significantly overexpressed in both MCF-7 and T47D tamoxifen-resistant (TamR) cells. Higher PDE4D expression predicted worse survival in tamoxifen-treated patients with breast cancer (n = 469, P = 0.0036 for DMFS; n = 561, P = 0.0229 for RFS) and remained an independent prognostic factor for RFS in multivariate analysis (n = 132, P = 0.049). Inhibition of PDE4D by either siRNAs or pharmacologic inhibitors (dipyridamole and Gebr-7b) restored tamoxifen sensitivity. Sensitization to tamoxifen is achieved via cAMP-mediated induction of unfolded protein response/ER stress pathway leading to activation of p38/JNK signaling and apoptosis. Remarkably, acetylsalicylic acid (aspirin) was predicted to be a tamoxifen sensitizer using a drug repositioning approach and was shown to reverse resistance by targeting PDE4D/cAMP/ER stress axis. Finally, combining PDE4D inhibitors and tamoxifen suppressed tumor growth better than individual groups in vivoConclusions: PDE4D plays a pivotal role in acquired tamoxifen resistance via blocking cAMP/ER stress/p38-JNK signaling and apoptosis. Clin Cancer Res; 24(8); 1987-2001. ©2018 AACR.

Loss of chromosome Y leads to down regulation of KDM5D and KDM6C epigenetic modifiers in clear cell renal cell carcinoma.

Recent genomic studies of sporadic clear cell renal cell carcinoma (ccRCC) have uncovered novel driver genes and pathways. Given the unequal incidence rates among men and women (male:female incidence ratio approaches 2:1), we compared the genome-wide distribution of the chromosomal abnormalities in both sexes. We observed a higher frequency for the somatic recurrent chromosomal copy number variations (CNVs) of autosomes in male subjects, whereas somatic loss of chromosome X was detected exclusively in female patients (17.1%). Furthermore, somatic loss of chromosome Y (LOY) was detected in about 40% of male subjects, while mosaic LOY was detected in DNA isolated from peripheral blood in 9.6% of them, and was the only recurrent CNV in constitutional DNA samples. LOY in constitutional DNA, but not in tumor DNA was associated with older age. Amongst Y-linked genes that were downregulated due to LOY, KDM5D and KDM6C epigenetic modifiers have functionally-similar X-linked homologs whose deficiency is involved in ccRCC progression. Our findings establish somatic LOY as a highly recurrent genetic defect in ccRCC that leads to downregulation of hitherto unsuspected epigenetic factors, and suggest that different mechanisms may underlie the somatic and mosaic LOY observed in tumors and peripheral blood, respectively.

Expression of DRD2 Is Increased in Human Pancreatic Ductal Adenocarcinoma and Inhibitors Slow Tumor Growth in Mice.

BACKGROUND & AIMS: Incidence of and mortality from pancreatic ductal adenocarcinoma (PDAC), the most common form of pancreatic cancer, are almost equivalent, so better treatments are needed. We studied gene expression profiles of PDACs and the functions of genes with altered expression to identify new therapeutic targets. METHODS: We performed microarray analysis to analyze gene expression profiles of 195 PDAC and 41 non-tumor pancreatic tissue samples. We undertook an extensive analysis of the PDAC transcriptome by superimposing interaction networks of proteins encoded by aberrantly expressed genes over signaling pathways associated with PDAC development to identify factors that might alter regulation of these pathways during tumor progression. We performed tissue microarray analysis to verify changes in expression of candidate protein using an independent set of 152 samples (40 nontumor pancreatic tissues, 63 PDAC sections, and 49 chronic pancreatitis samples). We validated the functional relevance of the candidate molecule using RNA interference or pharmacologic inhibitors in pancreatic cancer cell lines and analyses of xenograft tumors in mice. RESULTS: In an analysis of 38,276 human genes and loci, we identified 1676 genes that were significantly up-regulated and 1166 genes that were significantly down-regulated in PDAC compared with nontumor pancreatic tissues. One gene that was up-regulated and associated with multiple signaling pathways that are dysregulated in PDAC was G protein subunit αi2, which has not been previously associated with PDAC. G protein subunit αi2 mediates the effects of dopamine receptor D2 (DRD2) on cyclic adenosine monophosphate signaling; PDAC tissues had a slight but significant increase in DRD2 messenger RNA. Levels of DRD2 protein were substantially increased in PDACs, compared with non-tumor tissues, in tissue microarray analyses. RNA interference knockdown of DRD2 or inhibition with pharmacologic antagonists (pimozide and haloperidol) reduced proliferation of pancreatic cancer cells, induced endoplasmic reticulum stress and apoptosis, and reduced cell migration. RNA interference knockdown of DRD2 in pancreatic tumor cells reduced growth of xenograft tumors in mice, and administration of the DRD2 inhibitor haloperidol to mice with orthotopic xenograft tumors reduced final tumor size and metastasis. CONCLUSIONS: In gene expression profile analysis of PDAC samples, we found the DRD2 signaling pathway to be activated. Inhibition of DRD2 in pancreatic cancer cells reduced proliferation and migration, and slowed growth of xenograft tumors in mice. DRD2 antagonists routinely used for management of schizophrenia might be tested in patients with pancreatic cancer.

Early Epigenetic Downregulation of microRNA-192 Expression Promotes Pancreatic Cancer Progression.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by very early metastasis, suggesting the hypothesis that metastasis-associated changes may occur prior to actual tumor formation. In this study, we identified miR-192 as an epigenetically regulated suppressor gene with predictive value in this disease. miR-192 was downregulated by promoter methylation in both PDAC and chronic pancreatitis, the latter of which is a major risk factor for the development of PDAC. Functional studies in vitro and in vivo in mouse models of PDAC showed that overexpression of miR-192 was sufficient to reduce cell proliferation and invasion. Mechanistic analyses correlated changes in miR-192 promoter methylation and expression with epithelial-mesenchymal transition. Cell proliferation and invasion were linked to altered expression of the miR-192 target gene SERPINE1 that is encoding the protein plasminogen activator inhibitor-1 (PAI-1), an established regulator of these properties in PDAC cells. Notably, our data suggested that invasive capacity was altered even before neoplastic transformation occurred, as triggered by miR-192 downregulation. Overall, our results highlighted a role for miR-192 in explaining the early metastatic behavior of PDAC and suggested its relevance as a target to develop for early diagnostics and therapy. Cancer Res; 76(14); 4149-59. ©2016 AACR.

GHSR DNA hypermethylation is a common epigenetic alteration of high diagnostic value in a broad spectrum of cancers.

Identification of a single molecular trait that is determinant of common malignancies may serve as a powerful diagnostic supplement to cancer type-specific markers. Here, we report a DNA methylation mark that is characteristic of seven studied malignancies, namely cancers of lung, breast, prostate, pancreas, colorectum, glioblastoma and B cell chronic lymphocytic leukaemia (CLL) (n = 137). This mark was defined by substantial hypermethylation at the promoter and first exon of growth hormone secretagouge receptor (GHSR) through bisulfite pyrosequencing. The degree of aberrant methylation was capable of accurate discrimination between cancer and control samples. The highest sensitivity and specificity of cancer detection was achieved for cancers of pancreas, lung, breast and CLL yielding the area under the curve (AUC) values of 1.0000, 0.9952, 0.9800 and 0.9400, respectively. Narrowing to a single CpG site within the gene's promoter or four consecutive CpG units of the highest methylation levels within the first exon improved the detection power. GHSR hypermethylation was detected already at the early stage tumors. The accurate performance of this marker was further replicated in an independent set of pancreatic cancer and control samples (n = 78). These findings support the candidature of GHSR methylation as a highly accurate pan-cancer marker.

Diagnostic values of GHSR DNA methylation pattern in breast cancer.

DNA methylation patterns have been recognised as cancer-specific markers with high potential for clinical applications. We aimed at identifying methylation variations that differentiate between breast cancers and other breast tissue entities to establish a signature for diagnosis. Candidate genomic loci were analysed in 117 fresh-frozen breast specimens, which included cancer, benign and normal breast tissues from patients as well as material from healthy individuals. A cancer-specific DNA methylation signature was identified by microarray analysis in a test set of samples (n = 52, p < 2.1 × 10(-4)) and its performance was assessed through bisulphite pyrosequencing in an independent validation set (n = 65, p < 1.9 × 10(-7)). The signature is associated with SFRP2 and GHSR genes, and exhibited significant hypermethylation in cancers. Normal-appearing breast tissues from cancer patients were also methylated at these loci but to a markedly lower extent. This occurrence of methylated DNA in normal breast tissue of cancer patients is indicative of an epigenetic field defect. Concerning diagnosis, receiver operating characteristic curves and the corresponding area under the curve (AUC) analysis demonstrated a very high sensitivity and specificity of 89.3 and 100 %, respectively, for the GHSR methylation pattern (AUC >0.99). To date, this represents the DNA methylation marker of the highest sensitivity and specificity for breast cancer diagnosis. Functionally, ectopic expression of GHSR in a cell line model reduced breast cancer cell invasion without affecting cell viability upon stimulation of cells with ghrelin. Our data suggest a link between epigenetic down-regulation of GHSR and breast cancer cell invasion.

Correction of PCR-bias in quantitative DNA methylation studies by means of cubic polynomial regression.

DNA methylation profiling has become an important aspect of biomedical molecular analysis. Polymerase chain reaction (PCR) amplification of bisulphite-treated DNA is a processing step that is common to many currently used methods of quantitative methylation analysis. Preferential amplification of unmethylated alleles-known as PCR-bias-may significantly affect the accuracy of quantification. To date, no universal experimental approach has been reported to overcome the problem. This study presents an effective method of correcting biased methylation data. The procedure includes a calibration performed in parallel to the analysis of the samples under investigation. DNA samples with defined degrees of methylation are analysed. The observed deviation of the experimental results from the expected values is used for calculating a regression curve. The equation of the best-fitting curve is then used for correction of the data obtained from the samples of interest. The process can be applied irrespective of the locus interrogated and the number of sites analysed, avoiding an optimization of the amplification conditions for each individual locus.