5’isomiR-183-5p|+2 elicits tumor suppressor activity in a negative feedback loop with E2F1

Background

MicroRNAs (miRNAs) and isomiRs play important roles in tumorigenesis as essential regulators of gene expression. 5’isomiRs exhibit a shifted seed sequence compared to the canonical miRNA, resulting in different target spectra and thereby extending the phenotypic impact of the respective common pre-miRNA. However, for most miRNAs, expression and function of 5’isomiRs have not been studied in detail yet. Therefore, this study aims to investigate the functions of miRNAs and their 5’isomiRs.

Methods

The expression of 5’isomiRs was assessed in The Cancer Genome Atlas (TCGA) breast cancer patient dataset. Phenotypic effects of miR-183 overexpression in triple-negative breast cancer (TNBC) cell lines were investigated in vitro and in vivo by quantifying migration, proliferation, tumor growth and metastasis. Direct targeting of E2F1 by miR-183-5p|+2 was validated with a 3’UTR luciferase assay and linked to the phenotypes of isomiR overexpression.

Results

TCGA breast cancer patient data indicated that three variants of miR-183-5p are highly expressed and upregulated, namely miR-183-5p|0, miR-183-5p|+1 and miR-183-5p|+2. However, TNBC cell lines displayed reduced proliferation and invasion upon overexpression of pre-miR-183. While invasion was reduced individually by all three isomiRs, proliferation and cell cycle progression were specifically inhibited by overexpression of miR-183-5p|+2. Proteomic analysis revealed reduced expression of E2F target genes upon overexpression of this isomiR, which could be attributed to direct targeting of E2F1, specifically by miR-183-5p|+2. Knockdown of E2F1 partially phenocopied the effect of miR-183-5p|+2 overexpression on cell proliferation and cell cycle. Gene set enrichment analysis of TCGA and METABRIC patient data indicated that the activity of E2F strongly correlated with the expression of miR-183-5p, suggesting transcriptional regulation of the miRNA by a factor of the E2F family. Indeed, in vitro, expression of miR-183-5p was regulated by E2F1. Hence, miR-183-5p|+2 directly targeting E2F1 appears to be part of a negative feedback loop potentially fine-tuning its activity.

Conclusions

This study demonstrates that 5’isomiRs originating from the same arm of the same pre-miRNA (i.e. pre-miR-183-5p) may exhibit different functions and thereby collectively contribute to the same phenotype. Here, one of three isomiRs was shown to counteract expression of the pre-miRNA by negatively regulating a transcriptional activator (i.e. E2F1). We speculate that this might be part of a regulatory mechanism to prevent uncontrolled cell proliferation, which is disabled during cancer progression.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02380-8.

The antileukemic activity of decitabine upon PML/RARA-negative AML blasts is supported by all-trans retinoic acid: in vitro and in vivo evidence for cooperation

The prognosis of AML patients with adverse genetics, such as a complex, monosomal karyotype and TP53 lesions, is still dismal even with standard chemotherapy. DNA-hypomethylating agent monotherapy induces an encouraging response rate in these patients. When combined with decitabine (DAC), all-trans retinoic acid (ATRA) resulted in an improved response rate and longer overall survival in a randomized phase II trial (DECIDER; NCT00867672). The molecular mechanisms governing this in vivo synergism are unclear. We now demonstrate cooperative antileukemic effects of DAC and ATRA on AML cell lines U937 and MOLM-13. By RNA-sequencing, derepression of >1200 commonly regulated transcripts following the dual treatment was observed. Overall chromatin accessibility (interrogated by ATAC-seq) and, in particular, at motifs of retinoic acid response elements were affected by both single-agent DAC and ATRA, and enhanced by the dual treatment. Cooperativity regarding transcriptional induction and chromatin remodeling was demonstrated by interrogating the HIC1, CYP26A1, GBP4, and LYZ genes, in vivo gene derepression by expression studies on peripheral blood blasts from AML patients receiving DAC + ATRA. The two drugs also cooperated in derepression of transposable elements, more effectively in U937 (mutated TP53) than MOLM-13 (intact TP53), resulting in a “viral mimicry” response. In conclusion, we demonstrate that in vitro and in vivo, the antileukemic and gene-derepressive epigenetic activity of DAC is enhanced by ATRA.

GenomeTornadoPlot: a novel R package for CNV visualization and focality analysis

MOTIVATION

Analysis of focal copy number variations (CNVs) is highly relevant for cancer research, as they pinpoint driver genes. More specifically, due to selective pressure oncogenes and tumor suppressor genes are more often affected by these events than neighboring passengers. In cases where multiple candidates co-reside in a genomic locus, careful comparison is required to either identify multigenic minimally deleted regions of synergistic co-mutations, or the true single driver gene. The study of focal CNVs in large cancer genome cohorts requires specialized visualization and statistical analysis.

RESULTS

We developed the GenomeTornadoPlot R-package which generates gene-centric visualizations of CNV types, locations and lengths from cohortwise NGS data. Furthermore, the software enables the pairwise comparison of proximate genes to identify co-mutation patterns or driver-passenger hierarchies. The visual examination provided by GenomeTornadoPlot is further supported by adaptable local and global focality scoring. Integrated into the GenomeTornadoPlot R-Package is the comprehensive PCAWG database of CNVs, comprising 2976 cancer genome entities from 46 cohorts of the Pan-cancer Analysis of Whole Genomes project. The GenomeTornadoPlot R-package can be used to perform exploratory or hypothesis-driven analyses on the basis of the PCAWG data or in combination with data provided by the user.

AVAILABILITY AND IMPLEMENTATION

GenomeTornadoPlot is written in R script and released via github: <https://github.com/chenhong-dkfz/GenomeTornadoPlot/>. The package is under the license of GPL-3.0.

On the impact of batch effect correction in TCGA isomiR expression data

MicroRNAs (miRNAs) are small non-coding RNAs with diverse functions in post-transcriptional regulation of gene expression. Sequence and length variants of miRNAs are called isomiRs and can exert different functions compared to their canonical counterparts. The Cancer Genome Atlas (TCGA) provides isomiR-level expression data for patients of various cancer entities collected in a multi-center approach over several years. However, the impact of batch effects within individual cohorts has not been systematically investigated and corrected for before. Therefore, the aim of this study was to identify relevant cohort-specific batch variables and generate batch-corrected isomiR expression data for 16 TCGA cohorts. The main batch variables included sequencing platform, plate, sample purity and sequencing depth. Platform bias was related to certain length and sequence features of individual recurrently affected isomiRs. Furthermore, significant downregulation of reported tumor suppressive isomiRs in lung tumor tissue compared to normal samples was only observed after batch correction, highlighting the importance of working with corrected data. Batch-corrected datasets for all cohorts including quality control are provided as supplement. In summary, this study reveals that batch effects present in the TCGA dataset might mask biologically relevant effects and provides a valuable resource for research on isomiRs in cancer (accessible through GEO: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE164767).

Alternative lengthening of telomeres in childhood neuroblastoma from genome to proteome

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

The landscape of chromothripsis across adult cancer types

Chromothripsis is a recently identified mutational phenomenon, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Considered as an early event in tumour development, this form of genome instability plays a prominent role in tumour onset. Chromothripsis prevalence might have been underestimated when using low-resolution methods, and pan-cancer studies based on sequencing are rare. Here we analyse chromothripsis in 28 tumour types covering all major adult cancers (634 tumours, 316 whole-genome and 318 whole-exome sequences). We show that chromothripsis affects a substantial proportion of human cancers, with a prevalence of 49% across all cases. Chromothripsis generates entity-specific genomic alterations driving tumour development, including clinically relevant druggable fusions. Chromothripsis is linked with specific telomere patterns and univocal mutational signatures in distinct tumour entities. Longitudinal analysis of chromothriptic patterns in 24 matched tumour pairs reveals insights in the clonal evolution of tumours with chromothripsis.

The shattering of chromosomes is a dramatic early event in tumourigenesis and is termed chromothripsis. Here, the authors examine chromothripsis across 28 tumour types and show that 49% of cancers exhibit features of chromothripsis.

Genomic footprints of activated telomere maintenance mechanisms in cancer

Cancers require telomere maintenance mechanisms for unlimited replicative potential. They achieve this through TERT activation or alternative telomere lengthening associated with ATRX or DAXX loss. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, we dissect whole-genome sequencing data of over 2500 matched tumor-control samples from 36 different tumor types aggregated within the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium to characterize the genomic footprints of these mechanisms. While the telomere content of tumors with ATRX or DAXX mutations (ATRX/DAXXtrunc) is increased, tumors with TERT modifications show a moderate decrease of telomere content. One quarter of all tumor samples contain somatic integrations of telomeric sequences into non-telomeric DNA. This fraction is increased to 80% prevalence in ATRX/DAXXtrunc tumors, which carry an aberrant telomere variant repeat (TVR) distribution as another genomic marker. The latter feature includes enrichment or depletion of the previously undescribed singleton TVRs TTCGGG and TTTGGG, respectively. Our systematic analysis provides new insight into the recurrent genomic alterations associated with telomere maintenance mechanisms in cancer.

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.

Analyses of non-coding somatic drivers in 2,658 cancer whole genomes

The discovery of drivers of cancer has traditionally focused on protein-coding genes1-4. Here we present analyses of driver point mutations and structural variants in non-coding regions across 2,658 genomes from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium5 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). For point mutations, we developed a statistically rigorous strategy for combining significance levels from multiple methods of driver discovery that overcomes the limitations of individual methods. For structural variants, we present two methods of driver discovery, and identify regions that are significantly affected by recurrent breakpoints and recurrent somatic juxtapositions. Our analyses confirm previously reported drivers6,7, raise doubts about others and identify novel candidates, including point mutations in the 5' region of TP53, in the 3' untranslated regions of NFKBIZ and TOB1, focal deletions in BRD4 and rearrangements in the loci of AKR1C genes. We show that although point mutations and structural variants that drive cancer are less frequent in non-coding genes and regulatory sequences than in protein-coding genes, additional examples of these drivers will be found as more cancer genomes become available.

TelomereHunter - in silico estimation of telomere content and composition from cancer genomes

BACKGROUND

Establishment of telomere maintenance mechanisms is a universal step in tumor development to achieve replicative immortality. These processes leave molecular footprints in cancer genomes in the form of altered telomere content and aberrations in telomere composition. To retrieve these telomere characteristics from high-throughput sequencing data the available computational approaches need to be extended and optimized to fully exploit the information provided by large scale cancer genome data sets.

RESULTS

We here present TelomereHunter, a software for the detailed characterization of telomere maintenance mechanism footprints in the genome. The tool is implemented for the analysis of large cancer genome cohorts and provides a variety of diagnostic diagrams as well as machine-readable output for subsequent analysis. A novel key feature is the extraction of singleton telomere variant repeats, which improves the identification and subclassification of the alternative lengthening of telomeres phenotype. We find that whole genome sequencing-derived telomere content estimates strongly correlate with telomere qPCR measurements (r = 0.94). For the first time, we determine the correlation of in silico telomere content quantification from whole genome sequencing and whole genome bisulfite sequencing data derived from the same tumor sample (r = 0.78). An analogous comparison of whole exome sequencing data and whole genome sequencing data measured slightly lower correlation (r = 0.79). However, this is considerably improved by normalization with matched controls (r = 0.91).

CONCLUSIONS

TelomereHunter provides new functionality for the analysis of the footprints of telomere maintenance mechanisms in cancer genomes. Besides whole genome sequencing, whole exome sequencing and whole genome bisulfite sequencing are suited for in silico telomere content quantification, especially if matched control samples are available. The software runs under a GPL license and is available at https://www.dkfz.de/en/applied-bioinformatics/telomerehunter/telomerehunter.html .

Response to olaparib in a PALB2 germline mutated prostate cancer and genetic events associated with resistance

Prostate cancers harboring DNA repair gene alterations are particularly sensitive to PARP inhibitor treatment. We report a case of an advanced prostate cancer patient profiled within the NCT-MASTER (Molecularly Aided Stratification for Tumor Eradication Research) precision oncology program using next-generation sequencing. Comprehensive genomic and transcriptomic analysis identified a pathogenic germline PALB2 variant as well as a mutational signature associated with disturbed homologous recombination together with structural genomic rearrangements. A molecular tumor board identified a potential benefit of targeted therapy and recommended PARP inhibition and platinum-based chemotherapy. Single-agent treatment with the PARP inhibitor olaparib as well as subsequent combination with platinum-based chemotherapy resulted in disease stabilization and substantial improvement of clinical symptoms. Upon progression, we performed whole-exome and RNA sequencing of a liver metastasis, which demonstrated up-regulation of several genes characteristic for the neuroendocrine prostate cancer phenotype as well as a novel translocation resulting in an in-frame, loss-of-function fusion of RB1. We suggest that multidimensional genomic characterization of prostate cancer patients undergoing PARP inhibitor therapy will be necessary to capture and understand predictive biomarkers of PARP inhibitor sensitivity and resistance.

Genomic and transcriptomic changes complement each other in the pathogenesis of sporadic Burkitt lymphoma

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing.

Burkitt lymphoma (BL) is the most common pediatric B-cell lymphoma. Here, within the International Cancer Genome Consortium, the authors performed whole genome and transcriptome sequencing of 39 sporadic BL, describing the landscape of mutations, structural variants, and mutational processes that underpin this disease how alterations on different cellular levels cooperate in deregulating key pathways and complexes.

Molecular Evolution of Early-Onset Prostate Cancer Identifies Molecular Risk Markers and Clinical Trajectories

Identifying the earliest somatic changes in prostate cancer can give important insights into tumor evolution and aids in stratifying high- from low-risk disease. We integrated whole genome, transcriptome and methylome analysis of early-onset prostate cancers (diagnosis ≤55 years). Characterization across 292 prostate cancer genomes revealed age-related genomic alterations and a clock-like enzymatic-driven mutational process contributing to the earliest mutations in prostate cancer patients. Our integrative analysis identified four molecular subgroups, including a particularly aggressive subgroup with recurrent duplications associated with increased expression of ESRP1, which we validate in 12,000 tissue microarray tumors. Finally, we combined the patterns of molecular co-occurrence and risk-based subgroup information to deconvolve the molecular and clinical trajectories of prostate cancer from single patient samples.

Abstract 5333: Methylation profiling identifies two subclasses of cutaneous squamous cell carcinoma related to distinct cell types of origin

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer type and arises from keratinocytes. Most cSCC progress from a UV-induced precancerous lesion termed actinic keratosis (AK). Despite various efforts to characterize these lesions molecularly, the etiology of AK and its progression to cSCC remain only partially understood. Here we have used Infinium MethylationEPIC BeadChips to interrogate the DNA methylation status of about 850.000 CpGs in epidermis samples from healthy skin, AK and cSCC. Importantly, we found that the premalignant AK samples displayed classical features of cancer methylomes and were highly similar to cSCC methylomes. Further analysis identified typical features of stem cell methylomes, such as a reduced DNA methylation age, non-CpG methylation and stem cell-related keratin and enhancer methylation patterns. Interestingly, this signature was detected only in one half of the AK and cSCC samples, while the other half showed keratin and enhancer methylation patterns that were more closely related to the control epidermis. These findings suggest the existence of two distinct subclasses of AK and cSCC that originate from distinct keratinocyte differentiation stages.

Citation Format: Manuel Rodríguez Paredes, Felix Bormann, Günter Raddatz, Julian Gutekunst, Carlota Lucena-Porcel, Florian Köhler, Elisabeth Wurzer, Katrin Schmidt, Joachim Röwert-Huber, Evgeniya Denisova, Lars Feuerbach, Esther Herpel, Ingo Nindl, Thomas G. Hofmann, Marc Winnefeld, Frank Lyko. Methylation profiling identifies two subclasses of cutaneous squamous cell carcinoma related to distinct cell types of origin [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5333.

NRG1 Fusions in KRAS Wild-Type Pancreatic Cancer

We used whole-genome and transcriptome sequencing to identify clinically actionable genomic alterations in young adults with pancreatic ductal adenocarcinoma (PDAC). Molecular characterization of 17 patients with PDAC enrolled in a precision oncology program revealed gene fusions amenable to pharmacologic inhibition by small-molecule tyrosine kinase inhibitors in all patients with KRAS wild-type (KRAS^WT) tumors (4 of 17). These alterations included recurrent NRG1 rearrangements predicted to drive PDAC development through aberrant ERBB receptor-mediated signaling, and pharmacologic ERBB inhibition resulted in clinical improvement and remission of liver metastases in 2 patients with NRG1-rearranged tumors that had proved resistant to standard treatment. Our findings demonstrate that systematic screening of KRAS^WT tumors for oncogenic fusion genes will substantially improve the therapeutic prospects for a sizeable fraction of patients with PDAC.Significance: Advanced PDAC is a malignancy with few treatment options that lacks molecular mechanism-based therapies. Our study uncovers recurrent gene rearrangements such as NRG1 fusions as disease-driving events in KRAS^wt tumors, thereby providing novel insights into oncogenic signaling and new therapeutic options in this entity. Cancer Discov; 8(9); 1087-95. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

Methylation profiling identifies two subclasses of squamous cell carcinoma related to distinct cells of origin

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer and usually progresses from a UV-induced precancerous lesion termed actinic keratosis (AK). Despite various efforts to characterize these lesions molecularly, the etiology of AK and its progression to cSCC remain partially understood. Here, we use Infinium MethylationEPIC BeadChips to interrogate the DNA methylation status in healthy, AK and cSCC epidermis samples. Importantly, we show that AK methylation patterns already display classical features of cancer methylomes and are highly similar to cSCC profiles. Further analysis identifies typical features of stem cell methylomes, such as reduced DNA methylation age, non-CpG methylation, and stem cell-related keratin and enhancer methylation patterns. Interestingly, this signature is detected only in half of the samples, while the other half shows patterns more closely related to healthy epidermis. These findings suggest the existence of two subclasses of AK and cSCC emerging from distinct keratinocyte differentiation stages.

A comprehensive assessment of somatic mutation detection in cancer using whole-genome sequencing

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼ 100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.

Abstract B42: An integrated view on genetic and epigenetic mechanisms revealed aberrant DNA methylation as an important source for miRNA deregulation in prostate cancer

Prostate cancer (PCa) the most common malignant tumor in males and the third leading cause of cancer-related deaths in Western developed countries. The clinical spectrum of PCa ranges from indolent tumors requiring no therapy to highly aggressive and often metastatic diseases. MicroRNAs (miRNAs) have emerged as important regulators in human tumorigenesis and tumor progression during the last decade. miRNA expression is strongly deregulated in many human cancers including PCa. We investigated the contribution of genetic and epigenetic events to the deregulation of miRNAs in PCa by utilizing the profiling data of 19 Early-Onset Pca (EOPC) from International Cancer Genome Consortium (ICGC) (http://www.icgc.org). For 31% (154/491) of the downregulated miRNAs and for 23% (133/578) of the upregulated miRNAs, DNA hyper- and hypomethylation respectively of a promoter region was observed in at least one of the 13 PCa. Deletions were detected for 26% (126/491) of the downregulated miRNAs in at least one sample. We observed a single translocation event leading to disruption of miR-1244-2 from its putative regulatory region that might contribute to expression decrease. We did not detect any SNV within miRNA precursors. Only one miRNA precursor (miR-4307) was amplified in one patient, but the corresponding mature miRNA was not expressed at all. We demonstrated that, among deletions, methylation of miRNA promoters is one of the major mechanisms leading to miRNA silencing or activation, respectively in PCa. 37 miRNA were downregulated by hypermethylation of regulatory regions and 10 miRNAs were upregulated by hypomethylation in both the ICGC EOPC cohort and the validation dataset on 50 prostate cancer and 48 normal prostate tissues. 82% (9/11) of the miRNAs/miRNA clusters analyzed regulate genes that are involved in PI3K/AKT/PTEN signaling, which represents one of the PCa key pathways being deregulated in more than one third of primary PCa.

Note: This abstract was not presented at the conference.

Citation Format: Olga Bogatyrova, Daniela Wuttig, Lei Gu, Yassen Assenov, Ruprecht Kuner, Constance Baer, Lars Feuerbach, Clarissa Gerhäuser, Dieter Weichenhan, Thorsten Schlomm, Ronald Simon, Guido Sauter, Holger Sültmann. An integrated view on genetic and epigenetic mechanisms revealed aberrant DNA methylation as an important source for miRNA deregulation in prostate cancer. [abstract]. In: Proceedings of the AACR Special Conference on Chromatin and Epigenetics in Cancer; Jun 19-22, 2013; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2013;73(13 Suppl):Abstract nr B42.

DNA methylation patterns facilitate the identification of microRNA transcription start sites: a brain-specific study

Predicting the transcription start sites (TSSs) of microRNAs (miRNAs) is important for understanding how these small RNA molecules, known to regulate translation and stability of protein-coding genes, are regulated themselves. Previous approaches are primarily based on genetic features, trained on TSSs of protein-coding genes, and have low prediction accuracy. Recently, a support vector machine based technique has been proposed for miRNA TSS prediction that uses known miRNA TSS for training the classifier along with a set of existing and novel CpG island based features. Current progress in epigenetics research has provided genomewide and tissue-specific reports about various phenotypic traits. We hypothesize that incorporating epigenetic characteristics into statistical models may lead to better prediction of primary transcripts of human miRNAs. In this paper, we have tested our hypothesis on brain-specific miRNAs by using epigenetic as well as genetic features to predict the primary transcripts. For this, we have used a sophisticated feature selection technique and a robust classification model. Our prediction model achieves an accuracy of more than 80% and establishes the potential of epigenetic analysis for in silico prediction of TSSs.

Integrative genomic analyses reveal an androgen-driven somatic alteration landscape in early-onset prostate cancer

Early-onset prostate cancer (EO-PCA) represents the earliest clinical manifestation of prostate cancer. To compare the genomic alteration landscapes of EO-PCA with "classical" (elderly-onset) PCA, we performed deep sequencing-based genomics analyses in 11 tumors diagnosed at young age, and pursued comparative assessments with seven elderly-onset PCA genomes. Remarkable age-related differences in structural rearrangement (SR) formation became evident, suggesting distinct disease pathomechanisms. Whereas EO-PCAs harbored a prevalence of balanced SRs, with a specific abundance of androgen-regulated ETS gene fusions including TMPRSS2:ERG, elderly-onset PCAs displayed primarily non-androgen-associated SRs. Data from a validation cohort of > 10,000 patients showed age-dependent androgen receptor levels and a prevalence of SRs affecting androgen-regulated genes, further substantiating the activity of a characteristic "androgen-type" pathomechanism in EO-PCA.

Analyzing epigenome data in context of genome evolution and human diseases

This chapter describes bioinformatic tools for analyzing epigenome differences between species and in diseased versus normal cells. We illustrate the interplay of several Web-based tools in a case study of CpG island evolution between human and mouse. Starting from a list of orthologous genes, we use the Galaxy Web service to obtain gene coordinates for both species. These data are further analyzed in EpiGRAPH, a Web-based tool that identifies statistically significant epigenetic differences between genome region sets. Finally, we outline how the use of the statistical programming language R enables deeper insights into the epigenetics of human diseases, which are difficult to obtain without writing custom scripts. In summary, our tutorial describes how Web-based tools provide an easy entry into epigenome data analysis while also highlighting the benefits of learning a scripting language in order to unlock the vast potential of public epigenome datasets.

BiQ Analyzer HT: locus-specific analysis of DNA methylation by high-throughput bisulfite sequencing

Bisulfite sequencing is a widely used method for measuring DNA methylation in eukaryotic genomes. The assay provides single-base pair resolution and, given sufficient sequencing depth, its quantitative accuracy is excellent. High-throughput sequencing of bisulfite-converted DNA can be applied either genome wide or targeted to a defined set of genomic loci (e.g. using locus-specific PCR primers or DNA capture probes). Here, we describe BiQ Analyzer HT (http://biq-analyzer-ht.bioinf.mpi-inf.mpg.de/), a user-friendly software tool that supports locus-specific analysis and visualization of high-throughput bisulfite sequencing data. The software facilitates the shift from time-consuming clonal bisulfite sequencing to the more quantitative and cost-efficient use of high-throughput sequencing for studying locus-specific DNA methylation patterns. In addition, it is useful for locus-specific visualization of genome-wide bisulfite sequencing data.

Reconstructing the ancestral germ line methylation state of young repeats

One of the key objectives of comparative genomics is the characterization of the forces that shape genomes over the course of evolution. In the last decades, evidence has been accumulated that for vertebrate genomes also epigenetic modifications have to be considered in this context. Especially, the elevated mutation frequency of 5-methylcytosine (5mC) is assumed to facilitate the depletion of CpG dinucleotides in species that exhibit global DNA methylation. For instance, the underrepresentation of CpG dinucleotides in many mammalian genomes is attributed to this effect, which is only neutralized in so-called CpG islands (CGIs) that are preferentially unmethylated and thus partially protected from rapid CpG decay. For primate-specific CpG-rich transposable elements from the ALU family, it is unclear whether their elevated CpG frequency is caused by their small age or by the absence of DNA methylation. In consequence, these elements are often misclassified in CGI annotations. We present a method for the estimation of germ line methylation from pairwise ancestral-descendant alignments. The approach is validated in a simulation study and tested on DNA repeats from the AluSx family. We conclude that a predicted unmethylated state in the germ line is highly correlated with epigenetic activity of the respective genomic region. Thus, CpG-rich repeats can be facilitated as in silico probes for the epigenetic potential of their genomic neighborhood.