Detecting and visualizing gene fusions.

In recent years, gene fusions have gained significant recognition as biomarkers. They can assist treatment decisions, are seldom found in normal tissue and are detectable through Next-generation sequencing (NGS) of the transcriptome (RNA-seq). To transform the data provided by the sequencer into robust gene fusion detection several analysis steps are needed. Usually the first step is to map the sequenced transcript fragments (RNA-seq) to a reference genome. One standard application of this approach is to estimate expression and detect variants within known genes, e.g. SNPs and indels. In case of gene fusions, however, completely novel gene structures have to be detected. Here, we describe the detection of such gene fusion events based on our comprehensive transcript annotation (ElDorado). To demonstrate the utility of our approach, we extract gene fusion candidates from eight breast cancer cell lines, which we compare to experimentally verified gene fusions. We discuss several gene fusion events, like BCAS3-BCAS4 that was only detected in the breast cancer cell line MCF7. As supporting evidence we show that gene fusions occur more frequently in copy number enriched regions (CNV analysis). In addition, we present the Transcriptome Viewer (TViewer) a tool that allows to interactively visualize gene fusions. Finally, we support detected gene fusions through literature mining based annotations and network analyses. In conclusion, we present a platform that allows detecting gene fusions and supporting them through literature knowledge as well as rich visualization capabilities. This enables scientists to better understand molecular processes, biological functions and disease associations, which will ultimately lead to better biomedical knowledge for the development of biomarkers for diagnostics and therapies.

RNA-sequencing as useful screening tool in the combat against the misuse of anabolic agents.

The abuse of anabolic substances in animal husbandry is forbidden within the EU and well controlled by detecting substance residues in different matrices. The application of newly designed drugs or substance cocktails represents big problems. Therefore developing sensitive test methods is important. The analysis of physiological changes caused by the use of anabolic agents on the molecular level, for example, by quantifying gene expression response, is a new approach to develop such screening methods. A novel technology for holistic gene expression analysis is RNA sequencing. In this study, the potential of this high-throughput method for the identification of biomarkers was evaluated. The effect of trenbolone acetate plus estradiol on gene expression in liver from Nguni heifers was analyzed with RNA sequencing. The expression of 40 selected candidate genes was verified via RT-qPCR, whereby 20 of these genes were significantly regulated. To extract the intended information from these regulated genes, biostatistical tools for pattern recognition were applied and resulted in a clear separation of the treatment groups. Those candidate genes could be verified in boars and in calves treated with anabolic substances. These results show the potential of RNA sequencing to screen for biomarker candidates to detect the abuse of anabolics. The verification of these biomarkers in boars and calves leads to the assumption that gene expression biomarkers are independent of breed or even species and that biomarkers, identified in farm animals could also act as potential biomarker candidates to detect the abuse of anabolic substances in human sports.

Estrogen receptor alpha controls a gene network in luminal-like breast cancer cells comprising multiple transcription factors and microRNAs.

Luminal-like breast tumor cells express estrogen receptor alpha (ERalpha), a member of the nuclear receptor family of ligand-activated transcription factors that controls their proliferation, survival, and functional status. To identify the molecular determinants of this hormone-responsive tumor phenotype, a comprehensive genome-wide analysis was performed in estrogen stimulated MCF-7 and ZR-75.1 cells by integrating time-course mRNA expression profiling with global mapping of genomic ERalpha binding sites by chromatin immunoprecipitation coupled to massively parallel sequencing, microRNA expression profiling, and in silico analysis of transcription units and receptor binding regions identified. All 1270 genes that were found to respond to 17beta-estradiol in both cell lines cluster in 33 highly concordant groups, each of which showed defined kinetics of RNA changes. This hormone-responsive gene set includes several direct targets of ERalpha and is organized in a gene regulation cascade, stemming from ligand-activated receptor and reaching a large number of downstream targets via AP-2gamma, B-cell activating transcription factor, E2F1 and 2, E74-like factor 3, GTF2IRD1, hairy and enhancer of split homologue-1, MYB, SMAD3, RARalpha, and RXRalpha transcription factors. MicroRNAs are also integral components of this gene regulation network because miR-107, miR-424, miR-570, miR-618, and miR-760 are regulated by 17beta-estradiol along with other microRNAs that can target a significant number of transcripts belonging to one or more estrogen-responsive gene clusters.

RNA-sequencing reveals that STRN, ZNF484 and WNK1 add to the value of mitochondrial MT-COI and COX10 as markers of unstable coronary artery disease.

Markers in monocytes, precursors of macrophages, which are related to CAD, are largely unknown. Therefore, we aimed to identify genes in monocytes predictive of a new ischemic event in patients with CAD and/or discriminate between stable CAD and acute coronary syndrome. We included 66 patients with stable CAD, of which 24 developed a new ischemic event, and 19 patients with ACS. Circulating CD14+ monocytes were isolated with magnetic beads. RNA sequencing analysis in monocytes of patients with (n = 13) versus without (n = 11) ischemic event at follow-up and in patients with ACS (n = 12) was validated with qPCR (n = 85). MT-COI, STRN and COX10 predicted new ischemic events in CAD patients (power for separation at 1% error rate of 0.97, 0.90 and 0.77 respectively). Low MT-COI and high STRN were also related to shorter time between blood sampling and event. COX10 and ZNF484 together with MT-COI, STRN and WNK1 separated ACS completely from stable CAD patients. RNA expressions in monocytes of MT-COI, COX10, STRN, WNK1 and ZNF484 were independent of cholesterol lowering and antiplatelet treatment. They were independent of troponin T, a marker of myocardial injury. But, COX10 and ZNF484 in human plaques correlated to plaque markers of M1 macrophage polarization, reflecting vascular injury. Expression of MT-COI, COX10, STRN and WNK1, but not that of ZNF484, PBMCs paired with that in monocytes. The prospective study of relation of MT-COI, COX10, STRN, WNK1 and ZNF484 with unstable CAD is warranted.

Transcriptional profiling of summer wheat, grown under different realistic UV-B irradiation regimes.

There is limited information on the impact of present-day ultraviolet-B (UV-B) radiation on a reprogramming of gene expression in crops. Summer wheat was cultivated in controlled environmental facilities under simulated realistic climatic conditions. We investigated the effect of different regimes of UV-B radiation on summer wheat (Triticum aestivum L.) cultivars Nandu, Star and Turbo. Until recently, these were most important in Bavaria. Different cultivars of crops often show great differences in their sensitivity towards UV-B radiation. To identify genes that might be involved in UV-B defence mechanisms, we first analyzed selected genes known to be involved in plant defence mechanisms. RNA gel blot analysis of RNA isolated from the flag leaf of 84-day-old plants showed differences in transcript levels among the cultivars. Flag leaves are known to be important for grain development, which was completed at 84 days post-anthesis. Catalase 2 (Cat2) transcripts were elevated by increased UV irradiation in all cultivars with highest levels in cv. Nandu. Pathogenesis-related protein 1 (PR1) transcripts were elevated only in cv. Star. A minor influence on transcripts for phenylalanine ammonia-lyase (PAL) was observed in all three cultivars. This indicates different levels of acclimation to UV-B radiation in the wheat cultivars studied. To analyze these responses in more detail, UV-B-exposed flag leaves of 84-day-old wheat (cv. Nandu) were pooled to isolate cDNAs of induced genes by suppression-subtractive hybridization (SSH). Among the initially isolated cDNA clones, 13 were verified by RNA gel blot analysis showing an up-regulation at elevated levels of UV-B radiation. Functional classification revealed genes encoding proteins associated with protein assembly, chaperonins, programmed cell death and signal transduction. We also studied growth, flowering time, ear development and yield as more typical agricultural parameters. Plant growth of young plants was reduced at increased UV-B radiation. Flowering and ear development were delayed concomitantly, whereas total grain weight was not influenced at any of the UV-B irradiation regimes.

Oestrogen receptor ß regulates epigenetic patterns at specific genomic loci through interaction with thymine DNA glycosylase.

BACKGROUND: DNA methylation is one way to encode epigenetic information and plays a crucial role in regulating gene expression during embryonic development. DNA methylation marks are established by the DNA methyltransferases and, recently, a mechanism for active DNA demethylation has emerged involving the ten-eleven translocator proteins and thymine DNA glycosylase (TDG). However, so far it is not clear how these enzymes are recruited to, and regulate DNA methylation at, specific genomic loci. A number of studies imply that sequence-specific transcription factors are involved in targeting DNA methylation and demethylation processes. Oestrogen receptor beta (ERß) is a ligand-inducible transcription factor regulating gene expression in response to the female sex hormone oestrogen. Previously, we found that ERß deficiency results in changes in DNA methylation patterns at two gene promoters, implicating an involvement of ERß in DNA methylation. In this study, we set out to explore this involvement on a genome-wide level, and to investigate the underlying mechanisms of this function. RESULTS: Using reduced representation bisulfite sequencing, we compared genome-wide DNA methylation in mouse embryonic fibroblasts derived from wildtype and ERß knock-out mice, and identified around 8000 differentially methylated positions (DMPs). Validation and further characterisation of selected DMPs showed that differences in methylation correlated with changes in expression of the nearest gene. Additionally, re-introduction of ERß into the knock-out cells could reverse hypermethylation and reactivate expression of some of the genes. We also show that ERß is recruited to regions around hypermethylated DMPs. Finally, we demonstrate here that ERß interacts with TDG and that TDG binds ERß-dependently to hypermethylated DMPs. CONCLUSION: We provide evidence that ERß plays a role in regulating DNA methylation at specific genomic loci, likely as the result of its interaction with TDG at these regions. Our findings imply a novel function of ERß, beyond direct transcriptional control, in regulating DNA methylation at target genes. Further, they shed light on the question how DNA methylation is regulated at specific genomic loci by supporting a concept in which sequence-specific transcription factors can target factors that regulate DNA methylation patterns.

Deep sequencing of MYC DNA-binding sites in Burkitt lymphoma.

BACKGROUND: MYC is a key transcription factor involved in central cellular processes such as regulation of the cell cycle, histone acetylation and ribosomal biogenesis. It is overexpressed in the majority of human tumors including aggressive B-cell lymphoma. Especially Burkitt lymphoma (BL) is a highlight example for MYC overexpression due to a chromosomal translocation involving the c-MYC gene. However, no genome-wide analysis of MYC-binding sites by chromatin immunoprecipitation (ChIP) followed by next generation sequencing (ChIP-Seq) has been conducted in BL so far. METHODOLOGY/PRINCIPAL FINDINGS: ChIP-Seq was performed on 5 BL cell lines with a MYC-specific antibody giving rise to 7,054 MYC-binding sites after bioinformatics analysis of a total of approx. 19 million sequence reads. In line with previous findings, binding sites accumulate in gene sets known to be involved in the cell cycle, ribosomal biogenesis, histone acetyltransferase and methyltransferase complexes demonstrating a regulatory role of MYC in these processes. Unexpectedly, MYC-binding sites also accumulate in many B-cell relevant genes. To assess the functional consequences of MYC binding, the ChIP-Seq data were supplemented with siRNA- mediated knock-downs of MYC in BL cell lines followed by gene expression profiling. Interestingly, amongst others, genes involved in the B-cell function were up-regulated in response to MYC silencing. CONCLUSION/SIGNIFICANCE: The 7,054 MYC-binding sites identified by our ChIP-Seq approach greatly extend the knowledge regarding MYC binding in BL and shed further light on the enormous complexity of the MYC regulatory network. Especially our observations that (i) many B-cell relevant genes are targeted by MYC and (ii) that MYC down-regulation leads to an up-regulation of B-cell genes highlight an interesting aspect of BL biology.