3D chromatin connectivity underlies replication origin efficiency in mouse embryonic stem cells.

In mammalian cells, chromosomal replication starts at thousands of origins at which replisomes are assembled. Replicative stress triggers additional initiation events from 'dormant' origins whose genomic distribution and regulation are not well understood. In this study, we have analyzed origin activity in mouse embryonic stem cells in the absence or presence of mild replicative stress induced by aphidicolin, a DNA polymerase inhibitor, or by deregulation of origin licensing factor CDC6. In both cases, we observe that the majority of stress-responsive origins are also active in a small fraction of the cell population in a normal S phase, and stress increases their frequency of activation. In a search for the molecular determinants of origin efficiency, we compared the genetic and epigenetic features of origins displaying different levels of activation, and integrated their genomic positions in three-dimensional chromatin interaction networks derived from high-depth Hi-C and promoter-capture Hi-C data. We report that origin efficiency is directly proportional to the proximity to transcriptional start sites and to the number of contacts established between origin-containing chromatin fragments, supporting the organization of origins in higher-level DNA replication factories.

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.

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.

Chromatin regulation by Histone H4 acetylation at Lysine 16 during cell death and differentiation in the myeloid compartment.

Histone H4 acetylation at Lysine 16 (H4K16ac) is a key epigenetic mark involved in gene regulation, DNA repair and chromatin remodeling, and though it is known to be essential for embryonic development, its role during adult life is still poorly understood. Here we show that this lysine is massively hyperacetylated in peripheral neutrophils. Genome-wide mapping of H4K16ac in terminally differentiated blood cells, along with functional experiments, supported a role for this histone post-translational modification in the regulation of cell differentiation and apoptosis in the hematopoietic system. Furthermore, in neutrophils, H4K16ac was enriched at specific DNA repeats. These DNA regions presented an accessible chromatin conformation and were associated with the cleavage sites that generate the 50 kb DNA fragments during the first stages of programmed cell death. Our results thus suggest that H4K16ac plays a dual role in myeloid cells as it not only regulates differentiation and apoptosis, but it also exhibits a non-canonical structural role in poising chromatin for cleavage at an early stage of neutrophil cell death.

Bicycle: a bioinformatics pipeline to analyze bisulfite sequencing data.

Summary: High-throughput sequencing of bisulfite-converted DNA is a technique used to measure DNA methylation levels. Although a considerable number of computational pipelines have been developed to analyze such data, none of them tackles all the peculiarities of the analysis together, revealing limitations that can force the user to manually perform additional steps needed for a complete processing of the data. This article presents bicycle, an integrated, flexible analysis pipeline for bisulfite sequencing data. Bicycle analyzes whole genome bisulfite sequencing data, targeted bisulfite sequencing data and hydroxymethylation data. To show how bicycle overtakes other available pipelines, we compared them on a defined number of features that are summarized in a table. We also tested bicycle with both simulated and real datasets, to show its level of performance, and compared it to different state-of-the-art methylation analysis pipelines. Availability and implementation: Bicycle is publicly available under GNU LGPL v3.0 license at http://www.sing-group.org/bicycle. Users can also download a customized Ubuntu LiveCD including bicycle and other bisulfite sequencing data pipelines compared here. In addition, a docker image with bicycle and its dependencies, which allows a straightforward use of bicycle in any platform (e.g. Linux, OS X or Windows), is also available. Contact: ograna@cnio.es or dgpena@uvigo.es. Supplementary information: Supplementary data are available at Bioinformatics online.

A novel high-content analysis tool reveals Rab8-driven cytoskeletal reorganization through Rho GTPases, calpain and MT1-MMP.

Rab8 is a small Ras-related GTPase that regulates polarized membrane transport to the plasma membrane. Here, we developed a high-content analysis (HCA) tool to dissect Rab8-mediated actin and focal adhesion reorganization that revealed that Rab8 activation significantly induced Rac1 and Tiam1 to mediate cortical actin polymerization and RhoA-dependent stress fibre disassembly. Rab8 activation increased Rac1 activity, whereas its depletion activated RhoA, which led to reorganization of the actin cytoskeleton. Rab8 was also associated with focal adhesions, promoting their disassembly in a microtubule-dependent manner. This Rab8 effect involved calpain, MT1-MMP (also known as MMP14) and Rho GTPases. Moreover, we demonstrate the role of Rab8 in the cell migration process. Indeed, Rab8 is required for EGF-induced cell polarization and chemotaxis, as well as for the directional persistency of intrinsic cell motility. These data reveal that Rab8 drives cell motility by mechanisms both dependent and independent of Rho GTPases, thereby regulating the establishment of cell polarity, turnover of focal adhesions and actin cytoskeleton rearrangements, thus determining the directionality of cell migration.

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.

Transcriptome characterization by RNA sequencing identifies a major molecular and clinical subdivision in chronic lymphocytic leukemia.

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.

Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia.

Chronic lymphocytic leukaemia (CLL), the most frequent leukaemia in adults in Western countries, is a heterogeneous disease with variable clinical presentation and evolution. Two major molecular subtypes can be distinguished, characterized respectively by a high or low number of somatic hypermutations in the variable region of immunoglobulin genes. The molecular changes leading to the pathogenesis of the disease are still poorly understood. Here we performed whole-genome sequencing of four cases of CLL and identified 46 somatic mutations that potentially affect gene function. Further analysis of these mutations in 363 patients with CLL identified four genes that are recurrently mutated: notch 1 (NOTCH1), exportin 1 (XPO1), myeloid differentiation primary response gene 88 (MYD88) and kelch-like 6 (KLHL6). Mutations in MYD88 and KLHL6 are predominant in cases of CLL with mutated immunoglobulin genes, whereas NOTCH1 and XPO1 mutations are mainly detected in patients with unmutated immunoglobulins. The patterns of somatic mutation, supported by functional and clinical analyses, strongly indicate that the recurrent NOTCH1, MYD88 and XPO1 mutations are oncogenic changes that contribute to the clinical evolution of the disease. To our knowledge, this is the first comprehensive analysis of CLL combining whole-genome sequencing with clinical characteristics and clinical outcomes. It highlights the usefulness of this approach for the identification of clinically relevant mutations in 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.

A unique role of cohesin-SA1 in gene regulation and development.

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.

Genetic inactivation of Cdk7 leads to cell cycle arrest and induces premature aging due to adult stem cell exhaustion.

Cyclin-dependent kinase (Cdk)7, the catalytic subunit of the Cdk-activating kinase (CAK) complex has been implicated in the control of cell cycle progression and of RNA polymerase II (RNA pol II)-mediated transcription. Genetic inactivation of the Cdk7 locus revealed that whereas Cdk7 is completely dispensable for global transcription, is essential for the cell cycle via phosphorylation of Cdk1 and Cdk2. In vivo, Cdk7 is also indispensable for cell proliferation except during the initial stages of embryonic development. Interestingly, widespread elimination of Cdk7 in adult tissues with low proliferative indexes had no phenotypic consequences. However, ablation of conditional Cdk7 alleles in tissues with elevated cellular turnover led to the efficient repopulation of these tissues with Cdk7-expressing cells most likely derived from adult stem cells that may have escaped the inactivation of their targeted Cdk7 alleles. This process, a physiological attempt to maintain tissue homeostasis, led to the attrition of adult stem cell pools and to the appearance of age-related phenotypes, including telomere shortening and early death.

PLCG1 mutations in cutaneous T-cell lymphomas.

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.

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.

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.

Distinct DNA methylomes of newborns and centenarians.

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.

RUbioSeq: a suite of parallelized pipelines to automate exome variation and bisulfite-seq analyses.

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/.

Differential gene expression of medullary thyroid carcinoma reveals specific markers associated with genetic conditions.

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.

Combinatorial effects of microRNAs to suppress the Myc oncogenic pathway.

Many mammalian transcripts contain target sites for multiple miRNAs, although it is not clear to what extent miRNAs may coordinately regulate single genes. We have mapped the interactions between down-regulated miRNAs and overexpressed target protein-coding genes in murine and human lymphomas. Myc, one of the hallmark oncogenes in these lymphomas, stands out as the up-regulated gene with the highest number of genetic interactions with down-regulated miRNAs in mouse lymphomas. The regulation of Myc by several of these miRNAs is confirmed by cellular and reporter assays. The same approach identifies MYC and multiple Myc targets as a preferential target of down-regulated miRNAs in human Burkitt lymphoma, a pathology characterized by translocated MYC oncogenes. These results indicate that several miRNAs must be coordinately down-regulated to enhance critical oncogenes, such as Myc. Some of these Myc-targeting miRNAs are repressed by Myc, suggesting that these tumors are a consequence of the unbalanced activity of Myc versus miRNAs.

miRNA expression in diffuse large B-cell lymphoma treated with chemoimmunotherapy.

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.