Mutational topography reflects clinical neuroblastoma heterogeneity

Neuroblastoma is a pediatric solid tumor characterized by strong clinical heterogeneity. Although clinical risk-defining genomic alterations exist in neuroblastomas, the mutational processes involved in their generation remain largely unclear. By examining the topography and mutational signatures derived from all variant classes, we identified co-occurring mutational footprints, which we termed mutational scenarios. We demonstrate that clinical neuroblastoma heterogeneity is associated with differences in the mutational processes driving these scenarios, linking risk-defining pathognomonic variants to distinct molecular processes. Whereas high-risk MYCN-amplified neuroblastomas were characterized by signs of replication slippage and stress, homologous recombination-associated signatures defined high-risk non-MYCN-amplified patients. Non-high-risk neuroblastomas were marked by footprints of chromosome mis-segregation and TOP1 mutational activity. Furthermore, analysis of subclonal mutations uncovered differential activity of these processes through neuroblastoma evolution. Thus, clinical heterogeneity of neuroblastoma patients can be linked to differences in the mutational processes that are active in their tumors.

Tumor necrosis factor overcomes immune evasion in p53-mutant medulloblastoma

Many immunotherapies act by enhancing the ability of cytotoxic T cells to kill tumor cells. Killing depends on T cell recognition of antigens presented by class I major histocompatibility complex (MHC-I) proteins on tumor cells. In this study, we showed that medulloblastomas lacking the p53 tumor suppressor do not express surface MHC-I and are therefore resistant to immune rejection. Mechanistically, this is because p53 regulates expression of the peptide transporter Tap1 and the aminopeptidase Erap1, which are required for MHC-I trafficking to the cell surface. In vitro, tumor necrosis factor (TNF) or lymphotoxin-β receptor agonist can rescue expression of Erap1, Tap1 and MHC-I on p53-mutant tumor cells. In vivo, low doses of TNF prolong survival and synergize with immune checkpoint inhibitors to promote tumor rejection. These studies identified p53 as a key regulator of immune evasion and suggest that TNF could be used to enhance sensitivity of tumors to immunotherapy.

Allele-Specific Transcript Abundance: A Pilot Study in Healthy Centenarians

The genetic basis of healthy aging and longevity remains largely unexplained. One hypothesis as to why long-lived individuals do not appear to have a lower number of common-complex disease variants, is that despite carrying risk variants, they express disease-linked alleles at a lower level than the wild-type alleles. Allele-specific abundance (ASA) is the different transcript abundance of the two haplotypes of a diploid individual. We sequenced the transcriptomes of four healthy centenarians and four mid-life controls. CIBERSORT was used to estimate blood cell fractions: neutrophils were the most abundant source of RNA, followed by CD8+ T cells, resting NK cells, and monocytes. ASA variants were more common in noncoding than coding regions. Centenarians and controls had a comparable distribution of ASA variants by predicted effect, and we did not observe an overall bias in expression toward major or minor alleles. Immune pathways were most highly represented among the gene set that showed ASA. Although we found evidence of ASA in disease-associated genes and transcription factors, we did not observe any differences in the pattern of expression between centenarians and controls in this small pilot study.

Publisher Correction: Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Extrachromosomal circular DNA drives oncogenic genome remodeling in neuroblastoma

Extrachromosomal circularization of DNA is an important genomic feature in cancer. However, the structure, composition and genome-wide frequency of extrachromosomal circular DNA have not yet been profiled extensively. Here, we combine genomic and transcriptomic approaches to describe the landscape of extrachromosomal circular DNA in neuroblastoma, a tumor arising in childhood from primitive cells of the sympathetic nervous system. Our analysis identifies and characterizes a wide catalog of somatically acquired and undescribed extrachromosomal circular DNAs. Moreover, we find that extrachromosomal circular DNAs are an unanticipated major source of somatic rearrangements, contributing to oncogenic remodeling through chimeric circularization and reintegration of circular DNA into the linear genome. Cancer-causing lesions can emerge out of circle-derived rearrangements and are associated with adverse clinical outcome. It is highly probable that circle-derived rearrangements represent an ongoing mutagenic process. Thus, extrachromosomal circular DNAs represent a multihit mutagenic process, with important functional and clinical implications for the origins of genomic remodeling in cancer.

TNF superfamily cytokines overcome immune evasion in medulloblastoma

Many immunotherapies act by enhancing T cell killing of tumor cells. CD8+ cytotoxic T cells recognize antigens presented by class I major histocompatibility complex (MHC-I) proteins on tumor cells. Here we show that medulloblastomas lacking the p53 tumor suppressor do not express surface MHC-I and are therefore resistant to immune rejection. Mechanistically, this is because p53 regulates expression of the peptide transporter Tap1 and the aminopeptidase Erap1, which are required for MHC-I trafficking to the cell surface. Treatment with tumor necrosis factor (TNF) or lymphotoxin beta receptor agonist (LTβRag) rescues expression of Erap1, Tap1 and MHC-I on p53-mutant tumor cells. In vivo, TNF treatment prolongs survival and markedly augments the efficacy of the immune checkpoint inhibitor anti-PD-1. These studies identify p53 as a key regulator of immune evasion in vivo, and suggest that TNF could be used to enhance sensitivity of p53-mutant tumors to immunotherapy.

Homologous Recombination Deficiency and Platinum-Based Therapy Outcomes in Advanced Breast Cancer

Purpose: Recent studies have identified mutation signatures of homologous recombination deficiency (HRD) in over 20% of breast cancers, as well as pancreatic, ovarian, and gastric cancers. There is an urgent need to understand the clinical implications of HRD signatures. Whereas BRCA1/2 mutations confer sensitivity to platinum-based chemotherapies, it is not yet clear whether mutation signatures can independently predict platinum response.Experimental Design: In this observational study, we sequenced tumor whole genomes (100× depth) and matched normals (60×) of 93 advanced-stage breast cancers (33 platinum-treated). We computed a published metric called HRDetect, independently trained to predict BRCA1/2 status, and assessed its capacity to predict outcomes on platinum-based chemotherapies. Clinical endpoints were overall survival (OS), total duration on platinum-based therapy (TDT), and radiographic evidence of clinical improvement (CI).Results: HRDetect predicted BRCA1/2 status with an area under the curve (AUC) of 0.94 and optimal threshold of 0.7. Elevated HRDetect was also significantly associated with CI on platinum-based therapy (AUC = 0.89; P = 0.006) with the same optimal threshold, even after adjusting for BRCA1/2 mutation status and treatment timing. HRDetect scores over 0.7 were associated with a 3-month extended median TDT (P = 0.0003) and 1.3-year extended median OS (P = 0.04).Conclusions: Our findings not only independently validate HRDetect, but also provide the first evidence of its association with platinum response in advanced breast cancer. We demonstrate that HRD mutation signatures may offer clinically relevant information independently of BRCA1/2 mutation status and hope this work will guide the development of clinical trials. Clin Cancer Res; 23(24); 7521-30. ©2017 AACR.

Whole genome and whole transcriptome genomic profiling of a metastatic eccrine porocarcinoma

Eccrine porocarcinomas (EPs) are rare malignant tumours of the intraepidermic sweat gland duct and most often arise from benign eccrine poromas. Some recurrent somatic genomic events have been identified in these malignancies, but very little is known about the complexity of their molecular pathophysiology. We describe the whole genome and whole transcriptome genomic profiling of a metastatic EP in a 66-year-old male patient with a previous history of localized porocarcinoma of the scalp. Whole genome and whole transcriptome genomic profiling was performed on the metastatic EP. Whole genome sequencing was performed on blood-derived DNA in order to allow a comparison between germline and somatic events. We found somatic copy losses of several tumour suppressor genes including APC, PTEN and CDKN2A, CDKN2B and CDKN1A. We identified a somatic hemizygous CDKN2A pathogenic splice site variant. De novo transcriptome assembly revealed abnormal splicing of CDKN2A p14ARF and p16INK4a. Elevated expression of oncogenes EGFR and NOTCH1 was noted and no somatic mutations were found in these genes. Wnt pathway somatic alterations were also observed. In conclusion, our results suggest that the molecular pathophysiology of malignant EP features high complexity and subtle interactions of multiple key genes. Cell cycle dysregulation and CDKN2A loss of function was found to be a new potential driver in EP tumourigenesis. Moreover, the combination of somatic copy number variants and abnormal gene expression perhaps partly related to epigenetic mechanisms, all likely contribute to the development of this rare malignancy in our patient.

Abstract A184: Clinical application of whole genome and transcriptome sequencing in cancer care

Aim: To investigate feasibility and utility of whole genome and RNA sequencing in cancer care. The use of comprehensive personal genomic information to guide cancer treatment decisions has gained momentum in recent years. Here we present the Personalized Oncogenomics (POG) project launched at British Columbia Cancer Agency in 2012. This project uses paired tumor/normal whole genome and transcriptome sequencing to characterize a patient's tumor and inform treatment options within a clinically relevant time frame. In the past 5 years, around 600 patients with metastatic cancers (including 49 pediatric cases) have been sequenced and analyzed, representing over 50 cancer types. The breadth and depth of data in POG enable the detection of multiple types of alterations, from simple mutations, indels, and copy number changes, to more complex alterations such as gene fusion, mutation signatures, and homologous recombination deficiency score. Incorporation of gene expression data through transcriptome sequencing informs on the impact of observed genomic alterations, and provides information regarding specific diagnosis via expression comparison to other cancer subtypes. All genomic and transcriptomic variants are integrated to build a personalized tumor-specific molecular profile, identify actionable items supported by an in-house knowledgebase and publically available molecular oncology, and characterize each individual tumor by intensive pathway analysis and literature search. Such multidimensional data also impose challenges in interpretation and communication. We developed a pipeline to translate complex genomic data into clinically actionable and hypothetical recommendations for the treatment of individual patients. The pipeline produces two reports: a panel-like report that contains known SNVs and fusions with therapeutic relevance from a collection of more than 4000 events, and a second comprehensive and manually curated report to fully characterize the tumor using the whole genome and transcriptome datasets. Genomic data are presented and discussed at a multidisciplinary molecular tumor board consisting of medical oncologists, pathologists, bioinformaticians, geneticists, and biologists. This approach has enabled clinicians to make informed clinical decisions based on the genomic data integrated with other clinical features, as well as to form new treatment-related hypothesis. POG has shown that use of both whole genome and transcriptome sequencing allows identification of therapeutic targets in a significant proportion of patients. Almost 80% cases were found to have one or more actionable alterations (100% in pediatric cases), and almost one-third of these are defined only using RNA data. Based on molecular tumor board discussion, patients are directed to clinical trials, positioned to genomically informed standard-of-care options, or treated with off-label drugs. With demonstrated effectiveness, the integrative approach developed by POG not only provides molecular insight and treatment options into individual tumors, but also provides a rich resource of molecular data with matched clinical information that will aid our understanding of tumor biology and therapy response mechanisms to inform treatment strategies in the future.

Citation Format: Yaoqing Shen, Martin R. Jones, Erin Pleasance, Melika Bonakdar, Carolyn Ch'ng, Caralyn Reisle, Laura Williamson, Elisa Majounie, Greg Taylor, Simon Chan, Brandon Pierce, Wei Zhang, Amir Muhammadzadeh, Eric Y. Zhao, Dustin Bleile, Karen Mungall, Nina Thiessen, Eric Chuah, Tina Wong, Richard Corbett, Yussanne Ma, Richard A. Moore, Andrew J. Mungall, Yongjun Zhao, Stephen Yip, Anna F. Lee, Rod Rassekh, Rebecca Deyell, Howard Lim, Daniel Renouf, Robyn Roscoe, Steven J.M Jones, Janessa Laskin, Marco A. Marra. Clinical application of whole genome and transcriptome sequencing in cancer care [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2017 Oct 26-30; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Ther 2018;17(1 Suppl):Abstract nr A184.

The whole-genome landscape of medulloblastoma subtypes

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and ‘enhancer hijacking’ events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Genomic analysis of 491 medulloblastoma samples, including methylation profiling of 1,256 cases, effectively assigns candidate drivers to most tumours across all molecular subgroups.

Medulloblastomas are highly malignant brain tumours that develop during childhood. Paul Northcott and colleagues analysed the whole-genome sequences of 491 medulloblastomas in order to characterize the genomic landscape across tumours and identify new drivers and mutational signatures. Their integrative genomic analyses, including methylation profiling of 1,256 medulloblastomas, identifies subgroup-specific driver mutations and suggests additional tumour subtypes. The authors assign driver mutations to a high proportion of the less well characterized Group 3 and Group 4, which together contribute to more than 60% of all medulloblastomas.

Sporadic Early-Onset Diffuse Gastric Cancers Have High Frequency of Somatic CDH1 Alterations, but Low Frequency of Somatic RHOA Mutations Compared With Late-Onset Cancers

BACKGROUND & AIMS

Early-onset gastric cancer, which develops in patients younger than most gastric cancers, is usually detected at advanced stages, has diffuse histologic features, and occurs more frequently in women. We investigated somatic genomic alterations associated with the unique characteristics of sporadic diffuse gastric cancers (DGCs) from younger patients.

METHODS

We conducted whole exome and RNA sequence analyses of 80 resected DGC samples from patients 45 years old or younger in Korea. Patients with pathogenic germline mutations in CDH1, TP53, and ATM were excluded from the onset of this analysis, given our focus on somatic alterations. We used MutSig2CV to evaluate the significance of mutated genes. We recruited 29 additional early-onset Korean DGC samples and performed SNP6.0 array and targeted sequencing analyses of these 109 early-onset DGC samples (54.1% female, median age, 38 years). We compared the SNP6.0 array and targeted sequencing data of the 109 early-onset DGC samples with those from diffuse-type stomach tumor samples collected from 115 patients in Korea who were 46 years or older (late onset) at the time of diagnosis (controls; 29.6% female, median age, 67 years). We compared patient survival times among tumors from different subgroups and with different somatic mutations. We performed gene silencing of RHOA or CDH1 in DGC cells with small interfering RNAs for cell-based assays.

RESULTS

We identified somatic mutations in the following genes in a significant number of early-onset DGCs: the cadherin 1 gene (CDH1), TP53, ARID1A, KRAS, PIK3CA, ERBB3, TGFBR1, FBXW7, RHOA, and MAP2K1. None of 109 early-onset DGC cases had pathogenic germline CDH1 mutations. A higher proportion of early-onset DGCs had mutations in CDH1 (42.2%) or TGFBR1 (7.3%) compared with control DGCs (17.4% and 0.9%, respectively) (P < .001 and P = .014 for CDH1 and TGFBR1, respectively). In contrast, a smaller proportion of early-onset DGCs contained mutations in RHOA (9.2%) than control DGCs (19.1%) (P = .033). Late-onset DGCs in The Cancer Genome Atlas also contained less frequent mutations in CDH1 and TGFBR1 and more frequent RHOA mutations, compared with early-onset DGCs. Early-onset DGCs from women contained significantly more mutations in CDH1 or TGFBR1 than early-onset DGCs from men. CDH1 alterations, but not RHOA mutations, were associated with shorter survival times in patients with early-onset DGCs (hazard ratio, 3.4; 95% confidence interval, 1.5-7.7). RHOA activity was reduced by an R5W substitution-the RHOA mutation most frequently detected in early-onset DGCs. Silencing of CDH1, but not RHOA, increased migratory activity of DGC cells.

CONCLUSIONS

In an integrative genomic analysis, we found higher proportions of early-onset DGCs to contain somatic mutations in CDH1 or TGFBR1 compared with late-onset DGCs. However, a smaller proportion of early-onset DGCs contained somatic mutations in RHOA than late-onset DGCs. CDH1 alterations, but not RHOA mutations, were associated with shorter survival times of patients, which might account for the aggressive clinical course of early-onset gastric cancer. Female predominance in early-onset gastric cancer may be related to relatively high rates of somatic CDH1 and TGFBR1 mutations in this population.

Abstract 2473: Breast cancer whole genomes link homologous recombination deficiency (HRD) with therapeutic outcomes

Background: Homologous recombination deficiency (HRD) is common in cancer - germline BRCA1 & BRCA2 mutations account for 5-10% of breast cancers and confer 85% lifetime risk. HRD cancers exhibit genomic instability and sensitivity to platinum-based therapy and PARP inhibitors. While not all causes of HRD are known, recent sequencing efforts have revealed genome-wide somatic mutation signatures that characterize the HRD genomic instability phenotype, also known as “BRCA-ness”. This provides a promising new assay to predict sensitivity to platinum-based therapy. Here, we integrate two whole-genome sequencing metrics to assess their association with therapeutic outcomes in a breast cancer cohort.

Methods: Whole-genome sequencing of 47 breast cancer tumors (100x coverage) and matched normals (60x) was performed on an Illumina HiSeq. Alignment, assembly, SNV calling, and loss of heterozygosity (LOH) detection were performed with BWA, ABySS, Strelka, and APOLLOH respectively. SNV signatures were deciphered by non-negative matrix factorization with Monte Carlo resampling. An HRD score comprised of LOH, telomeric allelic imbalance (TAI), and large scale transition (LST) counts was computed. Clinical endpoints were obtained by retrospective review of treatment and imaging reports. Analysis is ongoing in an independent validation cohort of 62 sequenced cases.

Results: The HRD-linked SNV signature was significantly associated with radiographic clinical response (CR) to platinum-based therapy (p=0.015). Logistic regression demonstrated a 59% improved odds of CR to platinum-based therapy per 1000 somatic SNVs attributed to HRD (odds ratio 1.16-2.50). Tumors carried up to 10,246 such SNVs and all patients with CR were among the top quartile. The LOH-TAI-LST score was correlated with SNV signature (r=0.6, p=7×10-6) and associated with CR (p=0.025). Notably, elevated HRD signatures associated with CR were identified in tumors with wild-type BRCA1/BRCA2 or variants of unknown significance. Tumors with above median HRD signatures were associated with a 69-day longer time to treatment failure and an 18% daily decreased probability of treatment failure per 1000 HRD-attributed SNVs (hazard ratio 0.71-0.95, p = 0.007).

Discussion: We found that HRD mutation signatures are associated with clinical response and longer time to treatment failure with platinum-based therapy. While similar benefits were observed in patients with somatic bi-allelic loss of BRCA1/BRCA2, such cases are less common (8% of our cohort) compared to those with elevated HRD signature. Thus, mutation signature methods may identify patients who stand to benefit from platinum-based therapy missed by BRCA screening alone.

Citation Format: Eric Y. Zhao, Yaoqing Shen, Erin Pleasance, Katayoon Kasaian, Martin R. Jones, Carolyn Ch'ng, Caralyn Reisle, Peter Eirew, Karen Mungall, Nina Thiessen, Yussanne Ma, Alexandra Fok, Andrew J. Mungall, Yongjun Zhao, Richard Moore, Diego Villa, Tamara Shenkier, Caroline Lohrisch, Stephen Chia, Stephen Yip, Karen Gelmon, Howard Lim, Sophie Sun, Kasmintan A. Schrader, Sean Young, Aly Karsan, Robyn Roscoe, Janessa Laskin, Marco A. Marra, Steven J. Jones. Breast cancer whole genomes link homologous recombination deficiency (HRD) with therapeutic outcomes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2473. doi:10.1158/1538-7445.AM2017-2473

Genome-Wide Profiles of Extra-cranial Malignant Rhabdoid Tumors Reveal Heterogeneity and Dysregulated Developmental Pathways

Malignant rhabdoid tumors (MRTs) are rare lethal tumors of childhood that most commonly occur in the kidney and brain. MRTs are driven by SMARCB1 loss, but the molecular consequences of SMARCB1 loss in extra-cranial tumors have not been comprehensively described and genomic resources for analyses of extra-cranial MRT are limited. To provide such data, we used whole-genome sequencing, whole-genome bisulfite sequencing, whole transcriptome (RNA-seq) and microRNA sequencing (miRNA-seq), and histone modification profiling to characterize extra-cranial MRTs. Our analyses revealed gene expression and methylation subgroups and focused on dysregulated pathways, including those involved in neural crest development.

Comprehensive Molecular Characterization of Papillary Renal-Cell Carcinoma

BACKGROUND

Papillary renal-cell carcinoma, which accounts for 15 to 20% of renal-cell carcinomas, is a heterogeneous disease that consists of various 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, and no effective forms of therapy for advanced disease exist.

METHODS

We performed comprehensive molecular characterization of 161 primary papillary renal-cell carcinomas, using whole-exome sequencing, copy-number analysis, messenger RNA and microRNA sequencing, DNA-methylation analysis, and proteomic analysis.

RESULTS

Type 1 and type 2 papillary renal-cell carcinomas were shown to be different types of renal cancer characterized by specific genetic alterations, with type 2 further classified into three individual subgroups on the basis of molecular differences associated with patient survival. Type 1 tumors were associated with MET alterations, whereas type 2 tumors were characterized by CDKN2A silencing, SETD2 mutations, TFE3 fusions, and increased expression of the NRF2-antioxidant response element (ARE) pathway. A CpG island methylator phenotype (CIMP) was observed in a distinct subgroup of type 2 papillary renal-cell carcinomas that was characterized by poor survival and mutation of the gene encoding fumarate hydratase (FH).

CONCLUSIONS

Type 1 and type 2 papillary renal-cell carcinomas were shown to be clinically and biologically distinct. Alterations in the MET pathway were associated with type 1, and activation of the NRF2-ARE pathway was associated with type 2; CDKN2A loss and CIMP in type 2 conveyed a poor prognosis. Furthermore, type 2 papillary renal-cell carcinoma consisted of at least three subtypes based on molecular and phenotypic features. (Funded by the National Institutes of Health.).

Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas

BACKGROUND

Diffuse low-grade and intermediate-grade gliomas (which together make up the lower-grade gliomas, World Health Organization grades II and III) have highly variable clinical behavior that is not adequately predicted on the basis of histologic class. Some are indolent; others quickly progress to glioblastoma. The uncertainty is compounded by interobserver variability in histologic diagnosis. Mutations in IDH, TP53, and ATRX and codeletion of chromosome arms 1p and 19q (1p/19q codeletion) have been implicated as clinically relevant markers of lower-grade gliomas.

METHODS

We performed genomewide analyses of 293 lower-grade gliomas from adults, incorporating exome sequence, DNA copy number, DNA methylation, messenger RNA expression, microRNA expression, and targeted protein expression. These data were integrated and tested for correlation with clinical outcomes.

RESULTS

Unsupervised clustering of mutations and data from RNA, DNA-copy-number, and DNA-methylation platforms uncovered concordant classification of three robust, nonoverlapping, prognostically significant subtypes of lower-grade glioma that were captured more accurately by IDH, 1p/19q, and TP53 status than by histologic class. Patients who had lower-grade gliomas with an IDH mutation and 1p/19q codeletion had the most favorable clinical outcomes. Their gliomas harbored mutations in CIC, FUBP1, NOTCH1, and the TERT promoter. Nearly all lower-grade gliomas with IDH mutations and no 1p/19q codeletion had mutations in TP53 (94%) and ATRX inactivation (86%). The large majority of lower-grade gliomas without an IDH mutation had genomic aberrations and clinical behavior strikingly similar to those found in primary glioblastoma.

CONCLUSIONS

The integration of genomewide data from multiple platforms delineated three molecular classes of lower-grade gliomas that were more concordant with IDH, 1p/19q, and TP53 status than with histologic class. Lower-grade gliomas with an IDH mutation either had 1p/19q codeletion or carried a TP53 mutation. Most lower-grade gliomas without an IDH mutation were molecularly and clinically similar to glioblastoma. (Funded by the National Institutes of Health.).

Burden of Common Complex Disease Variants in the Exomes of Two Healthy Centenarian Brothers

BACKGROUND

It is not understood whether long-term good health is promoted by the absence of disease risk variants, the presence of protective variants, or both. We characterized the exomes of two exceptionally healthy centenarian brothers aged 106 and 109 years who had never been diagnosed with cancer, cardiovascular disease, diabetes, Alzheimer's disease, or major pulmonary disease.

OBJECTIVE

The aim of this study was to gain insight into whether exceptional health and longevity are a result of carrying fewer disease-associated variants than typical individuals.

METHODS

We compared the number of disease-associated alleles, and the proportion of alleles predicted to be functionally damaging, between the centenarian brothers and published population data. Mitochondrial sequence reads were extracted from the exome data in order to analyze mitochondrial variants.

RESULTS

The brothers carry a similar number of common disease-associated variants and predicted damaging variants compared to reference groups. They did not carry any high-penetrance clinically actionable variants. They carry mitochondrial haplogroup T, and one brother has a single heteroplasmic variant.

CONCLUSION

Although our small sample size does not allow for definitive conclusions, a healthy aging and longevity phenotype is not necessarily due to a decreased burden of common disease-associated variants. Instead, it may be rare 'positive' variants that play a role in this desirable phenotype.

Hippo signaling influences HNF4A and FOXA2 enhancer switching during hepatocyte differentiation

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.

JAGuaR: junction alignments to genome for RNA-seq reads

JAGuaR is an alignment protocol for RNA-seq reads that uses an extended reference to increase alignment sensitivity. It uses BWA to align reads to the genome and reference transcript models (including annotated exon-exon junctions) specifically allowing for the possibility of a single read spanning multiple exons. Reads aligned to the transcript models are then re-mapped on to genomic coordinates, transforming alignments that span multiple exons into large-gapped alignments on the genome. While JAGuaR does not detect novel junctions, we demonstrate how JAGuaR generates fast and accurate transcriptome alignments, which allows for both sensitive and specific SNV calling.

[Squamous cell carcinoma developing in oral lichen planus]

Oral lichen planus is a mucosal inflammatory disease whose pathogenesis is unclear. The chronic inflammation leads to development of a squamous cell carcinoma in 1-2% of the patients; we present an exemplary case.

Complete genomic landscape of a recurring sporadic parathyroid carcinoma

Parathyroid carcinoma is a rare endocrine malignancy with an estimated incidence of less than 1 per million population. Excessive secretion of parathyroid hormone, extremely high serum calcium level, and the deleterious effects of hypercalcaemia are the clinical manifestations of the disease. Up to 60% of patients develop multiple disease recurrences and although long-term survival is possible with palliative surgery, permanent remission is rarely achieved. Molecular drivers of sporadic parathyroid carcinoma have remained largely unknown. Previous studies, mostly based on familial cases of the disease, suggested potential roles for the tumour suppressor MEN1 and proto-oncogene RET in benign parathyroid tumourigenesis, while the tumour suppressor HRPT2 and proto-oncogene CCND1 may also act as drivers in parathyroid cancer. Here, we report the complete genomic analysis of a sporadic and recurring parathyroid carcinoma. Mutational landscapes of the primary and recurrent tumour specimens were analysed using high-throughput sequencing technologies. Such molecular profiling allowed for identification of somatic mutations never previously identified in this malignancy. These included single nucleotide point mutations in well-characterized cancer genes such as mTOR, MLL2, CDKN2C, and PIK3CA. Comparison of acquired mutations in patient-matched primary and recurrent tumours revealed loss of PIK3CA activating mutation during the evolution of the tumour from the primary to the recurrence. Structural variations leading to gene fusions and regions of copy loss and gain were identified at a single-base resolution. Loss of the short arm of chromosome 1, along with somatic missense and truncating mutations in CDKN2C and THRAP3, respectively, provides new evidence for the potential role of these genes as tumour suppressors in parathyroid cancer. The key somatic mutations identified in this study can serve as novel diagnostic markers as well as therapeutic targets.

Barnacle: detecting and characterizing tandem duplications and fusions in transcriptome assemblies

Background

Chimeric transcripts, including partial and internal tandem duplications (PTDs, ITDs) and gene fusions, are important in the detection, prognosis, and treatment of human cancers.

Results

We describe Barnacle, a production-grade analysis tool that detects such chimeras in de novo assemblies of RNA-seq data, and supports prioritizing them for review and validation by reporting the relative coverage of co-occurring chimeric and wild-type transcripts. We demonstrate applications in large-scale disease studies, by identifying PTDs in MLL, ITDs in FLT3, and reciprocal fusions between PML and RARA, in two deeply sequenced acute myeloid leukemia (AML) RNA-seq datasets.

Conclusions

Our analyses of real and simulated data sets show that, with appropriate filter settings, Barnacle makes highly specific predictions for three types of chimeric transcripts that are important in a range of cancers: PTDs, ITDs, and fusions. High specificity makes manual review and validation efficient, which is necessary in large-scale disease studies. Characterizing an extended range of chimera types will help generate insights into progression, treatment, and outcomes for complex diseases.

Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia

BACKGROUND

Many mutations that contribute to the pathogenesis of acute myeloid leukemia (AML) are undefined. The relationships between patterns of mutations and epigenetic phenotypes are not yet clear.

METHODS

We analyzed the genomes of 200 clinically annotated adult cases of de novo AML, using either whole-genome sequencing (50 cases) or whole-exome sequencing (150 cases), along with RNA and microRNA sequencing and DNA-methylation analysis.

RESULTS

AML genomes have fewer mutations than most other adult cancers, with an average of only 13 mutations found in genes. Of these, an average of 5 are in genes that are recurrently mutated in AML. A total of 23 genes were significantly mutated, and another 237 were mutated in two or more samples. Nearly all samples had at least 1 nonsynonymous mutation in one of nine categories of genes that are almost certainly relevant for pathogenesis, including transcription-factor fusions (18% of cases), the gene encoding nucleophosmin (NPM1) (27%), tumor-suppressor genes (16%), DNA-methylation-related genes (44%), signaling genes (59%), chromatin-modifying genes (30%), myeloid transcription-factor genes (22%), cohesin-complex genes (13%), and spliceosome-complex genes (14%). Patterns of cooperation and mutual exclusivity suggested strong biologic relationships among several of the genes and categories.

CONCLUSIONS

We identified at least one potential driver mutation in nearly all AML samples and found that a complex interplay of genetic events contributes to AML pathogenesis in individual patients. The databases from this study are widely available to serve as a foundation for further investigations of AML pathogenesis, classification, and risk stratification. (Funded by the National Institutes of Health.).

Identification and analysis of murine pancreatic islet enhancers

AIMS/HYPOTHESIS

The paucity of information on the epigenetic barriers that are blocking reprogramming protocols, and on what makes a beta cell unique, has hampered efforts to develop novel beta cell sources. Here, we aimed to identify enhancers in pancreatic islets, to understand their developmental ontologies, and to identify enhancers unique to islets to increase our understanding of islet-specific gene expression.

METHODS

We combined H3K4me1-based nucleosome predictions with pancreatic and duodenal homeobox 1 (PDX1), neurogenic differentiation 1 (NEUROD1), v-Maf musculoaponeurotic fibrosarcoma oncogene family, protein A (MAFA) and forkhead box A2 (FOXA2) occupancy data to identify enhancers in mouse islets.

RESULTS

We identified 22,223 putative enhancer loci in in vivo mouse islets. Our validation experiments suggest that nearly half of these loci are active in regulating islet gene expression, with the remaining regions probably poised for activity. We showed that these loci have at least nine developmental ontologies, and that islet enhancers predominately acquire H3K4me1 during differentiation. We next discriminated 1,799 enhancers unique to islets and showed that these islet-specific enhancers have reduced association with annotated genes, and identified a subset that are instead associated with novel islet-specific long non-coding RNAs (lncRNAs).

CONCLUSIONS/INTERPRETATIONS

Our results indicate that genes with islet-specific expression and function tend to have enhancers devoid of histone methylation marks or, less often, that are bivalent or repressed, in embryonic stem cells and liver. Further, we identify a subset of enhancers unique to islets that are associated with novel islet-specific genes and lncRNAs. We anticipate that these data will facilitate the development of novel sources of functional beta cell mass.

The genetic landscape of high-risk neuroblastoma

Neuroblastoma is a malignancy of the developing sympathetic nervous system that often presents with widespread metastatic disease, resulting in survival rates of less than 50%¹. To determine the spectrum of somatic mutation in high-risk neuroblastoma, we studied 240 cases using a combination of whole exome, genome and transcriptome sequencing as part of the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative. Here we report a low median exonic mutation frequency of 0.60 per megabase (0.48 non-silent), and remarkably few recurrently mutated genes in these tumors. Genes with significant somatic mutation frequencies included ALK (9.2% of cases), PTPN11 (2.9%), ATRX (2.5%, an additional 7.1% had focal deletions), MYCN (1.7%, a recurrent p.Pro44Leu alteration), and NRAS (0.83%). Rare, potentially pathogenic germline variants were significantly enriched in ALK, CHEK2, PINK1, and BARD1. The relative paucity of recurrent somatic mutations in neuroblastoma challenges current therapeutic strategies reliant upon frequently altered oncogenic drivers.

Bromodomain-containing protein 4 (BRD4) regulates RNA polymerase II serine 2 phosphorylation in human CD4+ T cells

Transcriptional elongation by RNA polymerase II (Pol II) is regulated by positive transcription elongation factor b (P-TEFb) in association with bromodomain-containing protein 4 (BRD4). We used genome-wide chromatin immunoprecipitation sequencing in primary human CD4+ T cells to reveal that BRD4 co-localizes with Ser-2-phosphorylated Pol II (Pol II Ser-2) at both enhancers and promoters of active genes. Disruption of bromodomain-histone acetylation interactions by JQ1, a small-molecule bromodomain inhibitor, resulted in decreased BRD4 binding, reduced Pol II Ser-2, and reduced expression of lineage-specific genes in primary human CD4+ T cells. A large number of JQ1-disrupted BRD4 binding regions exhibited diacetylated H4 (lysine 5 and -8) and H3K27 acetylation (H3K27ac), which correlated with the presence of histone acetyltransferases and deacetylases. Genes associated with BRD4/H3K27ac co-occupancy exhibited significantly higher activity than those associated with H3K27ac or BRD4 binding alone. Comparison of BRD4 binding in T cells and in human embryonic stem cells revealed that enhancer BRD4 binding sites were predominantly lineage-specific. Our findings suggest that BRD4-driven Pol II phosphorylation at serine 2 plays an important role in regulating lineage-specific gene transcription in human CD4+ T cells.

Concurrent CIC mutations, IDH mutations, and 1p/19q loss distinguish oligodendrogliomas from other cancers

Oligodendroglioma is characterized by unique clinical, pathological, and genetic features. Recurrent losses of chromosomes 1p and 19q are strongly associated with this brain cancer but knowledge of the identity and function of the genes affected by these alterations is limited. We performed exome sequencing on a discovery set of 16 oligodendrogliomas with 1p/19q co-deletion to identify new molecular features at base-pair resolution. As anticipated, there was a high rate of IDH mutations: all cases had mutations in either IDH1 (14/16) or IDH2 (2/16). In addition, we discovered somatic mutations and insertions/deletions in the CIC gene on chromosome 19q13.2 in 13/16 tumours. These discovery set mutations were validated by deep sequencing of 13 additional tumours, which revealed seven others with CIC mutations, thus bringing the overall mutation rate in oligodendrogliomas in this study to 20/29 (69%). In contrast, deep sequencing of astrocytomas and oligoastrocytomas without 1p/19q loss revealed that CIC alterations were otherwise rare (1/60; 2%). Of the 21 non-synonymous somatic mutations in 20 CIC-mutant oligodendrogliomas, nine were in exon 5 within an annotated DNA-interacting domain and three were in exon 20 within an annotated protein-interacting domain. The remaining nine were found in other exons and frequently included truncations. CIC mutations were highly associated with oligodendroglioma histology, 1p/19q co-deletion, and IDH1/2 mutation (p < 0.001). Although we observed no differences in the clinical outcomes of CIC mutant versus wild-type tumours, in a background of 1p/19q co-deletion, hemizygous CIC mutations are likely important. We hypothesize that the mutant CIC on the single retained 19q allele is linked to the pathogenesis of oligodendrogliomas with IDH mutation. Our detailed study of genetic aberrations in oligodendroglioma suggests a functional interaction between CIC mutation, IDH1/2 mutation, and 1p/19q co-deletion.

Cell of origin in AML: susceptibility to MN1-induced transformation is regulated by the MEIS1/AbdB-like HOX protein complex

Pathways defining susceptibility of normal cells to oncogenic transformation may be valuable therapeutic targets. We characterized the cell of origin and its critical pathways in MN1-induced leukemias. Common myeloid (CMP) but not granulocyte-macrophage progenitors (GMP) could be transformed by MN1. Complementation studies of CMP-signature genes in GMPs demonstrated that MN1-leukemogenicity required the MEIS1/AbdB-like HOX-protein complex. ChIP-sequencing identified common target genes of MN1 and MEIS1 and demonstrated identical binding sites for a large proportion of their chromatin targets. Transcriptional repression of MEIS1 targets in established MN1 leukemias demonstrated antileukemic activity. As MN1 relies on but cannot activate expression of MEIS1/AbdB-like HOX proteins, transcriptional activity of these genes determines cellular susceptibility to MN1-induced transformation and may represent a promising therapeutic target.

SNP discovery in black cottonwood (Populus trichocarpa) by population transcriptome resequencing

The western black cottonwood (Populus trichocarpa) was the first tree to have its genome fully sequenced and has emerged as the model species for the study of secondary growth and wood formation. It is also a good candidate species for the production of lignocellulosic biofuels. Here, we present and make available to the research community the results of the sequencing of the transcriptome of developing xylem in 20 accessions with high-throughput next generation sequencing technology. We found over 0.5 million putative single nucleotide polymorphisms (SNPs) in 26,595 genes that are expressed in developing secondary xylem. More than two-thirds of all SNPs were found in annotated exons, with 18% and 14% in regions of the genome annotated as introns and intergenic, respectively, where only 3% and 4% of sequence reads mapped. This suggests that the current annotation of the poplar genome is remarkably incomplete and that there are many transcripts and novel genes waiting to be annotated. We hope that this resource will stimulate further research in expression profiling, detection of alternative splicing and adaptive evolution in poplar.

Class proximity measures--dissimilarity-based classification and display of high-dimensional data

For two-class problems, we introduce and construct mappings of high-dimensional instances into dissimilarity (distance)-based Class-Proximity Planes. The Class Proximity Projections are extensions of our earlier relative distance plane mapping, and thus provide a more general and unified approach to the simultaneous classification and visualization of many-feature datasets. The mappings display all L-dimensional instances in two-dimensional coordinate systems, whose two axes represent the two distances of the instances to various pre-defined proximity measures of the two classes. The Class Proximity mappings provide a variety of different perspectives of the dataset to be classified and visualized. We report and compare the classification and visualization results obtained with various Class Proximity Projections and their combinations on four datasets from the UCI data base, as well as on a particular high-dimensional biomedical dataset.

SDHA loss-of-function mutations in KIT-PDGFRA wild-type gastrointestinal stromal tumors identified by massively parallel sequencing

Approximately 10%-15% of gastrointestinal stromal tumors (GISTs) in adults do not harbor any mutation in the KIT or PDGFRA genes (ie, KIT/PDGFRA wild-type GISTs). Recently, mutations in SDHB and SDHC (which encode succinate dehydrogenase subunits B and C, respectively) but not in SDHA and SDHD (which encode subunits A and D, respectively) were identified in KIT/PDGFRA wild-type GISTs. To search for novel pathogenic mutations, we sequenced the tumor transcriptome of two young adult patients who developed sporadic KIT/PDGFRA wild-type GISTs by using a massively parallel sequencing approach. The only variants identified as disease related by computational analysis were in SDHA. One patient carried the homozygous nonsense mutation p.Ser384X, the other patient was a compound heterozygote harboring a p.Arg31X nonsense mutation and a p.Arg589Trp missense mutation. The heterozygous nonsense mutations in both patients were present in germline DNA isolated from peripheral blood. Protein structure analysis indicates that all three mutations lead to functional inactivation of the protein. This is the first report, to our knowle dge, that identifies SDHA inactivation as a common oncogenic event in GISTs that lack a mutation in KIT and PDGFRA.

De novo assembly and analysis of RNA-seq data

We describe Trans-ABySS, a de novo short-read transcriptome assembly and analysis pipeline that addresses variation in local read densities by assembling read substrings with varying stringencies and then merging the resulting contigs before analysis. Analyzing 7.4 gigabases of 50-base-pair paired-end Illumina reads from an adult mouse liver poly(A) RNA library, we identified known, new and alternative structures in expressed transcripts, and achieved high sensitivity and specificity relative to reference-based assembly methods.

Conserved role of intragenic DNA methylation in regulating alternative promoters

While the methylation of DNA in 5′ promoters suppresses gene expression, the role of DNA methylation in gene bodies is unclear^1–5. In mammals, tissue- and cell type-specific methylation is present in a small percentage of 5′ CpG island (CGI) promoters, while a far greater proportion occurs across gene bodies, coinciding with highly conserved sequences^5–10. Tissue-specific intragenic methylation might reduce,³ or, paradoxically, enhance transcription elongation efficiency^1,2,4,5. Capped analysis of gene expression (CAGE) experiments also indicate that transcription commonly initiates within and between genes^11–15. To investigate the role of intragenic methylation, we generated a map of DNA methylation from human brain encompassing 24.7 million of the 28 million CpG sites. From the dense, high-resolution coverage of CpG islands, the majority of methylated CpG islands were revealed to be in intragenic and intergenic regions, while less than 3% of CpG islands in 5′ promoters were methylated. The CpG islands in all three locations overlapped with RNA markers of transcription initiation, and unmethylated CpG islands also overlapped significantly with trimethylation of H3K4, a histone modification enriched at promoters¹⁶. The general and CpG-island-specific patterns of methylation are conserved in mouse tissues. An in-depth investigation of the human SHANK3 locus^17,18 and its mouse homologue demonstrated that this tissue-specific DNA methylation regulates intragenic promoter activity in vitro and in vivo. These methylation-regulated, alternative transcripts are expressed in a tissue and cell type-specific manner, and are expressed differentially within a single cell type from distinct brain regions. These results support a major role for intragenic methylation in regulating cell context-specific alternative promoters in gene bodies.

Evolution of an adenocarcinoma in response to selection by targeted kinase inhibitors

Background

Adenocarcinomas of the tongue are rare and represent the minority (20 to 25%) of salivary gland tumors affecting the tongue. We investigated the utility of massively parallel sequencing to characterize an adenocarcinoma of the tongue, before and after treatment.

Results

In the pre-treatment tumor we identified 7,629 genes within regions of copy number gain. There were 1,078 genes that exhibited increased expression relative to the blood and unrelated tumors and four genes contained somatic protein-coding mutations. Our analysis suggested the tumor cells were driven by the RET oncogene. Genes whose protein products are targeted by the RET inhibitors sunitinib and sorafenib correlated with being amplified and or highly expressed. Consistent with our observations, administration of sunitinib was associated with stable disease lasting 4 months, after which the lung lesions began to grow. Administration of sorafenib and sulindac provided disease stabilization for an additional 3 months after which the cancer progressed and new lesions appeared. A recurring metastasis possessed 7,288 genes within copy number amplicons, 385 genes exhibiting increased expression relative to other tumors and 9 new somatic protein coding mutations. The observed mutations and amplifications were consistent with therapeutic resistance arising through activation of the MAPK and AKT pathways.

Conclusions

We conclude that complete genomic characterization of a rare tumor has the potential to aid in clinical decision making and identifying therapeutic approaches where no established treatment protocols exist. These results also provide direct in vivo genomic evidence for mutational evolution within a tumor under drug selection and potential mechanisms of drug resistance accrual.

System-level analysis of neuroblastoma tumor-initiating cells implicates AURKB as a novel drug target for neuroblastoma

PURPOSE

Neuroblastoma (NB) is an aggressive tumor of the developing peripheral nervous system that remains difficult to cure in the advanced stages. The poor prognosis for high-risk NB patients is associated with common disease recurrences that fail to respond to available therapies. NB tumor-initiating cells (TICs), isolated from metastases and primary tumors, may escape treatment and contribute to tumor relapse. New therapies that target the TICs may therefore prevent or treat tumor recurrences.

EXPERIMENTAL DESIGN

We undertook a system-level characterization of NB TICs to identify potential drug targets against recurrent NB. We used next-generation RNA sequencing and/or human exon arrays to profile the transcriptomes of 11 NB TIC lines from six NB patients, revealing genes that are highly expressed in the TICs compared with normal neural crest-like cells and unrelated cancer tissues. We used gel-free two-dimensional liquid chromatography coupled to shotgun tandem mass spectrometry to confirm the presence of proteins corresponding to the most abundant TIC-enriched transcripts, thereby providing validation to the gene expression result.

RESULTS

Our study revealed that genes in the BRCA1 signaling pathway are frequently misexpressed in NB TICs and implicated Aurora B kinase as a potential drug target for NB therapy. Treatment with a selective AURKB inhibitor was cytotoxic to NB TICs but not to the normal neural crest-like cells.

CONCLUSION

This work provides the first high-resolution system-level analysis of the transcriptomes of 11 primary human NB TICs and identifies a set of candidate NB TIC-enriched transcripts for further development as therapeutic targets.

Locus co-occupancy, nucleosome positioning, and H3K4me1 regulate the functionality of FOXA2-, HNF4A-, and PDX1-bound loci in islets and liver

The liver and pancreas share a common origin and coexpress several transcription factors. To gain insight into the transcriptional networks regulating the function of these tissues, we globally identify binding sites for FOXA2 in adult mouse islets and liver, PDX1 in islets, and HNF4A in liver. Because most eukaryotic transcription factors bind thousands of loci, many of which are thought to be inactive, methods that can discriminate functionally active binding events are essential for the interpretation of genome-wide transcription factor binding data. To develop such a method, we also generated genome-wide H3K4me1 and H3K4me3 localization data in these tissues. By analyzing our binding and histone methylation data in combination with comprehensive gene expression data, we show that H3K4me1 enrichment profiles discriminate transcription factor occupied loci into three classes: those that are functionally active, those that are poised for activation, and those that reflect pioneer-like transcription factor activity. Furthermore, we demonstrate that the regulated presence of H3K4me1-marked nucleosomes at transcription factor occupied promoters and enhancers controls their activity, implicating both tissue-specific transcription factor binding and nucleosome remodeling complex recruitment in determining tissue-specific gene expression. Finally, we apply these approaches to generate novel insights into how FOXA2, PDX1, and HNF4A cooperate to drive islet- and liver-specific gene expression.

Caenorhabditis elegans cisRED: a catalogue of conserved genomic elements

The availability of completely sequenced genomes from eight species of nematodes has provided an opportunity to identify novel cis-regulatory elements in the promoter regions of Caenorhabditis elegans transcripts using comparative genomics. We determined orthologues for C. elegans transcripts in C. briggsae, C. remanei, C. brenneri, C. japonica, Pristionchus pacificus, Brugia malayi and Trichinella spiralis using the WABA alignment algorithm. We pooled the upstream region of each transcript in C. elegans with the upstream regions of its orthologues and identified conserved DNA sequence elements by de novo motif discovery. In total, we discovered 158 017 novel conserved motifs upstream of 3847 C. elegans transcripts for which three or more orthologues were available, and identified 82% of 44 experimentally proven regulatory elements from ORegAnno. We annotated 26% of the motifs as similar to known binding sequences of transcription factors from ORegAnno, TRANSFAC and JASPAR. This is the first catalogue of annotated conserved upstream elements for nematodes and can be used to find putative regulatory elements, improve gene models, discover novel RNA genes, and understand the evolution of transcription factors and their binding sites in phylum Nematoda. The annotated motifs provide novel binding site candidates for both characterized transcription factors and orthologues of characterized mammalian transcription factors.

Genome-wide relationship between histone H3 lysine 4 mono- and tri-methylation and transcription factor binding

We characterized the relationship of H3K4me1 and H3K4me3 at distal and proximal regulatory elements by comparing ChIP-seq profiles for these histone modifications and for two functionally different transcription factors: STAT1 in the immortalized HeLa S3 cell line, with and without interferon-gamma (IFNG) stimulation; and FOXA2 in mouse adult liver tissue. In unstimulated and stimulated HeLa cells, respectively, we determined approximately 270,000 and approximately 301,000 H3K4me1-enriched regions, and approximately 54,500 and approximately 76,100 H3K4me3-enriched regions. In mouse adult liver, we determined approximately 227,000 and approximately 34,800 H3K4me1 and H3K4me3 regions. Seventy-five percent of the approximately 70,300 STAT1 binding sites in stimulated HeLa cells and 87% of the approximately 11,000 FOXA2 sites in mouse liver were distal to known gene TSS; in both cell types, approximately 83% of these distal sites were associated with at least one of the two histone modifications, and H3K4me1 was associated with over 96% of marked distal sites. After filtering against predicted transcription start sites, 50% of approximately 26,800 marked distal IFNG-stimulated STAT1 binding sites, but 95% of approximately 5800 marked distal FOXA2 sites, were associated with H3K4me1 only. Results for HeLa cells generated additional insights into transcriptional regulation involving STAT1. STAT1 binding was associated with 25% of all H3K4me1 regions in stimulated HeLa cells, suggesting that a single transcription factor can interact with an unexpectedly large fraction of regulatory regions. Strikingly, for a large majority of the locations of stimulated STAT1 binding, the dominant H3K4me1/me3 combinations were established before activation, suggesting mechanisms independent of IFNG stimulation and high-affinity STAT1 binding.

Global analysis of in vivo Foxa2-binding sites in mouse adult liver using massively parallel sequencing

Foxa2 (HNF3β) is a one of three, closely related transcription factors that are critical to the development and function of the mouse liver. We have used chromatin immunoprecipitation and massively parallel Illumina 1G sequencing (ChIP–Seq) to create a genome-wide profile of in vivo Foxa2-binding sites in the adult liver. More than 65% of the ∼11.5 k genomic sites associated with Foxa2 binding, mapped to extended gene regions of annotated genes, while more than 30% of intragenic sites were located within first introns. 20.5% of all sites were further than 50 kb from any annotated gene, suggesting an association with novel gene regions. QPCR analysis demonstrated a strong positive correlation between peak height and fold enrichment for Foxa2-binding sites. We measured the relationship between Foxa2 and liver gene expression by overlapping Foxa2-binding sites with a SAGE transcriptome profile, and found that 43.5% of genes expressed in the liver were also associated with Foxa2 binding. We also identified potential Foxa2-interacting transcription factors whose motifs were enriched near Foxa2-binding sites. Our comprehensive results for in vivo Foxa2-binding sites in the mouse liver will contribute to resolving transcriptional regulatory networks that are important for adult liver function.