CREBBP/EP300 bromodomains are critical to sustain the GATA1/MYC regulatory axis in proliferation.

BACKGROUND: The reported antitumor activity of the BET family bromodomain inhibitors has prompted the development of inhibitors against other bromodomains. However, the human genome encodes more than 60 different bromodomains and most of them remain unexplored. RESULTS: We report that the bromodomains of the histone acetyltransferases CREBBP/EP300 are critical to sustain the proliferation of human leukemia and lymphoma cell lines. EP300 is very abundant at super-enhancers in K562 and is coincident with sites of GATA1 and MYC occupancy. In accordance, CREBBP/EP300 bromodomain inhibitors interfere with GATA1- and MYC-driven transcription, causing the accumulation of cells in the G0/G1 phase of the cell cycle. The CREBBP/CBP30 bromodomain inhibitor CBP30 displaces CREBBP and EP300 from GATA1 and MYC binding sites at enhancers, resulting in a decrease in the levels of histone acetylation at these regulatory regions and consequently reduced gene expression of critical genes controlled by these transcription factors. CONCLUSIONS: Our data shows that inhibition of CREBBP/EP300 bromodomains can interfere with oncogene-driven transcriptional programs in cancer cells and consequently hold therapeutic potential.

A Braf kinase-inactive mutant induces lung adenocarcinoma.

The initiating oncogenic event in almost half of human lung adenocarcinomas is still unknown, a fact that complicates the development of selective targeted therapies. Yet these tumours harbour a number of alterations without obvious oncogenic function including BRAF-inactivating mutations. Inactivating BRAF mutants in lung predominate over the activating V600E mutant that is frequently observed in other tumour types. Here we demonstrate that the expression of an endogenous Braf(D631A) kinase-inactive isoform in mice (corresponding to the human BRAF(D594A) mutation) triggers lung adenocarcinoma in vivo, indicating that BRAF-inactivating mutations are initiating events in lung oncogenesis. Moreover, inactivating BRAF mutations have also been identified in a subset of KRAS-driven human lung tumours. Co-expression of Kras(G12V) and Braf(D631A) in mouse lung cells markedly enhances tumour initiation, a phenomenon mediated by Craf kinase activity, and effectively accelerates tumour progression when activated in advanced lung adenocarcinomas. We also report a key role for the wild-type Braf kinase in sustaining Kras(G12V)/Braf(D631A)-driven tumours. Ablation of the wild-type Braf allele prevents the development of lung adenocarcinoma by inducing a further increase in MAPK signalling that results in oncogenic toxicity; this effect can be abolished by pharmacological inhibition of Mek to restore tumour growth. However, the loss of wild-type Braf also induces transdifferentiation of club cells, which leads to the rapid development of lethal intrabronchiolar lesions. These observations indicate that the signal intensity of the MAPK pathway is a critical determinant not only in tumour development, but also in dictating the nature of the cancer-initiating cell and ultimately the resulting tumour phenotype.

Prediction of miRNA-mRNA Interactions Using miRGate.

miRGate ( http://mirgate.bioinfo.cnio.es /) is a freely available database that contains predicted and experimentally validated microRNA-messenger RNA (miRNA-mRNA) target pairs. This resource includes novel predictions from five well-established algorithms, but recalculated from a common and comprehensive sequence dataset. It includes all 3'-UTR sequences of all known genes of the three more widely employed genomes (human, mouse, and rat), and all annotated miRNA sequences from those genomes. Besides, it also contains predictions for all genes in human targeted by miRNA viruses such as Epstein-Barr and Kaposi sarcoma-associated herpes virus.The approach intends to circumvent one of the main drawbacks in this area, as diverse sequences and gene database versions cause poor overlap among different target prediction methods even with experimentally confirmed targets. As a result, miRGate predictions have been successfully validated using functional assays in several laboratories.This chapter describes how a user can access target information via miRGate's web interface. It also shows how automatically access the database through the programmatic interface based on representational state transfer services (REST), using the application programming interface (API) available at http://mirgate.bioinfo.cnio.es/API .

Vitamin D receptor expression and associated gene signature in tumour stromal fibroblasts predict clinical outcome in colorectal cancer.

OBJECTIVE: Colorectal cancer (CRC) is a major health concern. Vitamin D deficiency is associated with high CRC incidence and mortality, suggesting a protective effect of vitamin D against this disease. Given the strong influence of tumour stroma on cancer progression, we investigated the potential effects of the active vitamin D metabolite 1α,25-dihydroxyvitamin D3 (1,25(OH)2D3) on CRC stroma. DESIGN: Expression of vitamin D receptor (VDR) and two 1,25(OH)2D3 target genes was analysed in 658 patients with CRC with prolonged clinical follow-up. 1,25(OH)2D3 effects on primary cultures of patient-derived colon normal fibroblasts (NFs) and cancer-associated fibroblasts (CAFs) were studied using collagen gel contraction and migration assays and global gene expression analyses. Publicly available data sets (n=877) were used to correlate the 1,25(OH)2D3-associated gene signature in CAFs with CRC outcome. RESULTS: High VDR expression in tumour stromal fibroblasts was associated with better overall survival (OS) and progression-free survival in CRC, independently of its expression in carcinoma cells. 1,25(OH)2D3 inhibited the protumoural activation of NFs and CAFs and imposed in CAFs a 1,25(OH)2D3-associated gene signature that correlated with longer OS and disease-free survival in CRC. Furthermore, expression of two genes from the signature, CD82 and S100A4, correlated with stromal VDR expression and clinical outcome in our cohort of patients with CRC. CONCLUSIONS: 1,25(OH)2D3 has protective effects against CRC through the regulation of stromal fibroblasts. Accordingly, expression of VDR and 1,25(OH)2D3-associated gene signature in stromal fibroblasts predicts a favourable clinical outcome in CRC. Therefore, treatment of patients with CRC with VDR agonists could be explored even in the absence of VDR expression in carcinoma cells.

Deep sequencing reveals microRNAs predictive of antiangiogenic drug response.

The majority of metastatic renal cell carcinoma (RCC) patients are treated with tyrosine kinase inhibitors (TKI) in first-line treatment; however, a fraction are refractory to these antiangiogenic drugs. MicroRNAs (miRNAs) are regulatory molecules proven to be accurate biomarkers in cancer. Here, we identified miRNAs predictive of progressive disease under TKI treatment through deep sequencing of 74 metastatic clear cell RCC cases uniformly treated with these drugs. Twenty-nine miRNAs were differentially expressed in the tumors of patients who progressed under TKI therapy (P values from 6 × 10-9 to 3 × 10-3). Among 6 miRNAs selected for validation in an independent series, the most relevant associations corresponded to miR-1307-3p, miR-155-5p, and miR-221-3p (P = 4.6 × 10-3, 6.5 × 10-3, and 3.4 × 10-2, respectively). Furthermore, a 2 miRNA-based classifier discriminated individuals with progressive disease upon TKI treatment (AUC = 0.75, 95% CI, 0.64-0.85; P = 1.3 × 10-4) with better predictive value than clinicopathological risk factors commonly used. We also identified miRNAs significantly associated with progression-free survival and overall survival (P = 6.8 × 10-8 and 7.8 × 10-7 for top hits, respectively), and 7 overlapped with early progressive disease. In conclusion, this is the first miRNome comprehensive study, to our knowledge, that demonstrates a predictive value of miRNAs for TKI response and provides a new set of relevant markers that can help rationalize metastatic RCC treatment.

A genetic interaction between RAP1 and telomerase reveals an unanticipated role for RAP1 in telomere maintenance.

RAP1 is one of the components of shelterin, the capping complex at chromosome ends or telomeres, although its role in telomere length maintenance and protection has remained elusive. RAP1 also binds subtelomeric repeats and along chromosome arms, where it regulates gene expression and has been shown to function in metabolism control. Telomerase is the enzyme that elongates telomeres, and its deficiency causes a premature aging in humans and mice. We describe an unanticipated genetic interaction between RAP1 and telomerase. While RAP1 deficiency alone does not impact on mouse survival, mice lacking both RAP1 and telomerase show a progressively decreased survival with increasing mouse generations compared to telomerase single mutants. Telomere shortening is more pronounced in Rap1-/- Terc-/- doubly deficient mice than in the single-mutant Terc-/- counterparts, leading to an earlier onset of telomere-induced DNA damage and degenerative pathologies. Telomerase deficiency abolishes obesity and liver steatohepatitis provoked by RAP1 deficiency. Using genomewide ChIP sequencing, we find that progressive telomere shortening owing to telomerase deficiency leads to re-localization of RAP1 from telomeres and subtelomeric regions to extratelomeric sites in a genomewide manner. These findings suggest that although in the presence of sufficient telomere reserve RAP1 is not a key factor for telomere maintenance and protection, it plays a crucial role in the context of telomerase deficiency, thus in agreement with its evolutionary conservation as a telomere component from yeast to humans.

NSD2 contributes to oncogenic RAS-driven transcription in lung cancer cells through long-range epigenetic activation.

The histone methyltransferase NSD2/WHSC1/MMSET is overexpressed in a number of solid tumors but its contribution to the biology of these tumors is not well understood. Here, we describe that NSD2 contributes to the proliferation of a subset of lung cancer cell lines by supporting oncogenic RAS transcriptional responses. NSD2 knock down combined with MEK or BRD4 inhibitors causes co-operative inhibitory responses on cell growth. However, while MEK and BRD4 inhibitors converge in the downregulation of genes associated with cancer-acquired super-enhancers, NSD2 inhibition affects the expression of clusters of genes embedded in megabase-scale regions marked with H3K36me2 and that contribute to the RAS transcription program. Thus, combinatorial therapies using MEK or BRD4 inhibitors together with NSD2 inhibition are likely to be needed to ensure a more comprehensive inhibition of oncogenic RAS-driven transcription programs in lung cancers with NSD2 overexpression.

Analysis of Paired Primary-Metastatic Hormone-Receptor Positive Breast Tumors (HRPBC) Uncovers Potential Novel Drivers of Hormonal Resistance.

We sought to identify genetic variants associated with disease relapse and failure to hormonal treatment in hormone-receptor positive breast cancer (HRPBC). We analyzed a series of HRPBC with distant relapse, by sequencing pairs (n = 11) of tumors (primary and metastases) at >800X. Comparative genomic hybridization was performed as well. Top hits, based on the frequency of alteration and severity of the changes, were tested in the TCGA series. Genes determining the most parsimonious prognostic signature were studied for their functional role in vitro, by performing cell growth assays in hormonal-deprivation conditions, a setting that mimics treatment with aromatase inhibitors. Severe alterations were recurrently found in 18 genes in the pairs. However, only MYC, DNAH5, CSFR1, EPHA7, ARID1B, and KMT2C preserved an independent prognosis impact and/or showed a significantly different incidence of alterations between relapsed and non-relapsed cases in the TCGA series. The signature composed of MYC, KMT2C, and EPHA7 best discriminated the clinical course, (overall survival 90,7 vs. 144,5 months; p = 0.0001). Having an alteration in any of the genes of the signature implied a hazard ratio of death of 3.25 (p<0.0001), and early relapse during the adjuvant hormonal treatment. The presence of the D348N mutation in KMT2C and/or the T666I mutation in the kinase domain of EPHA7 conferred hormonal resistance in vitro. Novel inactivating mutations in KMT2C and EPHA7, which confer hormonal resistance, are linked to adverse clinical course in HRPBC.

Profiling of Sox4-dependent transcriptome in skin links tumour suppression and adult stem cell activation.

Adult stem cells (ASCs) reside in specific niches in a quiescent state in adult mammals. Upon specific cues they become activated and respond by self-renewing and differentiating into newly generated specialised cells that ensure appropriate tissue fitness. ASC quiescence also serves as a tumour suppression mechanism by hampering cellular transformation and expansion (White AC et al., 2014). Some genes restricted to early embryonic development and adult stem cell niches are often potent modulators of stem cell quiescence, and derailed expression of these is commonly associated to cancer (Vervoort SJ et al., 2013). Among them, it has been shown that recommissioned Sox4 expression facilitates proliferation, survival and migration of malignant cells. By generating a conditional Knockout mouse model in stratified epithelia (Sox4 (cKO) mice), we demonstrated a delayed plucking-induced Anagen in the absence of Sox4. Skin global transcriptome analysis revealed a prominent defect in the induction of transcriptional networks that control hair follicle stem cell (HFSC) activation such as those regulated by Wnt/Ctnnb1, Shh, Myc or Sox9, cell cycle and DNA damage response-associated pathways. Besides, Sox4 (cKO) mice are resistant to skin carcinogenesis, thus linking Sox4 to both normal and pathological HFSC activation (Foronda M et al., 2014). Here we provide additional details on the analysis of Sox4-regulated transcriptome in Telogen and Anagen skin. The raw and processed microarray data is deposited in GEO under GSE58155.

miRNA expression profiling of formalin-fixed paraffin-embedded (FFPE) hereditary breast tumors.

Hereditary breast cancer constitutes only 5-10% of all breast cancer cases and is characterized by strong family history of breast and/or other associated cancer types. Only ~ 25% of hereditary breast cancer cases carry a mutation in BRCA1 or BRCA2 gene, while mutations in other rare high and moderate-risk genes and common low penetrance variants may account for additional 20% of the cases. Thus the majority of cases are still unaccounted for and designated as BRCAX tumors. MicroRNAs are small non-coding RNAs that play important roles as regulators of gene expression and are deregulated in cancer. To characterize hereditary breast tumors based on their miRNA expression profiles we performed global microarray miRNA expression profiling on a retrospective cohort of 80 FFPE breast tissues, including 66 hereditary breast tumors (13 BRCA1, 10 BRCA2 and 43 BRCAX), 10 sporadic breast carcinomas and 4 normal breast tissues, using Exiqon miRCURY LNA™ microRNA Array v.11.0. Here we describe in detail the miRNA microarray expression data and tumor samples used for the study of BRCAX tumor heterogeneity (Tanic et al., 2013) and biomarkers associated with positive BRCA1/2 mutation status (Tanic et al., 2014). Additionally, we provide the R code for data preprocessing and quality control.

Non-coding recurrent mutations in chronic lymphocytic leukaemia.

Chronic lymphocytic leukaemia (CLL) is a frequent disease in which the genetic alterations determining the clinicobiological behaviour are not fully understood. Here we describe a comprehensive evaluation of the genomic landscape of 452 CLL cases and 54 patients with monoclonal B-lymphocytosis, a precursor disorder. We extend the number of CLL driver alterations, including changes in ZNF292, ZMYM3, ARID1A and PTPN11. We also identify novel recurrent mutations in non-coding regions, including the 3' region of NOTCH1, which cause aberrant splicing events, increase NOTCH1 activity and result in a more aggressive disease. In addition, mutations in an enhancer located on chromosome 9p13 result in reduced expression of the B-cell-specific transcription factor PAX5. The accumulative number of driver alterations (0 to ≥4) discriminated between patients with differences in clinical behaviour. This study provides an integrated portrait of the CLL genomic landscape, identifies new recurrent driver mutations of the disease, and suggests clinical interventions that may improve the management of this neoplasia.

The UBC-40 Urothelial Bladder Cancer cell line index: a genomic resource for functional studies.

BACKGROUND: Urothelial bladder cancer is a highly heterogeneous disease. Cancer cell lines are useful tools for its study. This is a comprehensive genomic characterization of 40 urothelial bladder carcinoma (UBC) cell lines including information on origin, mutation status of genes implicated in bladder cancer (FGFR3, PIK3CA, TP53, and RAS), copy number alterations assessed using high density SNP arrays, uniparental disomy (UPD) events, and gene expression. RESULTS: Based on gene mutation patterns and genomic changes we identify lines representative of the FGFR3-driven tumor pathway and of the TP53/RB tumor suppressor-driven pathway. High-density array copy number analysis identified significant focal gains (1q32, 5p13.1-12, 7q11, and 7q33) and losses (i.e. 6p22.1) in regions altered in tumors but not previously described as affected in bladder cell lines. We also identify new evidence for frequent regions of UPD, often coinciding with regions reported to be lost in tumors. Previously undescribed chromosome X losses found in UBC lines also point to potential tumor suppressor genes. Cell lines representative of the FGFR3-driven pathway showed a lower number of UPD events. CONCLUSIONS: Overall, there is a predominance of more aggressive tumor subtypes among the cell lines. We provide a cell line classification that establishes their relatedness to the major molecularly-defined bladder tumor subtypes. The compiled information should serve as a useful reference to the bladder cancer research community and should help to select cell lines appropriate for the functional analysis of bladder cancer genes, for example those being identified through massive parallel sequencing.

iMSRC: converting a standard automated microscope into an intelligent screening platform.

Microscopy in the context of biomedical research is demanding new tools to automatically detect and capture objects of interest. The few extant packages addressing this need, however, have enjoyed limited uptake due to complexity of use and installation. To overcome these drawbacks, we developed iMSRC, which combines ease of use and installation with high flexibility and enables applications such as rare event detection and high-resolution tissue sample screening, saving time and resources.

miRGate: a curated database of human, mouse and rat miRNA-mRNA targets.

MicroRNAs (miRNAs) are small non-coding elements involved in the post-transcriptional down-regulation of gene expression through base pairing with messenger RNAs (mRNAs). Through this mechanism, several miRNA-mRNA pairs have been described as critical in the regulation of multiple cellular processes, including early embryonic development and pathological conditions. Many of these pairs (such as miR-15 b/BCL2 in apoptosis or BART-6/BCL6 in diffuse large B-cell lymphomas) were experimentally discovered and/or computationally predicted. Available tools for target prediction are usually based on sequence matching, thermodynamics and conservation, among other approaches. Nevertheless, the main issue on miRNA-mRNA pair prediction is the little overlapping results among different prediction methods, or even with experimentally validated pairs lists, despite the fact that all rely on similar principles. To circumvent this problem, we have developed miRGate, a database containing novel computational predicted miRNA-mRNA pairs that are calculated using well-established algorithms. In addition, it includes an updated and complete dataset of sequences for both miRNA and mRNAs 3'-Untranslated region from human (including human viruses), mouse and rat, as well as experimentally validated data from four well-known databases. The underlying methodology of miRGate has been successfully applied to independent datasets providing predictions that were convincingly validated by functional assays. miRGate is an open resource available at http://mirgate.bioinfo.cnio.es. For programmatic access, we have provided a representational state transfer web service application programming interface that allows accessing the database at http://mirgate.bioinfo.cnio.es/API/ Database URL: http://mirgate.bioinfo.cnio.es

Nuclear DICKKOPF-1 as a biomarker of chemoresistance and poor clinical outcome in colorectal cancer.

Sporadic colorectal cancer (CRC) insurgence and progression depend on the activation of Wnt/ß-catenin signaling. Dickkopf (DKK)-1 is an extracellular inhibitor of Wnt/ß-catenin signaling that also has undefined ß-catenin-independent actions. Here we report for the first time that a proportion of DKK-1 locates within the nucleus of healthy small intestine and colon mucosa, and of CRC cells at specific chromatin sites of active transcription. Moreover, we show that DKK-1 regulates several cancer-related genes including the cancer stem cell marker aldehyde dehydrogenase 1A1 (ALDH1A1) and Ral-binding protein 1-associated Eps domain-containing 2 (REPS2), which are involved in detoxification of chemotherapeutic agents. Nuclear DKK-1 expression is lost along CRC progression; however, it remains high in a subset (15%) of CRC patients (n = 699) and associates with decreased progression-free survival (PFS) after chemotherapy administration and overall survival (OS) [adjusted HR, 1.65; 95% confidence interval (CI), 1.23-2.21; P = 0.002)]. Overexpression of ALDH1A1 and REPS2 associates with nuclear DKK-1 expression in tumors and correlates with decreased OS (P = 0.001 and 0.014) and PFS. In summary, our findings demonstrate a novel location of DKK-1 within the cell nucleus and support a role of nuclear DKK-1 as a predictive biomarker of chemoresistance in colorectal cancer.

The contribution of cohesin-SA1 to gene expression and chromatin architecture in two murine tissues.

Cohesin, which in somatic vertebrate cells consists of SMC1, SMC3, RAD21 and either SA1 or SA2, mediates higher-order chromatin organization. To determine how cohesin contributes to the establishment of tissue-specific transcriptional programs, we compared genome-wide cohesin distribution, gene expression and chromatin architecture in cerebral cortex and pancreas from adult mice. More than one third of cohesin binding sites differ between the two tissues and these show reduced overlap with CCCTC-binding factor (CTCF) and are enriched at the regulatory regions of tissue-specific genes. Cohesin/CTCF sites at active enhancers and promoters contain, at least, cohesin-SA1. Analyses of chromatin contacts at the Protocadherin (Pcdh) and Regenerating islet-derived (Reg) gene clusters, mostly expressed in brain and pancreas, respectively, revealed remarkable differences that correlate with the presence of cohesin. We could not detect significant changes in the chromatin contacts at the Pcdh locus when comparing brains from wild-type and SA1 null embryos. In contrast, reduced dosage of SA1 altered the architecture of the Reg locus and decreased the expression of Reg genes in the pancreas of SA1 heterozygous mice. Given the role of Reg proteins in inflammation, such reduction may contribute to the increased incidence of pancreatic cancer observed in these animals.

MicroRNA expression signatures for the prediction of BRCA1/2 mutation-associated hereditary breast cancer in paraffin-embedded formalin-fixed breast tumors.

Screening for germline mutations in breast cancer-associated genes BRCA1 and BRCA2 is indicated for patients with breast cancer from high-risk breast cancer families and influences both treatment options and clinical management. However, only 25% of selected patients test positive for BRCA1/2 mutation, indicating that additional diagnostic biomarkers are necessary. We analyzed 124 formalin-fixed paraffin-embedded (FFPE) tumor samples from patients with hereditary (104) and sporadic (20) invasive breast cancer, divided into two series (A and B). Microarray expression profiling of 829 human miRNAs was performed on 76 samples (Series A), and bioinformatics tool Prophet was used to develop and test a microarray classifier. Samples were stratified into a training set (n = 38) for microarray classifier generation and a test set (n = 38) for signature validation. A 35-miRNA microarray classifier was generated for the prediction of BRCA1/2 mutation status with a reported 95% (95% CI = 0.88-1.0) and 92% (95% CI: 0.84-1.0) accuracy in the training and the test set, respectively. Differential expression of 12 miRNAs between BRCA1/2 mutation carriers versus noncarriers was validated by qPCR in an independent tumor series B (n = 48). Logistic regression model based on the expression of six miRNAs (miR-142-3p, miR-505*, miR-1248, miR-181a-2*, miR-25* and miR-340*) discriminated between tumors from BRCA1/2 mutation carriers and noncarriers with 92% (95% CI: 0.84-0.99) accuracy. In conclusion, we identified miRNA expression signatures predictive of BRCA1/2 mutation status in routinely available FFPE breast tumor samples, which may be useful to complement current patient selection criteria for gene testing by identifying individuals with high likelihood of being BRCA1/2 mutation carriers.