The activities of pathways that regulate malignant transformation can be influenced by microRNAs (miRs). Recently, we showed that increased expression of five tumour-suppressor miRs, miR-508-3p, miR-508-5p, miR-509-3p, miR-509-5p and miR-130b-3p, correlate with improved clinical outcomes in human ovarian cancer patients, and that miR-509-3p attenuates invasion of ovarian cancer cell lines. Here, we investigate the mechanism underlying this reduced invasive potential by assessing the impact of these five miRs on the physical properties of cells. Human ovarian cancer cells (HEYA8, OVCAR8) that are transfected with miR mimics representing these five miRs exhibit decreased invasion through collagen matrices, increased cell size and reduced deformability as measured by microfiltration and microfluidic assays. To understand the molecular basis of altered invasion and deformability induced by these miRs, we use predicted and validated mRNA targets that encode structural and signalling proteins that regulate cell mechanical properties. Combined with analysis of gene transcripts by real-time PCR and image analysis of F-actin in single cells, our results suggest that these tumour-suppressor miRs may alter cell physical properties by regulating the actin cytoskeleton. Our findings provide biophysical insights into how tumour-suppressor miRs can regulate the invasive behaviour of ovarian cancer cells, and identify potential therapeutic targets that may be implicated in ovarian cancer progression.
cell deformability; actin cytoskeleton; microfluidics; microfiltration; tumour cell invasion
Chromosomal rearrangements are initiating events in acute lymphoblastic leukaemia (ALL). Here using RNA sequencing of 560 ALL cases, we identify rearrangements between MEF2D (myocyte enhancer factor 2D) and five genes (BCL9, CSF1R, DAZAP1, HNRNPUL1 and SS18) in 22 B progenitor ALL (B-ALL) cases with a distinct gene expression profile, the most common of which is MEF2D-BCL9. Examination of an extended cohort of 1,164 B-ALL cases identified 30 cases with MEF2D rearrangements, which include an additional fusion partner, FOXJ2; thus, MEF2D-rearranged cases comprise 5.3% of cases lacking recurring alterations. MEF2D-rearranged ALL is characterized by a distinct immunophenotype, DNA copy number alterations at the rearrangement sites, older diagnosis age and poor outcome. The rearrangements result in enhanced MEF2D transcriptional activity, lymphoid transformation, activation of HDAC9 expression and sensitive to histone deacetylase inhibitor treatment. Thus, MEF2D-rearranged ALL represents a distinct form of high-risk leukaemia, for which new therapeutic approaches should be considered.
Acute lymphoblastic leukaemia is characterized by chromosomal rearrangements. Here, the authors carry out RNA-sequencing on a large cohort of patients and identify recurrent rearrangements of MEF2D, which lead to increased transcriptional activity of the gene, and cellular transformation in vitro.
Aberrant activation or disruption of autophagy promotes tumorigenesis in various preclinical models of cancer, but whether the autophagy pathway is a target for recurrent molecular alteration in human cancer patient samples is unknown. To address this outstanding question, we surveyed 211 human autophagy-associated genes for tumor-related alterations to DNA sequence and RNA expression levels and examined their association with patient survival outcomes in multiple cancer types with sequence data from The Cancer Genome Atlas consortium. We found 3 (RB1CC1/FIP200, ULK4, WDR45/WIPI4) and one (ATG7) core autophagy genes to be under positive selection for somatic mutations in endometrial carcinoma and clear cell renal carcinoma, respectively, while 29 autophagy regulators and pathway interactors, including previously identified KEAP1, NFE2L2, and MTOR, were significantly mutated in 6 of the 11 cancer types examined. Gene expression analyses revealed that GABARAPL1 and MAP1LC3C/LC3C transcripts were less abundant in breast cancer and non-small cell lung cancers than in matched normal tissue controls; ATG4D transcripts were increased in lung squamous cell carcinoma, as were ATG16L2 transcripts in kidney cancer. Unsupervised clustering of autophagy-associated mRNA levels in tumors stratified patient overall survival in 3 of 9 cancer types (acute myeloid leukemia, clear cell renal carcinoma, and head and neck cancer). These analyses provide the first comprehensive resource of recurrently altered autophagy-associated genes in human tumors, and highlight cancer types and subtypes where perturbed autophagy may be relevant to patient overall survival.
AML; autophagy; cancer; gene expression; LAPTM4B; mutation; SNV; SQSTM1; TCGA; tumor
Papillary renal cell carcinoma, accounting for 15% of renal cell carcinoma, is a heterogeneous disease consisting of different types of renal cancer, including tumors with indolent, multifocal presentation and solitary tumors with an aggressive, highly lethal phenotype. Little is known about the genetic basis of sporadic papillary renal cell carcinoma; no effective forms of therapy for advanced disease exist.
We performed comprehensive molecular characterization utilizing whole-exome sequencing, copy number, mRNA, microRNA, methylation and proteomic analyses of 161 primary papillary renal cell carcinomas.
Type 1 and Type 2 papillary renal cell carcinomas were found to be different types of renal cancer characterized by specific genetic alterations, with Type 2 further classified into three individual subgroups based on molecular differences that influenced patient survival. MET alterations were associated with Type 1 tumors, whereas Type 2 tumors were characterized by CDKN2A silencing, SETD2 mutations, TFE3 fusions, and increased expression of the NRF2-ARE pathway. A CpG island methylator phenotype (CIMP) was found in a distinct subset of Type 2 papillary renal cell carcinoma characterized by poor survival and mutation of the fumarate hydratase (FH) gene.
Type 1 and Type 2 papillary renal cell carcinomas are clinically and biologically distinct. Alterations in the MET pathway are associated with Type 1 and activation of the NRF2-ARE pathway with Type 2; CDKN2A loss and CIMP in Type 2 convey a poor prognosis. Furthermore, Type 2 papillary renal cell carcinoma consists of at least 3 subtypes based upon molecular and phenotypic features.
The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon–driven, functional genomic mouse model of medulloblastoma with ‘humanized’ in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.
Human cancers, including breast cancers, are comprised of clones differing in mutation content. Clones evolve dynamically in space and time following principles of Darwinian evolution1,2, underpinning important emergent features such as drug resistance and metastasis3–7. Human breast cancer xenoengraftment is used as a means of capturing and studying tumour biology, and breast tumour xenografts are generally assumed to be reasonable models of the originating tumours8–10. However the consequences and reproducibility of engraftment and propagation on the genomic clonal architecture of tumours has not been systematically examined at single cell resolution. Here we show by both deep genome and single cell sequencing methods, the clonal dynamics of initial engraftment and subsequent serial propagation of primary and metastatic human breast cancers in immunodeficient mice. In all 15 cases examined, clonal selection on engraftment was observed in both primary and metastatic breast tumours, varying in degree from extreme selective engraftment of minor (<5% of starting population) clones to moderate, polyclonal engraftment. Furthermore, ongoing clonal dynamics during serial passaging is a feature of tumours experiencing modest initial selection. Through single cell sequencing, we show that major mutation clusters estimated from tumour population sequencing relate predictably to the most abundant clonal genotypes, even in clonally complex and rapidly evolving cases. Finally, we show that similar clonal expansion patterns can emerge in independent grafts of the same starting tumour population, indicating that genomic aberrations can be reproducible determinants of evolutionary trajectories. Our results show that measurement of genomically defined clonal population dynamics will be highly informative for functional studies utilizing patient-derived breast cancer xenoengraftment.
PMID: 25470049 CAMSID: cams5567
DDX3X encodes a DEAD-box family RNA helicase (DDX3) commonly mutated in medulloblastoma, a highly aggressive cerebellar tumor affecting both children and adults. Despite being implicated in several facets of RNA metabolism, the nature and scope of DDX3′s interactions with RNA remain unclear. Here, we show DDX3 collaborates extensively with the translation initiation machinery through direct binding to 5′UTRs of nearly all coding RNAs, specific sites on the 18S rRNA, and multiple components of the translation initiation complex. Impairment of translation initiation is also evident in primary medulloblastomas harboring mutations in DDX3X, further highlighting DDX3′s role in this process. Arsenite-induced stress shifts DDX3 binding from the 5′UTR into the coding region of mRNAs concomitant with a general reduction of translation, and both the shift of DDX3 on mRNA and decreased translation are blunted by expression of a catalytically-impaired, medulloblastoma-associated DDX3R534H variant. Furthermore, despite the global repression of translation induced by arsenite, translation is preserved on select genes involved in chromatin organization in DDX3R534H-expressing cells. Thus, DDX3 interacts extensively with RNA and ribosomal machinery to help remodel the translation landscape in response to stress, while cancer-related DDX3 variants adapt this response to selectively preserve translation.
medulloblastoma; DDX3X; DDX3; RNA helicase; CLIP-seq
Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.
While implicated in therapeutic resistance, malignant progenitor cell cycle kinetics have been difficult to quantify in real-time. We developed an efficient lentiviral bicistronic fluorescent, ubiquitination-based cell cycle indicator reporter (Fucci2BL) to image live single progenitors on a defined niche coupled with cell cycle gene expression analysis. We have identified key differences in cell cycle regulatory gene expression and transit times between normal and chronic myeloid leukemia progenitors that may inform cancer stem cell eradication strategies.
We report the most common single nucleotide substitution/deletion mutations in Favorable Histology Wilms Tumors (FHWT) to occur within SIX1/2 (7% of 534 tumors) and microRNA processing genes (miRNAPG) DGCR8 and DROSHA (15% of 534 tumors). Comprehensive analysis of 77 FHWTs indicates that tumors with SIX1/2 and/or miRNAPG mutations show a pre-induction metanephric mesenchyme gene expression pattern and are significantly associated with both perilobar nephrogenic rests and 11p15 imprinting aberrations. Significantly decreased expression of mature Let-7a and the miR-200 family (responsible for mesenchymal-to-epithelial transition) in miRNAPG-mutant tumors is associated with an undifferentiated blastemal histology. The combination of SIX and miRNAPG mutations in the same tumor is associated with evidence of RAS activation and a higher rate of relapse and death.
In an attempt to assess potential treatment options, whole-genome and transcriptome sequencing were performed on a patient with an unclassifiable small lymphoproliferative disorder. Variants from genome sequencing were prioritized using a combination of comparative variant distributions in a spectrum of lymphomas, and meta-analyses of gene expression profiling. In this patient, the molecular variants that we believe to be most relevant to the disease presentation most strongly resemble a diffuse large B-cell lymphoma (DLBCL), whereas the gene expression data are most consistent with a low-grade chronic lymphocytic leukemia (CLL). The variant of greatest interest was a predicted NOTCH2-truncating mutation, which has been recently reported in various lymphomas.
increase in B cell number; lymphocytosis
Medullary thyroid cancer (MTC) is a malignancy of the calcitonin-producing parafollicular cells of the thyroid gland. Surgery is the only curative treatment for this cancer. External beam radiation therapy is reserved for adjuvant treatment of MTC with aggressive features. Targeted therapeutics vandetanib and cabozantinib are approved for the treatment of aggressive and metastatic tumors that are not amenable to surgery. The use of these multikinase inhibitors are supported by the observed overactivation of the RET oncoprotein in a large subpopulation of MTCs. However, not all patients carry oncogenic alterations of this kinase. Hence, there is still a need for comprehensive molecular characterization of MTC utilizing whole-genome and transcriptome-sequencing methodologies with the aim of identifying targetable mutations. Here, we describe the genomic profiles of two medullary thyroid cancers and report the presence of a putative oncogenic BRAF fusion in one. Such alterations, previously observed in other malignancies and known targets of available drugs, can benefit patients who currently have no treatment options.
medullary thyroid carcinoma; neoplasm of the endocrine system
The MEF2 transcription factors have roles in muscle, cardiac, skeletal, vascular, neural, blood and immune system cell development through their effects on cell differentiation, proliferation, apoptosis, migration, shape and metabolism. Altered MEF2 activity plays a role in human diseases and has recently been implicated in the development of several cancer types. In particular, MEF2B, the most divergent and least studied protein of the MEF2 family, has a role unique from its paralogs in non-Hodgkin lymphomas. The use of genome-scale technologies has enabled comprehensive MEF2 target gene sets to be identified, contributing to our understanding of MEF2 proteins as nodes in complex regulatory networks. This review surveys the molecular interactions of MEF2 proteins and their effects on cellular and organismal phenotypes. We include a discussion of the emerging roles of MEF2 proteins as oncogenes and tumor suppressors of cancer. Throughout this article we highlight similarities and differences between the MEF2 family proteins, including a focus on functions of MEF2B.
MEF2; transcription factor; cancer; gene regulation; developmental biology
Wilms tumour is an embryonal tumour of childhood that closely resembles the developing kidney. Genomic changes responsible for the development of the majority of Wilms tumours remain largely unknown. Here we identify recurrent mutations within Wilms tumours that involve the highly conserved YEATS domain of MLLT1 (ENL), a gene known to be involved in transcriptional elongation during early development. The mutant MLLT1 protein shows altered binding to acetylated histone tails. Moreover, MLLT1-mutant tumours show an increase in MYC gene expression and HOX dysregulation. Patients with MLLT1-mutant tumours present at a younger age and have a high prevalence of precursor intralobar nephrogenic rests. These data support a model whereby activating MLLT1 mutations early in renal development result in the development of Wilms tumour.
Wilms tumour is a rare renal neoplasm that primarily affects children but the genomic changes responsible for its development are currently largely unknown. In this study, the authors identify somatic mutations of the MLLT1 gene that are potentially involved in the aetiology of a subset of Wilms tumours.
Anaplastic thyroid carcinoma is the most undifferentiated form of thyroid cancer and one of the deadliest of all adult solid malignancies. Here we report the first genomic and transcriptomic profile of anaplastic thyroid cancer including those of several unique cell lines and outline novel potential drivers of malignancy and targets of therapy.
We describe whole genomic and transcriptomic profiles of 1 primary anaplastic thyroid tumor and 3 authenticated cell lines. Those profiles augmented by the transcriptomes of 4 additional and unique cell lines were compared to 58 pairs of papillary thyroid carcinoma and matched normal tissue transcriptomes from The Cancer Genome Atlas study.
The most prevalent mutations were those of TP53 and BRAF; repeated alterations of the epigenetic machinery such as frame-shift deletions of HDAC10 and EP300, loss of SMARCA2 and fusions of MECP2, BCL11A and SS18 were observed. Sequence data displayed aneuploidy and large regions of copy loss and gain in all genomes. Common regions of gain were however evident encompassing chromosomes 5p and 20q. We found novel anaplastic gene fusions including MKRN1-BRAF, FGFR2-OGDH and SS18-SLC5A11, all expressed in-frame fusions involving a known proto-oncogene. Comparison of the anaplastic thyroid cancer expression datasets with the papillary thyroid cancer and normal thyroid tissue transcriptomes suggested several known drug targets such as FGFRs, VEGFRs, KIT and RET to have lower expression levels in anaplastic specimens compared with both papillary thyroid cancers and normal tissues, confirming the observed lack of response to therapies targeting these pathways. Further integrative data analysis identified the mTOR signaling pathway as a potential therapeutic target in this disease.
Anaplastic thyroid carcinoma possessed heterogeneous and unique profiles revealing the significance of detailed molecular profiling of individual tumors and the treatment of each as a unique entity; the cell line sequence data promises to facilitate the more accurate and intentional drug screening studies for anaplastic thyroid cancer.
Electronic supplementary material
The online version of this article (doi:10.1186/s12885-015-1955-9) contains supplementary material, which is available to authorized users.
Anaplastic thyroid carcinoma; cell line; whole genome and transcriptome sequencing; FGFR2-OGDH fusion; SS18-SLC5A11 fusion; MKRN1-BRAF fusion; epigenetic alterations; mTOR signaling pathway; therapy targets
Recurrence of the childhood brain cancer medulloblastoma is almost universally fatal. The current strategy for clinical trials is to test novel therapies at the time of recurrence, using targets discovered at diagnosis. To define the relationship between naïve and post-treatment medulloblastoma, we developed a functional genomic transposon-driven mouse-model. Deep sequencing of the commonly inserted genes in diagnostic and recurrent tumors revealed dramatic divergence post-therapy, with recurrences demonstrating a convergence on the TP53 pathway. Loss of function insertions in Tp53 are clonal in the recurrence, result in loss of expression of the Tp53 target gene P21, and arise through selection of minor subclones present before therapy. The majority of insertions identified at diagnosis are no longer present at recurrence, and therefore represent futile targets for therapy. To validate these findings, we studied 36 pairs of human primary/recurrent tumors by whole genome sequencing. When comparing SNVs, indels or structural aberrations, we found minimal genetic overlap between the dominant clones at diagnosis versus recurrence. Recurrent tumors have more genomic aberrations than their paired diagnostic sample, and commonly display alterations in the TP53 (SHH), or the TP53 and TGF-β pathways (Group 3/4 tumors). Clonal dynamics between primary and recurrent tumors was measured using ultra-deep sequencing (1,500X). Overall, recurrent tumors share little genetic identity with their diagnostic tumor, arise through clonal selection, target specific signaling pathways to avoid cell death, and assume clinical dominance after therapy. Our findings undermine the key assumption that diagnostic medulloblastoma samples provide rational targets for therapy in recurrent disease, and provide a simple reason for the failure of prior trials of targeted therapy. Repeated genetic convergence on specific signaling pathways at the time of recurrence suggests that it may be possible to develop ‘anticipatory therapy’ delivered up front to prevent the later emergence of minor resistant clones driving recurrence.
INTRODUCTION: Comprehensive, genome-wide profiling and next-generation sequencing based studies have dramatically improved our understanding of pediatric brain tumor biology in recent years. However, the vast majority of these studies rely on the assumption that single biopsies are representative of all areas within a tumor. Intratumor heterogeneity comprises a common phenomenon previously described in glioblastoma multiforme and other tumors. Highly disparate genetic profiles of spatially separated areas within the same tumor may preclude development of personalized, molecularly targeted therapies based on single tumor biopsies. MATERIAL AND METHODS: To assess the degree of intratumor heterogeneity, we conducted multiregion whole exome sequencing, high-resolution DNA copy number analysis (Cytoscan HD) and DNA methylation profiling (Infinium HumanMethylation450 BeadChip) on over 25 distinct pediatric and adult brain tumors with a median of six biopsies per tumor (range 4-9). Histological entities comprised ATRT (n = 2), DIPG (n = 2), ependymoma (n = 1), glioblastoma (n = 10), medulloblastoma (n = 10), and medulloepithelioma (n = 1). We elucidated the degree of intratumor heterogeneity and subgroup affiliation using integrated genomics and unsupervised hierarchical clustering algorithms. RESULTS: Epigenetic signatures were highly similar in multiregion biopsies from a single tumor. However, we identified up to 250,000 CpG dinucleotides that were differentially methylated when determining the intertumor heterogeneity of DNA methylation patterns even within disease subgroups. In addition, pediatric brain tumors displayed highly similar focal and broad DNA copy number alterations compared to their adult counterparts. Multiregion sequencing further reinforced the relatively higher degree of intratumor homogeneity in pediatric brain tumors. Lastly, we showed that subgroup affiliation was stable in all multiregion biopsies in medulloblastomas and glioblastomas. CONCLUSIONS: Our results reveal that single biopsies are representative of the tumor genomics landscape in pediatric brain tumors and that DNA methylation based subgrouping is geographically stable.
Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Here, we describe the genomic landscape of 496 PTCs. We observed a low frequency of somatic alterations (relative to other carcinomas) and extended the set of known PTC driver alterations to include EIF1AX, PPM1D and CHEK2 and diverse gene fusions. These discoveries reduced the fraction of PTC cases with unknown oncogenic driver from 25% to 3.5%. Combined analyses of genomic variants, gene expression, and methylation demonstrated that different driver groups lead to different pathologies with distinct signaling and differentiation characteristics. Similarly, we identified distinct molecular subgroups of BRAF-mutant tumors and multidimensional analyses highlighted a potential involvement of oncomiRs in less-differentiated subgroups. Our results propose a reclassification of thyroid cancers into molecular subtypes that better reflect their underlying signaling and differentiation properties, which has the potential to improve their pathological classification and better inform the management of the disease.
Cell fate acquisition is heavily influenced by direct interactions between master regulators and tissue-specific enhancers. However, it remains unclear how lineage-specifying transcription factors, which are often expressed in both progenitor and mature cell populations, influence cell differentiation. Using in vivo mouse liver development as a model, we identified thousands of enhancers that are bound by the master regulators HNF4A and FOXA2 in a differentiation-dependent manner, subject to chromatin remodeling, and associated with differentially expressed target genes. Enhancers exclusively occupied in the embryo were found to be responsive to developmentally regulated TEAD2 and coactivator YAP1. Our data suggest that Hippo signaling may affect hepatocyte differentiation by influencing HNF4A and FOXA2 interactions with temporal enhancers. In summary, transcription factor-enhancer interactions are not only tissue specific but also differentiation dependent, which is an important consideration for researchers studying cancer biology or mammalian development and/or using transformed cell lines.
Given the success of targeted agents in specific populations it is expected that some degree of molecular biomarker testing will become standard of care for many, if not all, cancers. To facilitate this, cancer centers worldwide are experimenting with targeted “panel” sequencing of selected mutations. Recent advances in genomic technology enable the generation of genome-scale data sets for individual patients. Recognizing the risk, inherent in panel sequencing, of failing to detect meaningful somatic alterations, we sought to establish processes to integrate data from whole-genome analysis (WGA) into routine cancer care. Between June 2012 and August 2014, 100 adult patients with incurable cancers consented to participate in the Personalized OncoGenomics (POG) study. Fresh tumor and blood samples were obtained and used for whole-genome and RNA sequencing. Computational approaches were used to identify candidate driver mutations, genes, and pathways. Diagnostic and drug information were then sought based on these candidate “drivers.” Reports were generated and discussed weekly in a multidisciplinary team setting. Other multidisciplinary working groups were assembled to establish guidelines on the interpretation, communication, and integration of individual genomic findings into patient care. Of 78 patients for whom WGA was possible, results were considered actionable in 55 cases. In 23 of these 55 cases, the patients received treatments motivated by WGA. Our experience indicates that a multidisciplinary team of clinicians and scientists can implement a paradigm in which WGA is integrated into the care of late stage cancer patients to inform systemic therapy decisions.
Next-generation sequencing (NGS) technologies and data have revolutionized cancer research and are increasingly being deployed to guide clinicians in treatment decision-making. NGS technologies have allowed us to take an “omics” approach to cancer in order to reveal genomic, transcriptomic, and epigenomic landscapes of individual malignancies. Integrative multi-platform analyses are increasingly used in large-scale projects that aim to fully characterize individual tumours as well as general cancer types and subtypes. In this review, we examine how NGS technologies in particular have contributed to “omics” approaches in cancer research, allowing for large-scale integrative analyses that consider hundreds of tumour samples. These types of studies have provided us with an unprecedented wealth of information, providing the background knowledge needed to make small-scale (including “N of 1”) studies informative and relevant. We also take a look at emerging opportunities provided by NGS and state-of-the-art third-generation sequencing technologies, particularly in the context of translational research. Cancer research and care are currently poised to experience significant progress catalyzed by accessible sequencing technologies that will benefit both clinical- and research-based efforts.
cancer; next-generation sequencing; translational research; precision oncology medicine
The reference human genome sequence set the stage for studies of genetic variation and its association with human disease, but a similar reference has lacked for epigenomic studies. To address this need, the NIH Roadmap Epigenomics Consortium generated the largest collection to-date of human epigenomes for primary cells and tissues. Here, we describe the integrative analysis of 111 reference human epigenomes generated as part of the program, profiled for histone modification patterns, DNA accessibility, DNA methylation, and RNA expression. We establish global maps of regulatory elements, define regulatory modules of coordinated activity, and their likely activators and repressors. We show that disease and trait-associated genetic variants are enriched in tissue-specific epigenomic marks, revealing biologically-relevant cell types for diverse human traits, and providing a resource for interpreting the molecular basis of human disease. Our results demonstrate the central role of epigenomic information for understanding gene regulation, cellular differentiation, and human disease.