Search

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.

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.

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 D 3 (1,25(OH) 2 D 3 ) on CRC stroma., Design: Expression of vitamin D receptor (VDR) and two 1,25(OH) 2 D 3 target genes was analysed in 658 patients with CRC with prolonged clinical follow-up. 1,25(OH) 2 D 3 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) 2 D 3 -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) 2 D 3 inhibited the protumoural activation of NFs and CAFs and imposed in CAFs a 1,25(OH) 2 D 3 -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) 2 D 3 has protective effects against CRC through the regulation of stromal fibroblasts. Accordingly, expression of VDR and 1,25(OH) 2 D 3 -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., (Published by the BMJ Publishing Group Limited. For permission to use (where not already granted under a licence) please go to http://www.bmj.com/company/products-services/rights-and-licensing/.)

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 .

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., (© 2016 The Authors. Aging Cell published by the Anatomical Society and John Wiley & Sons Ltd.)

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., Competing Interests: This work was supported by Eli Lilly and Company. B.H., H.B.B. and R.M.C. are employees of Eli Lilly and Company.

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.

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.

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.

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.

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.

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.

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, (© The Author(s) 2015. Published by Oxford University Press.)

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., (© The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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.

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., (© 2014 UICC.)

Molecular mechanisms responsible for hepatocellular carcinoma (HCC) remain largely unknown. Using genetically engineered mouse models, we show that hepatocyte-specific expression of unconventional prefoldin RPB5 interactor (URI) leads to a multistep process of HCC development, whereas its genetic reduction in hepatocytes protects against diethylnitrosamine (DEN)-induced HCC. URI inhibits aryl hydrocarbon (AhR)- and estrogen receptor (ER)-mediated transcription of enzymes implicated in L-tryptophan/kynurenine/nicotinamide adenine dinucleotide (NAD(+)) metabolism, thereby causing DNA damage at early stages of tumorigenesis. Restoring NAD(+) pools with nicotinamide riboside (NR) prevents DNA damage and tumor formation. Consistently, URI expression in human HCC is associated with poor survival and correlates negatively with L-tryptophan catabolism pathway. Our results suggest that boosting NAD(+) can be prophylactic or therapeutic in HCC., (Copyright © 2014 Elsevier Inc. All rights reserved.)

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.

Telomeric RNAs (TERRAs) are UUAGGG repeat-containing RNAs that are transcribed from the subtelomere towards the telomere. The precise genomic origin of TERRA has remained elusive. Using a whole-genome RNA-sequencing approach, we identify novel mouse transcripts arising mainly from the subtelomere of chromosome 18, and to a lesser extend chromosome 9, that resemble TERRA in several key aspects. Those transcripts contain UUAGGG-repeats and are heterogeneous in size, fluctuate in abundance in a TERRA-like manner during the cell cycle, are bound by TERRA RNA-binding proteins and are regulated in a manner similar to TERRA in response to stress and the induction of pluripotency. These transcripts are also found to associate with nearly all chromosome ends and downregulation of the transcripts that originate from chromosome 18 causes a reduction in TERRA abundance. Interestingly, downregulation of either chromosome 18 transcripts or TERRA results in increased number of telomere dysfunction-induced foci, suggesting a protective role at telomeres.

Sox4 expression is restricted in mammals to embryonic structures and some adult tissues, such as lymphoid organs, pancreas, intestine, and skin. During embryogenesis, Sox4 regulates mesenchymal and neural progenitor survival, as well as lymphocyte and myeloid differentiation, and contributes to pancreas, bone, and heart development. Aberrant Sox4 expression is linked to malignant transformation and metastasis in several types of cancer. To understand the role of Sox4 in the adult organism, we first generated mice with reduced whole-body Sox4 expression. These mice display accelerated aging and reduced cancer incidence. To specifically address a role for Sox4 in adult stem cells, we conditionally deleted Sox4 (Sox4(cKO)) in stratified epithelia. Sox4(cKO) mice show increased skin stem cell quiescence and resistance to chemical carcinogenesis concomitantly with downregulation of cell cycle, DNA repair, and activated hair follicle stem cell pathways. Altogether, these findings highlight the importance of Sox4 in regulating adult tissue homeostasis and cancer., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

Although common cancer hallmarks are well established, lineage-restricted oncogenes remain less understood. Here, we report an inherent dependency of melanoma cells on the small GTPase RAB7, identified within a lysosomal gene cluster that distinguishes this malignancy from over 35 tumor types. Analyses in human cells, clinical specimens, and mouse models demonstrated that RAB7 is an early-induced melanoma driver whose levels can be tuned to favor tumor invasion, ultimately defining metastatic risk. Importantly, RAB7 levels and function were independent of MITF, the best-characterized melanocyte lineage-specific transcription factor. Instead, we describe the neuroectodermal master modulator SOX10 and the oncogene MYC as RAB7 regulators. These results reveal a unique wiring of the lysosomal pathway that melanomas exploit to foster tumor progression., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a very rare disease that currently lacks genomic and genetic biomarkers to assist in its clinical management. We performed whole-exome sequencing (WES) of three BPDCN cases. Based on these data, we designed a resequencing approach to identify mutations in 38 selected genes in 25 BPDCN samples. WES revealed 37-99 deleterious gene mutations per exome with no common affected genes between patients, but with clear overlap in terms of molecular and disease pathways (hematological and dermatological disease). We identified for the first time deleterious mutations in IKZF3, HOXB9, UBE2G2 and ZEB2 in human leukemia. Target sequencing identified 29 recurring genes, ranging in prevalence from 36% for previously known genes, such as TET2, to 12-16% for newly identified genes, such as IKZF3 or ZEB2. Half of the tumors had mutations affecting either the DNA methylation or chromatin remodeling pathways. The clinical analysis revealed that patients with mutations in DNA methylation pathway had a significantly reduced overall survival (P=0.047). We provide the first mutational profiling of BPDCN. The data support the current WHO classification of the disease as a myeloid disorder and provide a biological rationale for the incorporation of epigenetic therapies for its treatment.

Cutaneous T-cell lymphoma (CTCL) is a heterogeneous group of primary cutaneous T-cell lymphoproliferative processes, mainly composed of mycosis fungoides and Sézary syndrome, the aggressive forms of which lack an effective treatment. The molecular pathogenesis of CTCL is largely unknown, although neoplastic cells show increased signaling from T-cell receptors (TCRs). DNAs from 11 patients with CTCL, both normal and tumoral, were target-enriched and sequenced by massive parallel sequencing for a selection of 524 TCR-signaling-related genes. Identified variants were validated by capillary sequencing. Multiple mutations were found that affected several signaling pathways, such as TCRs, nuclear factor κB, or Janus kinase/signal transducer and activator of transcription, but PLCG1 was found to be mutated in 3 samples, 2 of which featured a redundant mutation (c.1034T>C, S345F) in exon 11 that affects the PLCx protein catalytic domain. This mutation was further analyzed by quantitative polymerase chain reaction genotyping in a new cohort of 42 patients with CTCL, where it was found in 19% of samples. Immunohistochemical analysis for nuclear factor of activated T cells (NFAT) showed that PLCG1-mutated cases exhibited strong NFAT nuclear immunostaining. Functional studies demonstrated that PLCG1 mutants elicited increased downstream signaling toward NFAT activation, and inhibition of this pathway resulted in reduced CTCL cell proliferation and cell viability. Thus, increased proliferative and survival mechanisms in CTCL may partially depend on the acquisition of somatic mutations in PLCG1 and other genes that are essential for normal T-cell differentiation.

Chronic lymphocytic leukemia (CLL) has heterogeneous clinical and biological behavior. Whole-genome and -exome sequencing has contributed to the characterization of the mutational spectrum of the disease, but the underlying transcriptional profile is still poorly understood. We have performed deep RNA sequencing in different subpopulations of normal B-lymphocytes and CLL cells from a cohort of 98 patients, and characterized the CLL transcriptional landscape with unprecedented resolution. We detected thousands of transcriptional elements differentially expressed between the CLL and normal B cells, including protein-coding genes, noncoding RNAs, and pseudogenes. Transposable elements are globally derepressed in CLL cells. In addition, two thousand genes-most of which are not differentially expressed-exhibit CLL-specific splicing patterns. Genes involved in metabolic pathways showed higher expression in CLL, while genes related to spliceosome, proteasome, and ribosome were among the most down-regulated in CLL. Clustering of the CLL samples according to RNA-seq derived gene expression levels unveiled two robust molecular subgroups, C1 and C2. C1/C2 subgroups and the mutational status of the immunoglobulin heavy variable (IGHV) region were the only independent variables in predicting time to treatment in a multivariate analysis with main clinico-biological features. This subdivision was validated in an independent cohort of patients monitored through DNA microarrays. Further analysis shows that B-cell receptor (BCR) activation in the microenvironment of the lymph node may be at the origin of the C1/C2 differences.

Telomeres are nucleoprotein structures at the ends of eukaryotic chromosomes that protect them from degradation, end-to-end fusions, and fragility. In mammals, telomeres are composed of TTAGGG tandem repeats bound by a protein complex called shelterin, which has fundamental roles in the regulation of telomere protection and length. The telomeric repeat binding factor 1 (TERF1 or TRF1) is one of the components of shelterin and has been shown to be essential for telomere protection. Telomeric repeats can also be found throughout the genome, as Internal or Interstitial Telomeric Sequences (ITSs). Some of the components of shelterin have been described to bind to ITSs as well as other extra-telomeric regions, which in the case of RAP1 exert a key role in transcriptional regulation. Here, we set to address whether TRF1 can be found at extra-telomeric sites both under normal conditions and upon induction of telomere shortening. In particular, we performed a ChIP-sequencing technique to map TRF1 binding sites in MEFs wild-type and deficient for the telomerase RNA component (Terc(-/-)), with increasingly short telomeres. Our findings indicate that TRF1 is exclusively located at telomeres both under normal conditions, as well as under extreme telomere shortening. These results indicate that in mice not all members of shelterin have extra-telomeric roles as it was described for RAP1.

Urothelial bladder cancer (UBC) is heterogeneous at the clinical, pathological and genetic levels. Tumor invasiveness (T) and grade (G) are the main factors associated with outcome and determine patient management. A discovery exome sequencing screen (n = 17), followed by a prevalence screen (n = 60), identified new genes mutated in this tumor coding for proteins involved in chromatin modification (MLL2, ASXL2 and BPTF), cell division (STAG2, SMC1A and SMC1B) and DNA repair (ATM, ERCC2 and FANCA). STAG2, a subunit of cohesin, was significantly and commonly mutated or lost in UBC, mainly in tumors of low stage or grade, and its loss was associated with improved outcome. Loss of expression was often observed in chromosomally stable tumors, and STAG2 knockdown in bladder cancer cells did not increase aneuploidy. STAG2 reintroduction in non-expressing cells led to reduced colony formation. Our findings indicate that STAG2 is a new UBC tumor suppressor acting through mechanisms that are different from its role in preventing aneuploidy.

Motivation: RUbioSeq has been developed to facilitate the primary and secondary analysis of re-sequencing projects by providing an integrated software suite of parallelized pipelines to detect exome variants (single-nucleotide variants and copy number variations) and to perform bisulfite-seq analyses automatically. RUbioSeq's variant analysis results have been already validated and published., Availability: http://rubioseq.sourceforge.net/.

Replication stress (RS) is a source of DNA damage that has been linked to cancer and aging, which is suppressed by the ATR kinase. In mice, reduced ATR levels in a model of the ATR-Seckel syndrome lead to RS and accelerated aging. Similarly, ATR-Seckel embryonic fibroblasts (MEF) accumulate RS and undergo cellular senescence. We previously showed that senescence of ATR-Seckel MEF cannot be rescued by p53-deletion. Here, we show that the genetic ablation of the INK4a/Arf locus fully rescues senescence on ATR mutant MEF, but also that induced by other conditions that generate RS such as low doses of hydroxyurea or ATR inhibitors. In addition, we show that a persistent exposure to RS leads to increased levels of INK4a/Arf products, revealing that INK4a/ARF behaves as a bona fide RS checkpoint. Our data reveal an unknown role for INK4a/ARF in limiting the expansion of cells suffering from persistent replication stress, linking this well-known tumor suppressor to the maintenance of genomic integrity.

Urothelial bladder cancer (UBC) is heterogeneous at the clinical, pathological, genetic, and epigenetic levels. Exome sequencing has identified ARID1A as a novel tumor suppressor gene coding for a chromatin remodeling protein that is mutated in UBC. Here, we assess ARID1A alterations in two series of patients with UBC. In the first tumor series, we analyze exons 2-20 in 52 primary UBC and find that all mutant tumors belong to the aggressive UBC phenotype (high grade non-muscle invasive and muscle invasive tumors) (P = 0.05). In a second series (n = 84), we assess ARID1A expression using immunohistochemistry, a surrogate for mutation analysis, and find that loss of expression increases with higher stage/grade, it is inversely associated with FGFR3 overexpression (P = 0.03) but it is not correlated with p53 overexpression (P = 0.30). We also analyzed the expression of cytokeratins in the same set of tumor and find, using unsupervised clustering, that tumors with ARID1A loss of expression are generally KRT5/6-low. In this patient series, loss of ARID1A expression is also associated with worse prognosis, likely reflecting the higher prevalence of losses found in tumors of higher stage and grade. The independent findings in these two sets of patients strongly support the notion that ARID1A inactivation is a key player in bladder carcinogenesis occurring predominantly in FGFR3 wild type tumors.

Medullary thyroid carcinoma accounts for 2% to 5% of thyroid malignancies, of which 75% are sporadic and the remaining 25% are hereditary and related to multiple endocrine neoplasia type 2 syndrome. Despite a genotype-phenotype correlation with specific germline RET mutations, knowledge of pathways specifically associated with each mutation and with non-RET-mutated sporadic MTC remains lacking. Gene expression patterns have provided a tool for identifying molecular events related to specific tumor types and to different clinical features that could help identify novel therapeutic targets. Using transcriptional profiling of 49 frozen MTC specimens classified as RET mutation, we identified PROM1, LOXL2, GFRA1, and DKK4 as related to RET(M918T) and GAL as related to RET(634) mutation. An independent series of 19 frozen and 23 formalin-fixed, paraffin-embedded (FFPE) MTCs was used for validation by RT-qPCR. Two tissue microarrays containing 69 MTCs were available for IHC assays. According to pathway enrichment analysis and gene ontology biological processes, genes associated with the MTC(M918T) group were involved mainly in proliferative, cell adhesion, and general malignant metastatic effects and with Wnt, Notch, NFκB, JAK/Stat, and MAPK signaling pathways. Assays based on silencing of PROM1 by siRNAs performed in the MZ-CRC-1 cell line, harboring RET(M918T), caused an increase in apoptotic nuclei, suggesting that PROM1 is necessary for survival of these cells. This is the first report of PROM1 overexpression among primary tumors., (Copyright © 2013 American Society for Investigative Pathology. Published by Elsevier Inc. All rights reserved.)

ARF is a tumour suppressor activated by oncogenic stress, which stabilizes p53. Although p53 is a key component of the response to DNA damage, a similar function for ARF has not been ascribed. Here we show that primary mouse and human cells lacking the tumour suppressor BRCA2 accumulate DNA damage, which triggers checkpoint signalling and ARF activation. Furthermore, senescence induced by Brca2 deletion in primary mouse and human cells is reversed by the loss of ARF, a phenotype recapitulated in cells lacking RAD51. Surprisingly, ARF is not necessary for p53 accumulation per se but for altering the spectrum of genes activated by this transcription factor. Specifically, ARF enables p53 transcription of Dusp4 and Dusp7, which encode a pair of phosphatases known to inactivate the MAP kinases ERK1/2. Our results ascribe a previously unanticipated function to the ARF tumour suppressor in genome integrity, controlled by replicative stress and ATM/ATR-dependent checkpoint responses.

We have extensively characterized the DNA methylomes of 139 patients with chronic lymphocytic leukemia (CLL) with mutated or unmutated IGHV and of several mature B-cell subpopulations through the use of whole-genome bisulfite sequencing and high-density microarrays. The two molecular subtypes of CLL have differing DNA methylomes that seem to represent epigenetic imprints from distinct normal B-cell subpopulations. DNA hypomethylation in the gene body, targeting mostly enhancer sites, was the most frequent difference between naive and memory B cells and between the two molecular subtypes of CLL and normal B cells. Although DNA methylation and gene expression were poorly correlated, we identified gene-body CpG dinucleotides whose methylation was positively or negatively associated with expression. We have also recognized a DNA methylation signature that distinguishes new clinico-biological subtypes of CLL. We propose an epigenomic scenario in which differential methylation in the gene body may have functional and clinical implications in leukemogenesis.

Secondary diversification of antibodies through somatic hypermutation (SHM) and class switch recombination (CSR) is a critical component of the immune response. Activation-induced deaminase (AID) initiates both processes by deaminating cytosine residues in immunoglobulin genes. The resulting U:G mismatch can be processed by alternative pathways to give rise to a mutation (SHM) or a DNA double-strand break (CSR). Central to this processing is the activity of uracil-N-glycosylase (UNG), an enzyme normally involved in error-free base excision repair. We used next generation sequencing to analyze the contribution of UNG to the resolution of AID-induced lesions. Loss- and gain-of-function experiments showed that UNG activity can promote both error-prone and high fidelity repair of U:G lesions. Unexpectedly, the balance between these alternative outcomes was influenced by the sequence context of the deaminated cytosine, with individual hotspots exhibiting higher susceptibility to UNG-triggered error-free or error-prone resolution. These results reveal UNG as a new molecular layer that shapes the specificity of AID-induced mutations and may provide new insights into the role of AID in cancer development.

Human aging cannot be fully understood in terms of the constrained genetic setting. Epigenetic drift is an alternative means of explaining age-associated alterations. To address this issue, we performed whole-genome bisulfite sequencing (WGBS) of newborn and centenarian genomes. The centenarian DNA had a lower DNA methylation content and a reduced correlation in the methylation status of neighboring cytosine--phosphate--guanine (CpGs) throughout the genome in comparison with the more homogeneously methylated newborn DNA. The more hypomethylated CpGs observed in the centenarian DNA compared with the neonate covered all genomic compartments, such as promoters, exonic, intronic, and intergenic regions. For regulatory regions, the most hypomethylated sequences in the centenarian DNA were present mainly at CpG-poor promoters and in tissue-specific genes, whereas a greater level of DNA methylation was observed in CpG island promoters. We extended the study to a larger cohort of newborn and nonagenarian samples using a 450,000 CpG-site DNA methylation microarray that reinforced the observation of more hypomethylated DNA sequences in the advanced age group. WGBS and 450,000 analyses of middle-age individuals demonstrated DNA methylomes in the crossroad between the newborn and the nonagenarian/centenarian groups. Our study constitutes a unique DNA methylation analysis of the extreme points of human life at a single-nucleotide resolution level.

Besides its well-established role in sister chromatid cohesion, cohesin has recently emerged as major player in the organization of interphase chromatin. Such important function is related to its ability to entrap two DNA segments also in cis, thereby facilitating long-range DNA looping which is crucial for transcriptional regulation, organization of replication factories and V(D)J recombination. Vertebrate somatic cells have two different versions of cohesin, containing Smc1, Smc3, Rad21/Scc1 and either SA1 or SA2, but their functional specificity has been largely ignored. We recently generated a knockout mouse model for the gene encoding SA1, and found that this protein is essential to complete embryonic development. Cohesin-SA1 mediates cohesion at telomeres, which is required for their replication. Telomere defects in SA1- deficient cells provoke chromosome segregation errors resulting in aneuploidy despite robust centromere cohesion. This aneuploidy could explain why heterozygous animals have an earlier onset of tumorigenesis. In addition, the genome-wide distribution of cohesin changes dramatically in the absence of SA1, and the complex shows reduced accumulation at promoters and CTCF sites. As a consequence, gene expression is altered, leading to downregulation of biological processes related to a developmental disorder linked to cohesin function, the Cornelia de Lange Syndrome (CdLS). These results point out a prominent role of cohesin-SA1 in transcriptional regulation, with clear implications in the etiology of CdLS.

Vertebrates have two cohesin complexes that consist of Smc1, Smc3, Rad21/Scc1 and either SA1 or SA2, but their functional specificity is unclear. Mouse embryos lacking SA1 show developmental delay and die before birth. Comparison of the genome-wide distribution of cohesin in wild-type and SA1-null cells reveals that SA1 is largely responsible for cohesin accumulation at promoters and at sites bound by the insulator protein CTCF. As a consequence, ablation of SA1 alters transcription of genes involved in biological processes related to Cornelia de Lange syndrome (CdLS), a genetic disorder linked to dysfunction of cohesin. We show that the presence of cohesin-SA1 at the promoter of myc and of protocadherin genes positively regulates their expression, a task that cannot be assumed by cohesin-SA2. Lack of SA1 also alters cohesin-binding pattern along some gene clusters and leads to dysregulation of genes within. We hypothesize that impaired cohesin-SA1 function in gene expression underlies the molecular aetiology of CdLS.

Background: The transcription factor Pax8 is essential for the differentiation of thyroid cells. However, there are few data on genes transcriptionally regulated by Pax8 other than thyroid-related genes. To better understand the role of Pax8 in the biology of thyroid cells, we obtained transcriptional profiles of Pax8-silenced PCCl3 thyroid cells using whole genome expression arrays and integrated these signals with global cis-regulatory sequencing studies performed by ChIP-Seq analysis, Results: Exhaustive analysis of Pax8 immunoprecipitated peaks demonstrated preferential binding to intragenic regions and CpG-enriched islands, which suggests a role of Pax8 in transcriptional regulation of orphan CpG regions. In addition, ChIP-Seq allowed us to identify Pax8 partners, including proteins involved in tertiary DNA structure (CTCF) and chromatin remodeling (Sp1), and these direct transcriptional interactions were confirmed in vivo. Moreover, both factors modulate Pax8-dependent transcriptional activation of the sodium iodide symporter (Nis) gene promoter. We ultimately combined putative and novel Pax8 binding sites with actual target gene expression regulation to define Pax8-dependent genes. Functional classification suggests that Pax8-regulated genes may be directly involved in important processes of thyroid cell function such as cell proliferation and differentiation, apoptosis, cell polarity, motion and adhesion, and a plethora of DNA/protein-related processes., Conclusion: Our study provides novel insights into the role of Pax8 in thyroid biology, exerted through transcriptional regulation of important genes involved in critical thyrocyte processes. In addition, we found new transcriptional partners of Pax8, which functionally cooperate with Pax8 in the regulation of thyroid gene transcription. Besides, our data demonstrate preferential location of Pax8 in non-promoter CpG regions. These data point to an orphan CpG island-mediated mechanism that represents a novel role of Pax8 in the transcriptional output of the thyrocyte.

Accurate lymphoma diagnosis, prognosis and therapy still require additional markers. We explore the potential relevance of microRNA (miRNA) expression in a large series that included all major B-cell non-Hodgkin lymphoma (NHL) types. The data generated were also used to identify miRNAs differentially expressed in Burkitt lymphoma (BL) and diffuse large B-cell lymphoma (DLBCL) samples. A series of 147 NHL samples and 15 controls were hybridized on a human miRNA one-color platform containing probes for 470 human miRNAs. Each lymphoma type was compared against the entire set of NHLs. BL was also directly compared with DLBCL, and 43 preselected miRNAs were analyzed in a new series of routinely processed samples of 28 BLs and 43 DLBCLs using quantitative reverse transcription-polymerase chain reaction. A signature of 128 miRNAs enabled the characterization of lymphoma neoplasms, reflecting the lymphoma type, cell of origin and/or discrete oncogene alterations. Comparative analysis of BL and DLBCL yielded 19 differentially expressed miRNAs, which were confirmed in a second confirmation series of 71 paraffin-embedded samples. The set of differentially expressed miRNAs found here expands the range of potential diagnostic markers for lymphoma diagnosis, especially when differential diagnosis of BL and DLBCL is required.

Association studies are the choice approach in the discovery of the genomic basis of complex traits. To carry out such analysis, researchers frequently need to (1) select optimally informative sets of Single Nucleotide Polymorphisms (SNPs) in candidate regions and (2) annotate the results of associations found by means of genome-wide SNP arrays. These are complex tasks, since many criteria have to be considered, including the SNPs' functional properties, technological information and haplotype frequencies in given populations. SYSNPs implements algorithms that allow for efficient and simultaneous consideration of all the relevant criteria to obtain sets of SNPs that properly cover arbitrarily large lists of genes or genomic regions. Complementarily, SYSNPs allows for comprehensive functional annotation of SNPs linked to any given marker SNP. SYSNPs dramatically reduces the effort needed for SNP selection from days of searching various databases to a few minutes using a simple browser.

Chronic lymphocytic leukemia (CLL) is a heterogeneous disease without a well-defined genetic alteration responsible for the onset of the disease. Several lines of evidence coincide in identifying stimulatory and growth signals delivered by B-cell receptor (BCR), and co-receptors together with NFkB pathway, as being the driving force in B-cell survival in CLL. However, the molecular mechanism responsible for this activation has not been identified. Based on the hypothesis that BCR activation may depend on somatic mutations of the BCR and related pathways we have performed a complete mutational screening of 301 selected genes associated with BCR signaling and related pathways using massive parallel sequencing technology in 10 CLL cases. Four mutated genes in coding regions (KRAS, SMARCA2, NFKBIE and PRKD3) have been confirmed by capillary sequencing. In conclusion, this study identifies new genes mutated in CLL, all of them in cases with progressive disease, and demonstrates that next-generation sequencing technologies applied to selected genes or pathways of interest are powerful tools for identifying novel mutational changes.

Here we perform whole-exome sequencing of samples from 105 individuals with chronic lymphocytic leukemia (CLL), the most frequent leukemia in adults in Western countries. We found 1,246 somatic mutations potentially affecting gene function and identified 78 genes with predicted functional alterations in more than one tumor sample. Among these genes, SF3B1, encoding a subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP), is somatically mutated in 9.7% of affected individuals. Further analysis in 279 individuals with CLL showed that SF3B1 mutations were associated with faster disease progression and poor overall survival. This work provides the first comprehensive catalog of somatic mutations in CLL with relevant clinical correlates and defines a large set of new genes that may drive the development of this common form of leukemia. The results reinforce the idea that targeting several well-known genetic pathways, including mRNA splicing, could be useful in the treatment of CLL and other malignancies.

Transcription factors are common targets of epigenetic inactivation in human cancer. Promoter hypermethylation and subsequent silencing of transcription factors can lead to further deregulation of their targets. In this study, we explored the potential epigenetic deregulation in cancer of Ikaros family genes, which code for essential transcription factors in cell differentiation and exhibit genetic defects in hematologic neoplasias. Unexpectedly, our analysis revealed that Ikaros undergoes very specific promoter hypermethylation in colorectal cancer, including in all the cell lines studied and around 64% of primary colorectal adenocarcinomas, with increasing proportions in advanced Duke's stages. Ikaros hypermethylation occurred in the context of a novel long-range epigenetic silencing (LRES) region. Reintroduction of Ikaros in colorectal cancer cells, ChIP-chip analysis, and validation in primary samples led us to identify a number of direct targets that are possibly related with colorectal cancer progression. Our results not only provide the first evidence that LRES can have functional specific effects in cancer but also identify several deregulated Ikaros targets that may contribute to progression in colorectal adenocarcinoma.

Diffuse large B-cell lymphoma (DLBCL) prognostication requires additional biologic markers. miRNAs may constitute markers for cancer diagnosis, outcome, or therapy response. In the present study, we analyzed the miRNA expression profile in a retrospective multicenter series of 258 DLBCL patients uniformly treated with chemoimmunotherapy. Findings were correlated with overall survival (OS) and progression-free survival (PFS). miRNA and gene-expression profiles were studied using microarrays in an initial set of 36 cases. A selection of miRNAs associated with either DLBCL molecular subtypes (GCB/ABC) or clinical outcome were studied by multiplex RT-PCR in a test group of 240 cases with available formalin-fixed, paraffin-embedded (FFPE) diagnostic samples. The samples were divided into a training set (123 patients) and used to derive miRNA-based and combined (with IPI score) Cox regression models in an independent validation series (117 patients). Our model based on miRNA expression predicts OS and PFS and improves upon the predictions based on clinical variables. Combined models with IPI score identified a high-risk group of patients with a 2-year OS and a PFS probability of < 50%. In summary, a precise miRNA signature is associated with poor clinical outcome in chemoimmunotherapy-treated DLBCL patients. This information improves upon IPI-based predictions and identifies a subgroup of candidate patients for alternative therapeutic regimens.