BISCUIT: an efficient, standards-compliant tool suite for simultaneous genetic and epigenetic inference in bulk and single-cell studies.

Data from both bulk and single-cell whole-genome DNA methylation experiments are under-utilized in many ways. This is attributable to inefficient mapping of methylation sequencing reads, routinely discarded genetic information, and neglected read-level epigenetic and genetic linkage information. We introduce the BISulfite-seq Command line User Interface Toolkit (BISCUIT) and its companion R/Bioconductor package, biscuiteer, for simultaneous extraction of genetic and epigenetic information from bulk and single-cell DNA methylation sequencing. BISCUIT's performance, flexibility and standards-compliant output allow large, complex experimental designs to be characterized on clinical timescales. BISCUIT is particularly suited for processing data from single-cell DNA methylation assays, with its excellent scalability, efficiency, and ability to greatly enhance mappability, a key challenge for single-cell studies. We also introduce the epiBED format for single-molecule analysis of coupled epigenetic and genetic information, facilitating the study of cellular and tissue heterogeneity from DNA methylation sequencing.

Comprehensive molecular and clinical characterization of NUP98 fusions in pediatric acute myeloid leukemia.

NUP98 fusions comprise a family of rare recurrent alterations in AML, associated with adverse outcomes. In order to define the underlying biology and clinical implications of this family of fusions, we performed comprehensive transcriptome, epigenome, and immunophenotypic profiling of 2,235 children and young adults with AML and identified 160 NUP98 rearrangements (7.2%), including 108 NUP98-NSD1 (4.8%), 32 NUP98-KDM5A (1.4%) and 20 NUP98-X cases (0.9%) with 13 different fusion partners. Fusion partners defined disease characteristics and biology; patients with NUP98-NSD1 or NUP98-KDM5A had distinct immunophenotypic, transcriptomic, and epigenomic profiles. Unlike the two most prevalent NUP98 fusions, NUP98-X variants are typically not cryptic. Furthermore, NUP98-X cases are associated with WT1 mutations, and have epigenomic profiles that resemble either NUP98-NSD1 or NUP98-KDM5A. Cooperating FLT3-ITD and WT1 mutations define NUP98-NSD1, and chromosome 13 aberrations are highly enriched in NUP98-KDM5A. Importantly, we demonstrate that NUP98 fusions portend dismal overall survival, with the noteworthy exception of patients bearing abnormal chromosome 13 (clinicaltrials gov. Identifiers: NCT00002798, NCT00070174, NCT00372593, NCT01371981).

Integrated DNA Copy Number and Expression Profiling Identifies IGF1R as a Prognostic Biomarker in Pediatric Osteosarcoma.

Osteosarcoma is a primary malignant bone tumor arising from bone-forming mesenchymal cells in children and adolescents. Despite efforts to understand the biology of the disease and identify novel therapeutics, the survival of osteosarcoma patients remains dismal. We have concurrently profiled the copy number and gene expression of 226 osteosarcoma samples as part of the Strategic Partnering to Evaluate Cancer Signatures (SPECS) initiative. Our results demonstrate the heterogeneous landscape of osteosarcoma in younger populations by showing the presence of genome-wide copy number abnormalities occurring both recurrently among samples and in a high frequency. Insulin growth factor receptor 1 (IGF1R) is a receptor tyrosine kinase which binds IGF1 and IGF2 to activate downstream pathways involved in cell apoptosis and proliferation. We identify prevalent amplification of IGF1R corresponding with increased gene expression in patients with poor survival outcomes. Our results substantiate previously tenuously associated copy number abnormalities identified in smaller datasets (13q34+, 20p13+, 4q35-, 20q13.33-), and indicate the significance of high fibroblast growth factor receptor 2 (FGFR2) expression in distinguishing patients with poor prognosis. FGFR2 is involved in cellular proliferation processes such as division, growth and angiogenesis. In summary, our findings demonstrate the prognostic significance of several genes associated with osteosarcoma pathogenesis.

Non-Invasive Profiling of Advanced Prostate Cancer via Multi-Parametric Liquid Biopsy and Radiomic Analysis.

Integrating liquid biopsies of circulating tumor cells (CTCs) and cell-free DNA (cfDNA) with other minimally invasive measures may yield more comprehensive disease profiles. We evaluated the feasibility of concurrent cellular and molecular analysis of CTCs and cfDNA combined with radiomic analysis of CT scans from patients with metastatic castration-resistant PC (mCRPC). CTCs from 22 patients were enumerated, stained for PC-relevant markers, and clustered based on morphometric and immunofluorescent features using machine learning. DNA from single CTCs, matched cfDNA, and buffy coats was sequenced using a targeted amplicon cancer hotspot panel. Radiomic analysis was performed on bone metastases identified on CT scans from the same patients. CTCs were detected in 77% of patients and clustered reproducibly. cfDNA sequencing had high sensitivity (98.8%) for germline variants compared to WBC. Shared and unique somatic variants in PC-related genes were detected in cfDNA in 45% of patients (MAF > 0.1%) and in CTCs in 92% of patients (MAF > 10%). Radiomic analysis identified a signature that strongly correlated with CTC count and plasma cfDNA level. Integration of cellular, molecular, and radiomic data in a multi-parametric approach is feasible, yielding complementary profiles that may enable more comprehensive non-invasive disease modeling and prediction.

Integrative genomic analysis of matched primary and metastatic pediatric osteosarcoma.

Despite being the most common childhood bone tumor, the genomic characterization of osteosarcoma remains incomplete. In particular, very few osteosarcoma metastases have been sequenced to date, critical to better understand mechanisms of progression and evolution in this tumor. We performed an integrated whole genome and exome sequencing analysis of paired primary and metastatic pediatric osteosarcoma specimens to identify recurrent genomic alterations. Sequencing of 13 osteosarcoma patients including 13 primary, 10 metastatic, and 3 locally recurring tumors revealed a highly heterogeneous mutational landscape, including cases of hypermutation and microsatellite instability positivity, but with virtually no recurrent alterations except for mutations involving the tumor suppressor genes RB1 and TP53. At the germline level, we detected alterations in multiple cancer related genes in the majority of the cohort, including those potentially disrupting DNA damage response pathways. Metastases retained only a minimal number of short variants from their corresponding primary tumors, while copy number alterations showed higher conservation. One recurrently amplified gene, KDR, was highly expressed in advanced cases and associated with poor prognosis. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

SeSAMe: reducing artifactual detection of DNA methylation by Infinium BeadChips in genomic deletions.

We report a new class of artifacts in DNA methylation measurements from Illumina HumanMethylation450 and MethylationEPIC arrays. These artifacts reflect failed hybridization to target DNA, often due to germline or somatic deletions and manifest as incorrectly reported intermediate methylation. The artifacts often survive existing preprocessing pipelines, masquerade as epigenetic alterations and can confound discoveries in epigenome-wide association studies and studies of methylation-quantitative trait loci. We implement a solution, P-value with out-of-band (OOB) array hybridization (pOOBAH), in the R package SeSAMe. Our method effectively masks deleted and hyperpolymorphic regions, reducing or eliminating spurious reports of epigenetic silencing at oft-deleted tumor suppressor genes such as CDKN2A and RB1 in cases with somatic deletions. Furthermore, our method substantially decreases technical variation whilst retaining biological variation, both within and across HM450 and EPIC platform measurements. SeSAMe provides a light-weight, modular DNA methylation data analysis suite, with a performant implementation suitable for efficient analysis of thousands of samples.

GLI3 repressor determines Hedgehog pathway activation and is required for response to SMO antagonist glasdegib in AML.

The Hedgehog (Hh) signaling pathway is activated in many cancers and is a promising target for therapeutic development. Deletions in the receptor Patched (PTCH) or activating mutations in Smoothened (SMO) have been reported in basal cell carcinoma and medulloblastoma, but are largely absent in most tumor types. Therefore, the mechanism of pathway activation in most cancers, including hematological malignancies, remains unknown. In normal tissues, Hh pathway activation via PTCH/SMO causes an increase in the downstream transcriptional activator GLI1 and a decrease in the GLI3 transcriptional repressor (GLI3R). In this article, we confirm that the Hh pathway is active in acute myeloid leukemia (AML), however, this activity is largely independent of SMO. Epigenetic and gene expression analysis of The Cancer Genome Atlas AML data set reveals that GLI3 expression is silenced in most AML patient samples, and the GLI3 locus is abnormally methylated. We show that GLI3R is required for the therapeutic effect of SMO antagonists in AML samples and restoration of GLI3R suppresses the growth of AML. We additionally demonstrate that GLI3R represses AML growth by downregulating AKT expression. In summary, this study provides the first evidence that GLI3R plays an essential role in SMO-independent Hh signaling in AML, and suggests that GLI3R could serve as a potential biomarker for patient selection in SMO antagonist clinical trials. Furthermore, these data support rational combinations of hypomethylating agents with SMO antagonists in clinical trials.

Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi.

Summary: The minfi package is widely used for analyzing Illumina DNA methylation array data. Here we describe modifications to the minfi package required to support the HumanMethylationEPIC ('EPIC') array from Illumina. We discuss methods for the joint analysis and normalization of data from the HumanMethylation450 ('450k') and EPIC platforms. We introduce the single-sample Noob ( ssNoob ) method, a normalization procedure suitable for incremental preprocessing of individual methylation arrays and conclude that this method should be used when integrating data from multiple generations of Infinium methylation arrays. We show how to use reference 450k datasets to estimate cell type composition of samples on EPIC arrays. The cumulative effect of these updates is to ensure that minfi provides the tools to best integrate existing and forthcoming Illumina methylation array data. Availability and Implementation: The minfi package version 1.19.12 or higher is available for all platforms from the Bioconductor project. Contact: khansen@jhsph.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Reprogramming of the human intestinal epigenome by surgical tissue transposition.

Extracellular cues play critical roles in the establishment of the epigenome during development and may also contribute to epigenetic perturbations found in disease states. The direct role of the local tissue environment on the post-development human epigenome, however, remains unclear due to limitations in studies of human subjects. Here, we use an isogenic human ileal neobladder surgical model and compare global DNA methylation levels of intestinal epithelial cells pre- and post-neobladder construction using the Infinium HumanMethylation450 BeadChip. Our study is the first to quantify the effect of environmental cues on the human epigenome and show that the local tissue environment directly modulates DNA methylation patterns in normal differentiated cells in vivo. In the neobladder, the intestinal epithelial cells lose their tissue-specific epigenetic landscape in a time-dependent manner following the tissue's exposure to a bladder environment. We find that de novo methylation of many intestine-specific enhancers occurs at the rate of 0.41% per month (P < 0.01, Pearson = 0.71), while demethylation of primarily non-intestine-specific transcribed regions occurs at the rate of -0.37% per month (P < 0.01, Pearson = -0.57). The dynamic resetting of the DNA methylome in the neobladder not only implicates local environmental cues in the shaping and maintenance of the epigenome but also illustrates an unexpected cross-talk between the epigenome and the cellular environment.

The genomic landscape of the Ewing Sarcoma family of tumors reveals recurrent STAG2 mutation.

The Ewing sarcoma family of tumors (EFT) is a group of highly malignant small round blue cell tumors occurring in children and young adults. We report here the largest genomic survey to date of 101 EFT (65 tumors and 36 cell lines). Using a combination of whole genome sequencing and targeted sequencing approaches, we discover that EFT has a very low mutational burden (0.15 mutations/Mb) but frequent deleterious mutations in the cohesin complex subunit STAG2 (21.5% tumors, 44.4% cell lines), homozygous deletion of CDKN2A (13.8% and 50%) and mutations of TP53 (6.2% and 71.9%). We additionally note an increased prevalence of the BRCA2 K3326X polymorphism in EFT patient samples (7.3%) compared to population data (OR 7.1, p = 0.006). Using whole transcriptome sequencing, we find that 11% of tumors pathologically diagnosed as EFT lack a typical EWSR1 fusion oncogene and that these tumors do not have a characteristic Ewing sarcoma gene expression signature. We identify samples harboring novel fusion genes including FUS-NCATc2 and CIC-FOXO4 that may represent distinct small round blue cell tumor variants. In an independent EFT tissue microarray cohort, we show that STAG2 loss as detected by immunohistochemistry may be associated with more advanced disease (p = 0.15) and a modest decrease in overall survival (p = 0.10). These results significantly advance our understanding of the genomic and molecular underpinnings of Ewing sarcoma and provide a foundation towards further efforts to improve diagnosis, prognosis, and precision therapeutics testing.

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

Insulin-like Growth Factor 2 Gene Expression Molecularly Differentiates Pleuropulmonary Blastoma and Embryonal Rhabdomyosarcoma.

The sarcomatous element in pleuropulmonary blastoma (PPB) is often histologically indistinguishable from embryonal rhabdomyosarcoma (ERMS). A diagnosis of PPB is often made after definitive surgical resection based on pathologic features, most notably the presence of hamartomatous pulmonary elements. Samples from seven PPB patients were obtained from the rhabdomyosarcomatous portion of the tumor by macrodissection. Representative ERMS tumor tissue was selected from 21 ERMS patient samples. Formalin-fixed paraffin-embedded tissue scrolls from each sample were analyzed using the Affymetrix Human Exon arrays. All PPB patients and 7 of 21 ERMS patients were 3 years old and younger. Twenty transcripts (10 annotated, 10 noncoding RNAs) were significantly differentially expressed in ERMS when compared with PPB samples. Insulin-like growth factor 2 (IGF2) was uniformly overexpressed in ERMS (19/21>400) but was expressed at low levels in PPB (P<0.001). Two ERMS cases that had low level IGF2 expression were 3 years and younger of age. No other differences between the 2 approached this degree of significance, despite a common rhabdomyogenic phenotype in the sarcomatous areas of PPB. PPB, unlike most ERMS, appears not to be driven by autocrine IGF2 signaling.

Low-level processing of Illumina Infinium DNA Methylation BeadArrays.

We propose a novel approach to background correction for Infinium HumanMethylation data to account for technical variation in background fluorescence signal. Our approach capitalizes on a new use for the Infinium I design bead types to measure non-specific fluorescence in the colour channel opposite of their design (Cy3/Cy5). This provides tens of thousands of features for measuring background instead of the much smaller number of negative control probes on the platforms (n = 32 for HumanMethylation27 and n = 614 for HumanMethylation450, respectively). We compare the performance of our methods with existing approaches, using technical replicates of both mixture samples and biological samples, and demonstrate that within- and between-platform artefacts can be substantially reduced, with concomitant improvement in sensitivity, by the proposed methods.

Genome-wide characterization of circulating tumor cells identifies novel prognostic genomic alterations in systemic melanoma metastasis.

BACKGROUND: Circulating tumor cells (CTC) have been found in patients with metastatic melanoma and are associated with advanced melanoma stage and poor patient outcome. We hypothesize that CTC harbor genomic changes critical in the development of distant systemic metastasis. Here, we present the first genome-wide copy-number aberration (CNA) and loss of heterozygosity (LOH)-based characterization of melanoma CTC. METHODS: CTC were isolated from peripheral blood monocytes of 13 melanoma patients with regional metastasis stage IIIB/C using antibodies against melanoma-associated cell surface gangliosides. RESULTS: We characterized 251 CNA in CTC. Comparative analysis demonstrated >90% concordance in single-nucleotide polymorphism profiles between paired CTC and tumor metastases. In particular, there were notable recurring CNA across patients. In exploratory studies, the presence of several top CTC-associated CNA was verified in distant metastasis (stage IV) from 27 patients, suggesting that certain genomic changes are propagated from regional metastasis to CTC and to distant systemic metastases. Lastly, an exploratory biomarker panel derived from 5 CTC-associated CNA [CSMD2 (CUB and Sushi multiple domains 2), 1p35.1; CNTNAP5 (contactin associated protein-like 5), 2q14.3; NRDE2 (NRDE-2, necessary for RNA interference, domain containing), 14q32.11; ADAM6 (ADAM metallopeptidase domain 6, pseudogene), 14q32.33; and TRPM2 (transient receptor potential cation channel, subfamily m, member 2), 21q22.3] conferred prognostic utility for melanoma recurrence [hazard ratio (HR), 1.14; CI, 1.00-1.44; P = 0.0471] and death (HR, 2.86; CI, 1.23-14.42; P = 0.0014) in 35 patients with stage IIIB/C melanoma, with a 5-year disease-free survival of 13% vs 69% (P = 0.0006) and overall survival of 28% vs 94% between high-risk and low-risk groups defined by the biomarker panel, respectively. CONCLUSIONS: This study provides the first detailed CNA-based profile of melanoma CTC and illustrates how CTC may be used as a novel approach for identification of systemic metastasis.

EGFR-dependent endocytosis of Wnt9a and Fzd9b promotes ß-catenin signaling during hematopoietic stem cell development in zebrafish.

Cell-to-cell communication through secreted Wnt ligands that bind to members of the Frizzled (Fzd) family of transmembrane receptors is critical for development and homeostasis. Wnt9a signals through Fzd9b, the co-receptor LRP5 or LRP6 (LRP5/6), and the epidermal growth factor receptor (EGFR) to promote early proliferation of zebrafish and human hematopoietic stem cells during development. Here, we developed fluorescently labeled, biologically active Wnt9a and Fzd9b fusion proteins to demonstrate that EGFR-dependent endocytosis of the ligand-receptor complex was required for signaling. In human cells, the Wnt9a-Fzd9b complex was rapidly endocytosed and trafficked through early and late endosomes, lysosomes, and the endoplasmic reticulum. Using small-molecule inhibitors and genetic and knockdown approaches, we found that Wnt9a-Fzd9b endocytosis required EGFR-mediated phosphorylation of the Fzd9b tail, caveolin, and the scaffolding protein EGFR protein substrate 15 (EPS15). LRP5/6 and the downstream signaling component AXIN were required for Wnt9a-Fzd9b signaling but not for endocytosis. Knockdown or loss of EPS15 impaired hematopoietic stem cell development in zebrafish. Other Wnt ligands do not require endocytosis for signaling activity, implying that specific modes of endocytosis and trafficking may represent a method by which Wnt-Fzd specificity is established.

EZH2-driven immune evasion defines high-risk pediatric AML with t(16;21) FUS::ERG gene fusion.

Despite decades of research, acute myeloid leukemia (AML) remains a remarkably lethal malignancy. While pediatric AML (pAML) carries a more favorable prognosis than adult AML, the past 25 years of large clinical trials have produced few improvements in pAML survival. Nowhere is this more evident than in patients carrying a t(16;21)(p11;q22) translocation, which yields the FUS::ERG fusion transcript. Patients with FUS::ERG-positive AML are often primary refractory, and most responders quickly relapse. In COG clinical trials, allogeneic stem cell transplantation was of no benefit to FUS::ERG pAML patients; 100% of transplanted patients succumbed to their disease. Expression of major histocompatibility complex (MHC) class I & II and costimulatory molecules is absent at diagnosis in FUS::ERG AML, mirroring the epigenetic mechanism of post-transplant relapse seen in adult AML and its associated dismal outcomes. Here we show that this class-defining immune-repressive phenotype is driven by overexpression of the EZH2 histone lysine methyltransferase in vitro and in multiple clinical cohorts. We show that treatment with the FDA-approved EZH2 inhibitor tazemetostat along with IFN-γ reverses this phenotype, re-establishes MHC presentation, and severely impairs the viability of FUS::ERG AML cells. EZH2 inhibitors may thus provide the first targeted therapeutic option for patients with this high-risk subtype of pAML, with particular benefit as a bridge to successful allogeneic stem cell transplantation.

T lymphocyte abnormalities in juvenile systemic sclerosis patients.

Multi-center evaluations of pediatric patients with juvenile systemic sclerosis (jSSc) have suggested that the pathogenesis of jSSc may differ from that of systemic sclerosis (SSc) in adult patients. Therefore, we undertook to identify abnormalities in the T lymphocytes of jSSc patients and to determine if they differed from the abnormalities reported in the T lymphocytes of adult SSc patients. We identified decreases in the frequency of resting regulatory T lymphocytes and an increased frequency of CD45RA expressing effector memory (EMRA) CD4 T lymphocytes, which were characterized by an increased frequency of CCR7 protein expressing cells. Neither the increases in the EMRA subpopulation nor the increased CCR7 protein expression have been reported in adult SSc patients. The decrease in resting regulatory T lymphocytes in jSSc patients may permit the expansion of the disease initiating CD4 T lymphocytes present in the CCR7 expressing EMRA CD4 T lymphocyte subpopulation.