Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Longitudinal molecular trajectories of diffuse glioma in adults.

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear1,2. Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of diffuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specific gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at different rates across the glioma subtypes, and hypermutation was not associated with differences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner.

Quantitative Assessment of Preanalytic Variables on Clinical Evaluation of PI3/AKT/mTOR Signaling Activity in Diffuse Glioma.

Biomarker-driven therapeutic clinical trials require the implementation of standardized, evidence-based practices for sample collection. In diffuse glioma, phosphatidylinositol 3 (PI3)-kinase/AKT/mTOR (PI3/AKT/mTOR) signaling is an attractive therapeutic target for which window-of-opportunity clinical trials could facilitate the identification of promising new agents. Yet, the relevant preanalytic variables and optimal tumor sampling methods necessary to measure pathway activity are unknown. To address this, we used a murine model for isocitrate dehydrogenase (IDH)-wildtype glioblastoma (GBM) and human tumor tissue, including IDH-wildtype GBM and IDH-mutant diffuse glioma. First, we determined the impact of delayed time-to-formalin fixation, or cold ischemia time (CIT), on the quantitative assessment of cellular expression of 6 phosphoproteins that are readouts of PI3K/AK/mTOR activity (phosphorylated-proline-rich Akt substrate of 40 kDa (p-PRAS40, T246), -mechanistic target of rapamycin (p-mTOR; S2448); -AKT (p-AKT, S473); -ribosomal protein S6 (p-RPS6, S240/244 and S235/236), and -eukaryotic initiation factor 4E-binding protein 1 (p-4EBP1, T37/46). With CITs ≥ 2 hours, typical of routine clinical handling, all had reduced or altered expression with p-RPS6 (S240/244) exhibiting relatively greater stability. A similar pattern was observed using patient tumor samples from the operating room with p-4EBP1 more sensitive to delayed fixation than p-RPS6 (S240/244). Many clinical trials utilize unstained slides for biomarker evaluation. Thus, we evaluated the impact of slide storage conditions on the detection of p-RPS6 (S240/244), p-4EBP1, and p-AKT. After 5 months, storage at -80°C was required to preserve the expression of p-4EBP1 and p-AKT, whereas p-RPS6 (240/244) expression was not stable regardless of storage temperature. Biomarker heterogeneity impacts optimal tumor sampling. Quantification of p-RPS6 (240/244) expression in multiple regionally distinct human tumor samples from 8 patients revealed significant intratumoral heterogeneity. Thus, the accurate assessment of PI3K/AKT/mTOR signaling in diffuse glioma must overcome intratumoral heterogeneity and multiple preanalytic factors, including time-to-formalin fixation, slide storage conditions, and phosphoprotein of interest.

Cancer-specific loss of TERT activation sensitizes glioblastoma to DNA damage.

Most glioblastomas (GBMs) achieve cellular immortality by acquiring a mutation in the telomerase reverse transcriptase (TERT) promoter. TERT promoter mutations create a binding site for a GA binding protein (GABP) transcription factor complex, whose assembly at the promoter is associated with TERT reactivation and telomere maintenance. Here, we demonstrate increased binding of a specific GABPB1L-isoform-containing complex to the mutant TERT promoter. Furthermore, we find that TERT promoter mutant GBM cells, unlike wild-type cells, exhibit a critical near-term dependence on GABPB1L for proliferation, notably also posttumor establishment in vivo. Up-regulation of the protein paralogue GABPB2, which is normally expressed at very low levels, can rescue this dependence. More importantly, when combined with frontline temozolomide (TMZ) chemotherapy, inducible GABPB1L knockdown and the associated TERT reduction led to an impaired DNA damage response that resulted in profoundly reduced growth of intracranial GBM tumors. Together, these findings provide insights into the mechanism of cancer-specific TERT regulation, uncover rapid effects of GABPB1L-mediated TERT suppression in GBM maintenance, and establish GABPB1L inhibition in combination with chemotherapy as a therapeutic strategy for TERT promoter mutant GBM.

miRNA-independent function of long noncoding pri-miRNA loci.

Among the large, diverse set of mammalian long noncoding RNAs (lncRNAs), long noncoding primary microRNAs (lnc-pri-miRNAs) are those that host miRNAs. Whether lnc-pri-miRNA loci have important biological function independent of their cognate miRNAs is poorly understood. From a genome-scale lncRNA screen, lnc-pri-miRNA loci were enriched for function in cell proliferation, and in glioblastoma (i.e., GBM) cells with DGCR8 or DROSHA knockdown, lnc-pri-miRNA screen hits still regulated cell growth. To molecularly dissect the function of a lnc-pri-miRNA locus, we studied LOC646329 (also known as MIR29HG), which hosts the miR-29a/b1 cluster. In GBM cells, LOC646329 knockdown reduced miR-29a/b1 levels, and these cells exhibited decreased growth. However, genetic deletion of the miR-29a/b1 cluster (LOC646329-miR29Δ) did not decrease cell growth, while knockdown of LOC646329-miR29Δ transcripts reduced cell proliferation. The miR-29a/b1-independent activity of LOC646329 corresponded to enhancer-like activation of a neighboring oncogene (MKLN1), regulating cell propagation. The LOC646329 locus interacts with the MKLN1 promoter, and antisense oligonucleotide knockdown of the lncRNA disrupts these interactions and reduces the enhancer-like activity. More broadly, analysis of genome-wide data from multiple human cell types showed that lnc-pri-miRNA loci are significantly enriched for DNA looping interactions with gene promoters as well as genomic and epigenetic characteristics of transcriptional enhancers. Functional studies of additional lnc-pri-miRNA loci demonstrated cognate miRNA-independent enhancer-like activity. Together, these data demonstrate that lnc-pri-miRNA loci can regulate cell biology via both miRNA-dependent and miRNA-independent mechanisms.

"De novo replication repair deficient glioblastoma, IDH-wildtype" is a distinct glioblastoma subtype in adults that may benefit from immune checkpoint blockade.

Glioblastoma is a clinically and molecularly heterogeneous disease, and new predictive biomarkers are needed to identify those patients most likely to respond to specific treatments. Through prospective genomic profiling of 459 consecutive primary treatment-naïve IDH-wildtype glioblastomas in adults, we identified a unique subgroup (2%, 9/459) defined by somatic hypermutation and DNA replication repair deficiency due to biallelic inactivation of a canonical mismatch repair gene. The deleterious mutations in mismatch repair genes were often present in the germline in the heterozygous state with somatic inactivation of the remaining allele, consistent with glioblastomas arising due to underlying Lynch syndrome. A subset of tumors had accompanying proofreading domain mutations in the DNA polymerase POLE and resultant "ultrahypermutation". The median age at diagnosis was 50 years (range 27-78), compared with 63 years for the other 450 patients with conventional glioblastoma (p < 0.01). All tumors had histologic features of the giant cell variant of glioblastoma. They lacked EGFR amplification, lacked combined trisomy of chromosome 7 plus monosomy of chromosome 10, and only rarely had TERT promoter mutation or CDKN2A homozygous deletion, which are hallmarks of conventional IDH-wildtype glioblastoma. Instead, they harbored frequent inactivating mutations in TP53, NF1, PTEN, ATRX, and SETD2 and recurrent activating mutations in PDGFRA. DNA methylation profiling revealed they did not align with known reference adult glioblastoma methylation classes, but instead had unique globally hypomethylated epigenomes and mostly classified as "Diffuse pediatric-type high grade glioma, RTK1 subtype, subclass A". Five patients were treated with immune checkpoint blockade, four of whom survived greater than 3 years. The median overall survival was 36.8 months, compared to 15.5 months for the other 450 patients (p < 0.001). We conclude that "De novo replication repair deficient glioblastoma, IDH-wildtype" represents a biologically distinct subtype in the adult population that may benefit from prospective identification and treatment with immune checkpoint blockade.

Multiplatform molecular analyses refine classification of gliomas arising in patients with neurofibromatosis type 1.

Gliomas arising in the setting of neurofibromatosis type 1 (NF1) are heterogeneous, occurring from childhood through adulthood, can be histologically low-grade or high-grade, and follow an indolent or aggressive clinical course. Comprehensive profiling of genetic alterations beyond NF1 inactivation and epigenetic classification of these tumors remain limited. Through next-generation sequencing, copy number analysis, and DNA methylation profiling of gliomas from 47 NF1 patients, we identified 2 molecular subgroups of NF1-associated gliomas. The first harbored biallelic NF1 inactivation only, occurred primarily during childhood, followed a more indolent clinical course, and had a unique epigenetic signature for which we propose the terminology "pilocytic astrocytoma, arising in the setting of NF1". The second subgroup harbored additional oncogenic alterations including CDKN2A homozygous deletion and ATRX mutation, occurred primarily during adulthood, followed a more aggressive clinical course, and was epigenetically diverse, with most tumors aligning with either high-grade astrocytoma with piloid features or various subclasses of IDH-wildtype glioblastoma. Several patients were treated with small molecule MEK inhibitors that resulted in stable disease or tumor regression when used as a single agent, but only in the context of those tumors with NF1 inactivation lacking additional oncogenic alterations. Together, these findings highlight recurrently altered pathways in NF1-associated gliomas and help inform targeted therapeutic strategies for this patient population.

Clonal expansion and epigenetic reprogramming following deletion or amplification of mutant IDH1.

IDH1 mutation is the earliest genetic alteration in low-grade gliomas (LGGs), but its role in tumor recurrence is unclear. Mutant IDH1 drives overproduction of the oncometabolite d-2-hydroxyglutarate (2HG) and a CpG island (CGI) hypermethylation phenotype (G-CIMP). To investigate the role of mutant IDH1 at recurrence, we performed a longitudinal analysis of 50 IDH1 mutant LGGs. We discovered six cases with copy number alterations (CNAs) at the IDH1 locus at recurrence. Deletion or amplification of IDH1 was followed by clonal expansion and recurrence at a higher grade. Successful cultures derived from IDH1 mutant, but not IDH1 wild type, gliomas systematically deleted IDH1 in vitro and in vivo, further suggestive of selection against the heterozygous mutant state as tumors progress. Tumors and cultures with IDH1 CNA had decreased 2HG, maintenance of G-CIMP, and DNA methylation reprogramming outside CGI. Thus, while IDH1 mutation initiates gliomagenesis, in some patients mutant IDH1 and 2HG are not required for later clonal expansions.

The clonal and mutational evolution spectrum of primary triple-negative breast cancers.

Primary triple-negative breast cancers (TNBCs), a tumour type defined by lack of oestrogen receptor, progesterone receptor and ERBB2 gene amplification, represent approximately 16% of all breast cancers. Here we show in 104 TNBC cases that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some having only a handful of coding somatic aberrations in a few pathways, whereas others contain hundreds of coding somatic mutations. High-throughput RNA sequencing (RNA-seq) revealed that only approximately 36% of mutations are expressed. Using deep re-sequencing measurements of allelic abundance for 2,414 somatic mutations, we determine for the first time-to our knowledge-in an epithelial tumour subtype, the relative abundance of clonal frequencies among cases representative of the population. We show that TNBCs vary widely in their clonal frequencies at the time of diagnosis, with the basal subtype of TNBC showing more variation than non-basal TNBC. Although p53 (also known as TP53), PIK3CA and PTEN somatic mutations seem to be clonally dominant compared to other genes, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal, cell shape and motility proteins occurred at lower clonal frequencies, suggesting that they occurred later during tumour progression. Taken together, our results show that understanding the biology and therapeutic responses of patients with TNBC will require the determination of individual tumour clonal genotypes.

Probing the phosphatidylinositol 3-kinase/mammalian target of rapamycin pathway in gliomas: A phase 2 study of everolimus for recurrent adult low-grade gliomas.

BACKGROUND: Activation of the phosphatidylinositol 3-kinase (PI3K)/mammalian target of rapamycin (mTOR) pathway is common in patients with low-grade gliomas (LGGs), but agents that inhibit this pathway, including mTOR inhibitors, have not been studied in this population. METHODS: Fifty-eight patients with pathologic evidence of recurrence after they had initially been diagnosed with World Health Organization (WHO) grade II gliomas were enrolled into a prospective phase 2 clinical trial and received daily everolimus (RAD001) for 1 year or until progression. Tissue at the time of enrollment was analyzed for markers of PI3K/mTOR pathway activation. Thirty-eight patients underwent serial multiparametric magnetic resonance imaging, with the tumor volume and the perfusion metrics (the fractional blood volume [fBV] for capillary density and the transfer coefficient [Kps ] for vascular permeability) measured during treatment. The primary endpoint was progression-free survival at 6 months (PFS-6) in patients with WHO II disease at enrollment. RESULTS: For patients with WHO II gliomas at enrollment, the PFS-6 rate was 84%, and this met the primary endpoint (P < .001 for an improvement from the historical rate of 17%). Evidence of PI3K/mTOR activation by immunohistochemistry for phosphorylated ribosomal S6Ser240/244 (p-S6Ser240/244 ) was associated with worse progression-free survival (PFS; hazard ratio [HR], 3.03; P = .004) and overall survival (HR, 12.7; P = .01). Tumor perfusion decreased after 6 months (median decrease in fBV, 15%; P = .03; median decrease in Kps , 12%; P = .09), with greater decreases associated with improved PFS (HR for each 10% fBV decrease, 0.71; P = .01; HR for each 10% Kps decrease, 0.82; P = .04). CONCLUSIONS: Patients with recurrent LGGs demonstrated a high degree of disease stability during treatment with everolimus. PI3K/mTOR activation, as measured by immunohistochemistry for p-S6, was associated with a worse prognosis. Tumor vascular changes were observed that were consistent with the antiangiogenic effects of mTOR inhibitors. These results support further study of everolimus for LGGs. Cancer 2017;123:4631-4639. © 2017 American Cancer Society.

Recurrent epimutations activate gene body promoters in primary glioblastoma.

Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

Estimating absolute methylation levels at single-CpG resolution from methylation enrichment and restriction enzyme sequencing methods.

Recent advancements in sequencing-based DNA methylation profiling methods provide an unprecedented opportunity to map complete DNA methylomes. These include whole-genome bisulfite sequencing (WGBS, MethylC-seq, or BS-seq), reduced-representation bisulfite sequencing (RRBS), and enrichment-based methods such as MeDIP-seq, MBD-seq, and MRE-seq. These methods yield largely comparable results but differ significantly in extent of genomic CpG coverage, resolution, quantitative accuracy, and cost, at least while using current algorithms to interrogate the data. None of these existing methods provides single-CpG resolution, comprehensive genome-wide coverage, and cost feasibility for a typical laboratory. We introduce methylCRF, a novel conditional random fields-based algorithm that integrates methylated DNA immunoprecipitation (MeDIP-seq) and methylation-sensitive restriction enzyme (MRE-seq) sequencing data to predict DNA methylation levels at single-CpG resolution. Our method is a combined computational and experimental strategy to produce DNA methylomes of all 28 million CpGs in the human genome for a fraction (<10%) of the cost of whole-genome bisulfite sequencing methods. methylCRF was benchmarked for accuracy against Infinium arrays, RRBS, WGBS sequencing, and locus-specific bisulfite sequencing performed on the same human embryonic stem cell line. methylCRF transformation of MeDIP-seq/MRE-seq was equivalent to a biological replicate of WGBS in quantification, coverage, and resolution. We used conventional bisulfite conversion, PCR, cloning, and sequencing to validate loci where our predictions do not agree with whole-genome bisulfite data, and in 11 out of 12 cases, methylCRF predictions of methylation level agree better with validated results than does whole-genome bisulfite sequencing. Therefore, methylCRF transformation of MeDIP-seq/MRE-seq data provides an accurate, inexpensive, and widely accessible strategy to create full DNA methylomes.

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm.

DNA methylation plays key roles in diverse biological processes such as X chromosome inactivation, transposable element repression, genomic imprinting, and tissue-specific gene expression. Sequencing-based DNA methylation profiling provides an unprecedented opportunity to map and compare complete DNA methylomes. This includes one of the most widely applied technologies for measuring DNA methylation: methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq), coupled with a complementary method, methylation-sensitive restriction enzyme sequencing (MRE-seq). A computational approach that integrates data from these two different but complementary assays and predicts methylation differences between samples has been unavailable. Here, we present a novel integrative statistical framework M&M (for integration of MeDIP-seq and MRE-seq) that dynamically scales, normalizes, and combines MeDIP-seq and MRE-seq data to detect differentially methylated regions. Using sample-matched whole-genome bisulfite sequencing (WGBS) as a gold standard, we demonstrate superior accuracy and reproducibility of M&M compared to existing analytical methods for MeDIP-seq data alone. M&M leverages the complementary nature of MeDIP-seq and MRE-seq data to allow rapid comparative analysis between whole methylomes at a fraction of the cost of WGBS. Comprehensive analysis of nineteen human DNA methylomes with M&M reveals distinct DNA methylation patterns among different tissue types, cell types, and individuals, potentially underscoring divergent epigenetic regulation at different scales of phenotypic diversity. We find that differential DNA methylation at enhancer elements, with concurrent changes in histone modifications and transcription factor binding, is common at the cell, tissue, and individual levels, whereas promoter methylation is more prominent in reinforcing fundamental tissue identities.

Conserved role of intragenic DNA methylation in regulating alternative promoters.

Although it is known that the methylation of DNA in 5' promoters suppresses gene expression, the role of DNA methylation in gene bodies is unclear. In mammals, tissue- and cell type-specific methylation is present in a small percentage of 5' CpG island (CGI) promoters, whereas a far greater proportion occurs across gene bodies, coinciding with highly conserved sequences. Tissue-specific intragenic methylation might reduce, or, paradoxically, enhance transcription elongation efficiency. Capped analysis of gene expression (CAGE) experiments also indicate that transcription commonly initiates within and between genes. To investigate the role of intragenic methylation, we generated a map of DNA methylation from the 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 shown to be in intragenic and intergenic regions, whereas 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 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.

Loss of WNK2 expression by promoter gene methylation occurs in adult gliomas and triggers Rac1-mediated tumour cell invasiveness.

The gene encoding protein kinase WNK2 was recently identified to be silenced by promoter hypermethylation in gliomas and meningiomas, suggesting a tumour-suppressor role in these brain tumours. Following experimental depletion in cell lines, WNK2 was further found to control GTP-loading of Rac1, a signalling guanosine triphosphatase involved in cell migration and motility. Here we show that WNK2 promoter methylation also occurs in 17.5% (29 out of 166) of adult gliomas, whereas it is infrequent in its paediatric forms (1.6%; 1 out of 66). Re-expression of WNK2 in glioblastoma cells presenting WNK2 gene silencing reduced cell proliferation in vitro, tumour growth in vivo and also cell migration and invasion, an effect correlated with reduced activation of Rac1. In contrast, when endogenous WNK2 was depleted from glioblastoma cells with unmethylated WNK2 promoter, changes in cell morphology, an increase in invasion and activation of Rac1 were observed. Together, these results validate the WNK2 gene as a recurrent target for epigenetic silencing in glia-derived brain tumours and provide first mechanistic evidence for a tumour-suppressing role of WNK2 that is related to Rac1 signalling and tumour cell invasion and proliferation.

Spark: a navigational paradigm for genomic data exploration.

Biologists possess the detailed knowledge critical for extracting biological insight from genome-wide data resources, and yet they are increasingly faced with nontrivial computational analysis challenges posed by genome-scale methodologies. To lower this computational barrier, particularly in the early data exploration phases, we have developed an interactive pattern discovery and visualization approach, Spark, designed with epigenomic data in mind. Here we demonstrate Spark's ability to reveal both known and novel epigenetic signatures, including a previously unappreciated binding association between the YY1 transcription factor and the corepressor CTBP2 in human embryonic stem cells.

methylC Track: visual integration of single-base resolution DNA methylation data on the WashU EpiGenome Browser.

SUMMARY: We present methylC track, an efficient mechanism for visualizing single-base resolution DNA methylation data on a genome browser. The methylC track dynamically integrates the level of methylation, the position and context of the methylated cytosine (i.e. CG, CHG and CHH), strand and confidence level (e.g. read coverage depth in the case of whole-genome bisulfite sequencing data). Investigators can access and integrate these information visually at specific locus or at the genome-wide level on the WashU EpiGenome Browser in the context of other rich epigenomic datasets. AVAILABILITY AND IMPLEMENTATION: The methylC track is part of the WashU EpiGenome Browser, which is open source and freely available at http://epigenomegateway.wustl.edu/browser/. The most up-to-date instructions and tools for preparing methylC track are available at http://epigenomegateway.wustl.edu/+/cmtk. CONTACT: twang@genetics.wustl.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.