Complex exon-intron marking by histone modifications is not determined solely by nucleosome distribution

It has recently been shown that nucleosome distribution, histone modifications and RNA polymerase II (Pol II) occupancy show preferential association with exons (“exon-intron marking”), linking chromatin structure and function to co-transcriptional splicing in a variety of eukaryotes. Previous ChIP-sequencing studies suggested that these marking patterns reflect the nucleosomal landscape. By analyzing ChIP-chip datasets across the human genome in three cell types, we have found that this marking system is far more complex than previously observed. We show here that a range of histone modifications and Pol II are preferentially associated with exons. However, there is noticeable cell-type specificity in the degree of exon marking by histone modifications and, surprisingly, this is also reflected in some histone modifications patterns showing biases towards introns. Exon-intron marking is laid down in the absence of transcription on silent genes, with some marking biases changing or becoming reversed for genes expressed at different levels. Furthermore, the relationship of this marking system with splicing is not simple, with only some histone modifications reflecting exon usage/inclusion, while others mirror patterns of exon exclusion. By examining nucleosomal distributions in all three cell types, we demonstrate that these histone modification patterns cannot solely be accounted for by differences in nucleosome levels between exons and introns. In addition, because of inherent differences between ChIP-chip array and ChIP-sequencing approaches, these platforms report different nucleosome distribution patterns across the human genome. Our findings confound existing views and point to active cellular mechanisms which dynamically regulate histone modification levels and account for exon-intron marking. We believe that these histone modification patterns provide links between chromatin accessibility, Pol II movement and co-transcriptional splicing

FusionFinder: A Software Tool to Identify Expressed Gene Fusion Candidates from RNA-Seq Data

The hallmarks of many haematological malignancies and solid tumours are chromosomal translocations, which may lead to gene fusions. Recently, next-generation sequencing techniques at the transcriptome level (RNA-Seq) have been used to verify known and discover novel transcribed gene fusions. We present FusionFinder, a Perl-based software designed to automate the discovery of candidate gene fusion partners from single-end (SE) or paired-end (PE) RNA-Seq read data. FusionFinder was applied to data from a previously published analysis of the K562 chronic myeloid leukaemia (CML) cell line. Using FusionFinder we successfully replicated the findings of this study and detected additional previously unreported fusion genes in their dataset, which were confirmed experimentally. These included two isoforms of a fusion involving the genes BRK1 and VHL, whose co-deletion has previously been associated with the prevalence and severity of renal-cell carcinoma. FusionFinder is made freely available for non-commercial use and can be downloaded from the project website (http://bioinformatics.childhealthresearch.org.au/software/fusionfinder/)

Natural Killer Cell Mediated Cytotoxic Responses in the Tasmanian Devil

The Tasmanian devil (Sarcophilus harrisii), the world's largest marsupial carnivore, is under threat of extinction following the emergence of an infectious cancer. Devil facial tumour disease (DFTD) is spread between Tasmanian devils during biting. The disease is consistently fatal and devils succumb without developing a protective immune response. The aim of this study was to determine if Tasmanian devils were capable of forming cytotoxic antitumour responses and develop antibodies against DFTD cells and foreign tumour cells. The two Tasmanian devils immunised with irradiated DFTD cells did not form cytotoxic or humoral responses against DFTD cells, even after multiple immunisations. However, following immunisation with xenogenic K562 cells, devils did produce cytotoxic responses and antibodies against this foreign tumour cell line. The cytotoxicity appeared to occur through the activity of natural killer (NK) cells in an antibody dependent manner. Classical NK cell responses, such as innate killing of DFTD and foreign cancer cells, were not observed. Cells with an NK-like phenotype comprised approximately 4 percent of peripheral blood mononuclear cells. The results of this study suggest that Tasmanian devils have NK cells with functional cytotoxic pathways. Although devil NK cells do not directly recognise DFTD cancer cells, the development of antibody dependent cell-mediated cytotoxicity presents a potential pathway to induce cytotoxic responses against the disease. These findings have positive implications for future DFTD vaccine research

Knockdown of ZNF268, which Is Transcriptionally Downregulated by GATA-1, Promotes Proliferation of K562 Cells

The human ZNF268 gene encodes a typical KRAB-C2H2 zinc finger protein that may participate in hematopoiesis and leukemogenesis. A recent microarray study revealed that ZNF268 expression continuously decreases during erythropoiesis. However, the molecular mechanisms underlying regulation of ZNF268 during hematopoiesis are not well understood. Here we found that GATA-1, a master regulator of erythropoiesis, repressed the promoter activity and transcription of ZNF268. Electrophoretic mobility shift assays and chromatin immunoprecipitation assays showed that GATA-1 directly bound to a GATA binding site in the ZNF268 promoter in vitro and in vivo. Knockdown of ZNF268 in K562 erythroleukemia cells with specific siRNA accelerated cellular proliferation, suppressed apoptosis, and reduced expression of erythroid-specific developmental markers. It also promoted growth of subcutaneous K562-derived tumors in nude mice. These results suggest that ZNF268 is a crucial downstream target and effector of GATA-1. They also suggest the downregulation of ZNF268 by GATA-1 is important in promoting the growth and suppressing the differentiation of K562 erythroleukemia cells

TMEM8 – a non-globin gene entrapped in the globin web

For more than 30 years it was believed that globin gene domains included only genes encoding globin chains. Here we show that in chickens, the domain of α-globin genes also harbor the non-globin gene TMEM8. It was relocated to the vicinity of the α-globin cluster due to inversion of an ∼170-kb genomic fragment. Although in humans TMEM8 is preferentially expressed in resting T-lymphocytes, in chickens it acquired an erythroid-specific expression profile and is upregulated upon terminal differentiation of erythroblasts. This correlates with the presence of erythroid-specific regulatory elements in the body of chicken TMEM8, which interact with regulatory elements of the α-globin genes. Surprisingly, TMEM8 is not simply recruited to the α-globin gene domain active chromatin hub. An alternative chromatin hub is assembled, which includes some of the regulatory elements essential for the activation of globin gene expression. These regulatory elements should thus shuttle between two different chromatin hubs

Involvement of cell surface TG2 in the aggregation of K562 cells triggered by gluten

Gluten-induced aggregation of K562 cells represents an in vitro model reproducing the early steps occurring in the small bowel of celiac patients exposed to gliadin. Despite the clear involvement of TG2 in the activation of the antigen-presenting cells, it is not yet clear in which compartment it occurs. Herein we study the calcium-dependent aggregation of these cells, using either cell-permeable or cell-impermeable TG2 inhibitors. Gluten induces efficient aggregation when calcium is absent in the extracellular environment, while TG2 inhibitors do not restore the full aggregating potential of gluten in the presence of calcium. These findings suggest that TG2 activity is not essential in the cellular aggregation mechanism. We demonstrate that gluten contacts the cells and provokes their aggregation through a mechanism involving the A-gliadin peptide 31-43. This peptide also activates the cell surface associated extracellular TG2 in the absence of calcium. Using a bioinformatics approach, we identify the possible docking sites of this peptide on the open and closed TG2 structures. Peptide docks with the closed TG2 structure near to the GTP/GDP site, by establishing molecular interactions with the same amino acids involved in stabilization of GTP binding. We suggest that it may occur through the displacement of GTP, switching the TG2 structure from the closed to the active open conformation. Furthermore, docking analysis shows peptide binding with the β-sandwich domain of the closed TG2 structure, suggesting that this region could be responsible for the different aggregating effects of gluten shown in the presence or absence of calcium. We deduce from these data a possible mechanism of action by which gluten makes contact with the cell surface, which could have possible implications in the celiac disease onset

The Recognition of N-Glycans by the Lectin ArtinM Mediates Cell Death of a Human Myeloid Leukemia Cell Line

ArtinM, a d-mannose-binding lectin from Artocarpus heterophyllus (jackfruit), interacts with N-glycosylated receptors on the surface of several cells of hematopoietic origin, triggering cell migration, degranulation, and cytokine release. Because malignant transformation is often associated with altered expression of cell surface glycans, we evaluated the interaction of ArtinM with human myelocytic leukemia cells and investigated cellular responses to lectin binding. The intensity of ArtinM binding varied across 3 leukemia cell lines: NB4>K562>U937. The binding, which was directly related to cell growth suppression, was inhibited in the presence of Manα1-3(Manα1-6)Manβ1, and was reverted in underglycosylated NB4 cells. ArtinM interaction with NB4 cells induced cell death (IC50 = 10 µg/mL), as indicated by cell surface exposure of phosphatidylserine and disruption of mitochondrial membrane potential unassociated with caspase activation or DNA fragmentation. Moreover, ArtinM treatment of NB4 cells strongly induced reactive oxygen species generation and autophagy, as indicated by the detection of acidic vesicular organelles in the treated cells. NB4 cell death was attributed to ArtinM recognition of the trimannosyl core of N-glycans containing a ß1,6-GlcNAc branch linked to α1,6-mannose. This modification correlated with higher levels of N-acetylglucosaminyltransferase V transcripts in NB4 cells than in K562 or U937 cells. Our results provide new insights into the potential of N-glycans containing a β1,6-GlcNAc branch linked to α1,6-mannose as a novel target for anti-leukemia treatment

Prediction of RNA Polymerase II recruitment, elongation and stalling from histone modification data

<p>Abstract</p> <p>Background</p> <p>Initiation and elongation of RNA polymerase II (RNAPII) transcription is regulated by both DNA sequence and chromatin signals. Recent breakthroughs make it possible to measure the chromatin state and activity of core promoters genome-wide, but dedicated computational strategies are needed to progress from descriptive annotation of data to quantitative, predictive models.</p> <p>Results</p> <p>Here, we describe a computational framework which with high accuracy can predict the locations of core promoters, the amount of recruited RNAPII at the promoter, the amount of elongating RNAPII in the gene body, the mRNA production originating from the promoter and finally also the stalling characteristics of RNAPII by considering both quantitative and spatial features of histone modifications around the transcription start site (TSS).</p> <p>As the model framework can also pinpoint the signals that are the most influential for prediction, it can be used to infer underlying regulatory biology. For example, we show that the H3K4 di- and tri- methylation signals are strongly predictive for promoter location while the acetylation marks H3K9 and H3K27 are highly important in estimating the promoter usage. All of these four marks are found to be necessary for recruitment of RNAPII but not sufficient for the elongation. We also show that the spatial distributions of histone marks are almost as predictive as the signal strength and that a set of histone marks immediately downstream of the TSS is highly predictive of RNAPII stalling.</p> <p>Conclusions</p> <p>In this study we introduce a general framework to accurately predict the level of RNAPII recruitment, elongation, stalling and mRNA expression from chromatin signals. The versatility of the method also makes it ideally suited to investigate other genomic data.</p

Small Molecule Amiloride Modulates Oncogenic RNA Alternative Splicing to Devitalize Human Cancer Cells

Alternative splicing involves differential exon selection of a gene transcript to generate mRNA and protein isoforms with structural and functional diversity. Abnormal alternative splicing has been shown to be associated with malignant phenotypes of cancer cells, such as chemo-resistance and invasive activity. Screening small molecules and drugs for modulating RNA splicing in human hepatocellular carcinoma cell line Huh-7, we discovered that amiloride, distinct from four pH-affecting amiloride analogues, could “normalize” the splicing of BCL-X, HIPK3 and RON/MISTR1 transcripts. Our proteomic analyses of amiloride-treated cells detected hypo-phosphorylation of splicing factor SF2/ASF, and decreased levels of SRp20 and two un-identified SR proteins. We further observed decreased phosphorylation of AKT, ERK1/2 and PP1, and increased phosphorylation of p38 and JNK, suggesting that amiloride treatment down-regulates kinases and up-regulates phosphatases in the signal pathways known to affect splicing factor protein phosphorylation. These amiloride effects of “normalized” oncogenic RNA splicing and splicing factor hypo-phosphorylation were both abrogated by pre-treatment with a PP1 inhibitor. Global exon array of amiloride-treated Huh-7 cells detected splicing pattern changes involving 584 exons in 551 gene transcripts, many of which encode proteins playing key roles in ion transport, cellular matrix formation, cytoskeleton remodeling, and genome maintenance. Cellular functional analyses revealed subsequent invasion and migration defects, cell cycle disruption, cytokinesis impairment, and lethal DNA degradation in amiloride-treated Huh-7 cells. Other human solid tumor and leukemic cells, but not a few normal cells, showed similar amiloride-altered RNA splicing with devitalized consequence. This study thus provides mechanistic underpinnings for exploiting small molecule modulation of RNA splicing for cancer therapeutics