Stefan T. Börno

Genome-wide DNA Methylation Events in TMPRSS2–ERG Fusion-Negative Prostate Cancers Implicate an EZH2-Dependent Mechanism with miR-26a Hypermethylation

UNLABELLED Prostate cancer is the second most common cancer among men worldwide. Alterations in the DNA methylation pattern can be one of the leading causes for prostate cancer formation. This study is the first high-throughput sequencing study investigating genome-wide DNA methylation patterns in a large cohort of 51 tumor and 53 benign prostate samples using methylated DNA immunoprecipitation sequencing. Comparative analyses identified more than 147,000 cancer-associated epigenetic alterations. In addition, global methylation patterns show significant differences based on the TMPRSS2-ERG rearrangement status. We propose the hypermethylation of miR-26a as an alternative pathway of ERG rearrangement-independent EZH2 activation. The observed increase in differential methylation events in fusion-negative tumors can explain the tumorigenic process in the absence of genomic rearrangements. SIGNIFICANCE In contrast to TMPRSS2-ERG -rearranged tumors, the pathomechanism for gene fusion-negative tumors is completely unclear. Using a sequencing-based approach, our work uncovers significant global epigenetic alterations in TMPRSS2-ERG gene fusion-negative tumors and provides a mechanistic explanation for the tumor formation process.

High-Throughput miRNA and mRNA Sequencing of Paired Colorectal Normal, Tumor and Metastasis Tissues and Bioinformatic Modeling of miRNA-1 Therapeutic Applications

MiRNAs are discussed as diagnostic and therapeutic molecules. However, effective miRNA drug treatments with miRNAs are, so far, hampered by the complexity of the miRNA networks. To identify potential miRNA drugs in colorectal cancer, we profiled miRNA and mRNA expression in matching normal, tumor and metastasis tissues of eight patients by Illumina sequencing. We validated six miRNAs in a large tissue screen containing 16 additional tumor entities and identified miRNA-1, miRNA-129, miRNA-497 and miRNA-215 as constantly de-regulated within the majority of cancers. Of these, we investigated miRNA-1 as representative in a systems-biology simulation of cellular cancer models implemented in PyBioS and assessed the effects of depletion as well as overexpression in terms of miRNA-1 as a potential treatment option. In this system, miRNA-1 treatment reverted the disease phenotype with different effectiveness among the patients. Scoring the gene expression changes obtained through mRNA-Seq from the same patients we show that the combination of deep sequencing and systems biological modeling can help to identify patient-specific responses to miRNA treatments. We present this data as guideline for future pre-clinical assessments of new and personalized therapeutic options.

The bromodomain protein BRD4 regulates the KEAP1/NRF2-dependent oxidative stress response

The epigenetic sensor BRD4 (bromodomain protein 4) is a potent target for anti-cancer therapies. To study the transcriptional impact of BRD4 in cancer, we generated an expression signature of BRD4 knockdown cells and found oxidative stress response genes significantly enriched. We integrated the RNA-Seq results with DNA-binding sites of BRD4 generated by chromatin immunoprecipitations, correlated these with gene expressions from human prostate cancers and identified 21 top BRD4 candidate genes among which the oxidative stress pathway genes KEAP1, SESN3 and HDAC6 are represented. Knock down of BRD4 or treatment with the BRD4 inhibitor JQ1 resulted in decreased reactive oxygen species (ROS) production and increased cell viability under H2O2 exposure. Consistently, a deregulation of BRD4 diminished the KEAP1/NRF2 axis and led to a disturbed regulation of the inducible heme oxygenase 1 (HMOX1). Without exogenous stress induction, we also found BRD4 directly targeting the HMOX1 promoter over the SP1-binding sites. Our findings provide insight into the transcriptional regulatory network of BRD4 and highlight BRD4 as signal transducer of the cellular response to oxidative stress.

The metabolic background is a global player in Saccharomyces gene expression epistasis

The regulation of gene expression in response to nutrient availability is fundamental to the genotype–phenotype relationship. The metabolic–genetic make-up of the cell, as reflected in auxotrophy, is hence likely to be a determinant of gene expression. Here, we address the importance of the metabolic–genetic background by monitoring transcriptome, proteome and metabolome in a repertoire of 16 Saccharomyces cerevisiae laboratory backgrounds, combinatorially perturbed in histidine, leucine, methionine and uracil biosynthesis. The metabolic background affected up to 85% of the coding genome. Suggesting widespread confounding, these transcriptional changes show, on average, 83% overlap between unrelated auxotrophs and 35% with previously published transcriptomes generated for non-metabolic gene knockouts. Background-dependent gene expression correlated with metabolic flux and acted, predominantly through masking or suppression, on 88% of transcriptional interactions epistatically. As a consequence, the deletion of the same metabolic gene in a different background could provoke an entirely different transcriptional response. Propagating to the proteome and scaling up at the metabolome, metabolic background dependencies reveal the prevalence of metabolism-dependent epistasis at all regulatory levels. Urging a fundamental change of the prevailing laboratory practice of using auxotrophs and nutrient supplemented media, these results reveal epistatic intertwining of metabolism with gene expression on the genomic scale.

ERG Induces Epigenetic Activation of Tudor Domain-Containing Protein 1 (TDRD1) in ERG Rearrangement-Positive Prostate Cancer

Background Overexpression of ERG transcription factor due to genomic ERG-rearrangements defines a separate molecular subtype of prostate tumors. One of the consequences of ERG accumulation is modulation of the cell’s gene expression profile. Tudor domain-containing protein 1 gene (TDRD1) was reported to be differentially expressed between TMPRSS2:ERG-negative and TMPRSS2:ERG-positive prostate cancer. The aim of our study was to provide a mechanistic explanation for the transcriptional activation of TDRD1 in ERG rearrangement-positive prostate tumors. Methodology/Principal Findings Gene expression measurements by real-time quantitative PCR revealed a remarkable co-expression of TDRD1 and ERG (r2 = 0.77) but not ETV1 (r2<0.01) in human prostate cancer in vivo. DNA methylation analysis by MeDIP-Seq and bisulfite sequencing showed that TDRD1 expression is inversely correlated with DNA methylation at the TDRD1 promoter in vitro and in vivo (ρ = −0.57). Accordingly, demethylation of the TDRD1 promoter in TMPRSS2:ERG-negative prostate cancer cells by DNA methyltransferase inhibitors resulted in TDRD1 induction. By manipulation of ERG dosage through gene silencing and forced expression we show that ERG governs loss of DNA methylation at the TDRD1 promoter-associated CpG island, leading to TDRD1 overexpression. Conclusions/Significance We demonstrate that ERG is capable of disrupting a tissue-specific DNA methylation pattern at the TDRD1 promoter. As a result, TDRD1 becomes transcriptionally activated in TMPRSS2:ERG-positive prostate cancer. Given the prevalence of ERG fusions, TDRD1 overexpression is a common alteration in human prostate cancer which may be exploited for diagnostic or therapeutic procedures.

miR-22 and miR-29a Are Members of the Androgen Receptor Cistrome Modulating LAMC1 and Mcl-1 in Prostate Cancer

The normal prostate as well as early stages and advanced prostate cancer (PCa) require a functional androgen receptor (AR) for growth and survival. The recent discovery of microRNAs (miRNAs) as novel effector molecules of AR disclosed the existence of an intricate network between AR, miRNAs and downstream target genes. In this study DUCaP cells, characterized by high content of wild-type AR and robust AR transcriptional activity, were chosen as the main experimental model. By integrative analysis of chromatin immunoprecipitation-sequencing (ChIP-seq) and microarray expression profiling data, miRNAs putatively bound and significantly regulated by AR were identified. A direct AR regulation of miR-22, miR-29a, and miR-17-92 cluster along with their host genes was confirmed. Interestingly, endogenous levels of miR-22 and miR-29a were found to be reduced in PCa cells expressing AR. In primary tumor samples, miR-22 and miR-29a were less abundant in the cancerous tissue compared with the benign counterpart. This specific expression pattern was associated with a differential DNA methylation of the genomic AR binding sites. The identification of laminin gamma 1 (LAMC1) and myeloid cell leukemia 1 (MCL1) as direct targets of miR-22 and miR-29a, respectively, suggested a tumor-suppressive role of these miRNAs. Indeed, transfection of miRNA mimics in PCa cells induced apoptosis and diminished cell migration and viability. Collectively, these data provide additional information regarding the complex regulatory machinery that guides miRNAs activity in PCa, highlighting an important contribution of miRNAs in the AR signaling.

QSEA – modelling of genome-wide DNA methylation from sequencing enrichment experiments

Abstract Genome-wide enrichment of methylated DNA followed by sequencing (MeDIP-seq) offers a reasonable compromise between experimental costs and genomic coverage. However, the computational analysis of these experiments is complex, and quantification of the enrichment signals in terms of absolute levels of methylation requires specific transformation. In this work, we present QSEA, Quantitative Sequence Enrichment Analysis, a comprehensive workflow for the modelling and subsequent quantification of MeDIP-seq data. As the central part of the workflow we have developed a Bayesian statistical model that transforms the enrichment read counts to absolute levels of methylation and, thus, enhances interpretability and facilitates comparison with other methylation assays. We suggest several calibration strategies for the critical parameters of the model, either using additional data or fairly general assumptions. By comparing the results with bisulfite sequencing (BS) validation data, we show the improvement of QSEA over existing methods. Additionally, we generated a clinically relevant benchmark data set consisting of methylation enrichment experiments (MeDIP-seq), BS-based validation experiments (Methyl-seq) as well as gene expression experiments (RNA-seq) derived from non-small cell lung cancer patients, and show that the workflow retrieves well-known lung tumour methylation markers that are causative for gene expression changes, demonstrating the applicability of QSEA for clinical studies. QSEA is implemented in R and available from the Bioconductor repository 3.4 (www.bioconductor.org/packages/qsea).

Odd skipped-related 1 identifies a population of embryonic fibro-adipogenic progenitors regulating myogenesis during limb development

Fibro-adipogenic progenitors (FAPs) are an interstitial cell population in adult skeletal muscle that support muscle regeneration. During development, interstitial muscle connective tissue (MCT) cells support proper muscle patterning, however the underlying molecular mechanisms are not well understood and it remains unclear whether adult FAPs and embryonic MCT cells share a common lineage. We show here that mouse embryonic limb MCT cells expressing the transcription factor Osr1, differentiate into fibrogenic and adipogenic cells in vivo and in vitro defining an embryonic FAP-like population. Genetic lineage tracing shows that developmental Osr1+ cells give rise to a subset of adult FAPs. Loss of Osr1 function leads to a reduction of myogenic progenitor proliferation and survival resulting in limb muscle patterning defects. Transcriptome and functional analyses reveal that Osr1+ cells provide a critical pro-myogenic niche via the production of MCT specific extracellular matrix components and secreted signaling factors.Fibro-adipogenic progenitors (FAPs) form part of interstitial muscle connective tissue (MCT) in adults but the origin of this non-myogenic lineage is unclear. Here, the authors show that Odd skipped related 1 (Osr1) in mice marks embryonic MCT, giving rise to FAPs, and loss of Osr1 in the limb causes muscle defects.

A dual transcript-discovery approach to improve the delimitation of gene features from RNA-seq data in the chicken model

ABSTRACT The sequence of the chicken genome, like several other draft genome sequences, is presently not fully covered. Gaps, contigs assigned with low confidence and uncharacterized chromosomes result in gene fragmentation and imprecise gene annotation. Transcript abundance estimation from RNA sequencing (RNA-seq) data relies on read quality, library complexity and expression normalization. In addition, the quality of the genome sequence used to map sequencing reads, and the gene annotation that defines gene features, must also be taken into account. A partially covered genome sequence causes the loss of sequencing reads from the mapping step, while an inaccurate definition of gene features induces imprecise read counts from the assignment step. Both steps can significantly bias interpretation of RNA-seq data. Here, we describe a dual transcript-discovery approach combining a genome-guided gene prediction and a de novo transcriptome assembly. This dual approach enabled us to increase the assignment rate of RNA-seq data by nearly 20% as compared to when using only the chicken reference annotation, contributing therefore to a more accurate estimation of transcript abundance. More generally, this strategy could be applied to any organism with partial genome sequence and/or lacking a manually-curated reference annotation in order to improve the accuracy of gene expression studies. Summary: We combined genome-guided gene prediction and whole transcriptome assembly from RNA sequencing data to improve the chicken genome annotation. This method may be also applicable to other imperfectly annotated genomes.

Transcriptomic and proteomic insight into the effects of a defined European mistletoe extract in Ewing sarcoma cells reveals cellular stress responses

BackgroundThe hydrophobic triterpenes, oleanolic and betulinic acid as well as the hydrophilic mistletoe lectins and viscotoxins possess anticancer properties. They do all occur in combination in European mistletoe (Viscum album L.). Commercial Viscum album L. extracts are aqueous, excluding the insoluble triterpenes. We have previously shown that mistletoe lectins and triterpene acids are effective against Ewing sarcoma in vitro, ex vivo and in vivo.MethodsWe recreated a total mistletoe effect (viscumTT) by combining an aqueous extract (viscum) and a triterpene extract (TT) solubilised with cyclodextrins and analysed the effects of viscumTT and the single extracts on TC-71 Ewing sarcoma cells in vitro by transcriptomic and proteomic profiling.ResultsTreatment with the extracts strongly impacted Ewing sarcoma cell gene and protein expression. Apoptosis-associated and stress-activated genes were upregulated, proteasomal protein abundance enhanced and ribosomal and spliceosomal proteins downregulated. The mechanism of action of viscum, TT and viscumTT in TC-71 and MHH-ES-1 cells suggests the involvement of the unfolded protein response. While viscum and viscumTT extract treatment indicate response to oxidative stress and activation of stress-mediated MAPK signalling, TT extract treatment suggests the involvement of TLR signalling and autophagy.ConclusionsSince the combinatory extract viscumTT exerts highly effective pro-apoptotic effects on Ewing sarcoma cells in vitro, this phytopolychemotherapy could be a promising adjuvant therapeutic option for paediatric patients with Ewing sarcoma.