The forkhead transcription factor FOXK2 premarks lineage-specific genes in human embryonic stem cells for activation during differentiation.

Enhancers play important roles in controlling gene expression in a choreographed spatial and temporal manner during development. However, it is unclear how these regulatory regions are established during differentiation. Here we investigated the genome-wide binding profile of the forkhead transcription factor FOXK2 in human embryonic stem cells (ESCs) and downstream cell types. This transcription factor is bound to thousands of regulatory regions in human ESCs, and binding at many sites is maintained as cells differentiate to mesendodermal and neural precursor cell (NPC) types, alongside the emergence of new binding regions. FOXK2 binding is generally associated with active histone marks in any given cell type. Furthermore newly acquired, or retained FOXK2 binding regions show elevated levels of activating histone marks following differentiation to NPCs. In keeping with this association with activating marks, we demonstrate a role for FOXK transcription factors in gene activation during NPC differentiation. FOXK2 occupancy in ESCs is therefore an early mark for delineating the regulatory regions, which become activated in later lineages.

Complexities in the role of acetylation dynamics in modifying inducible gene activation parameters.

High levels of histone acetylation are associated with the regulatory elements of active genes, suggesting a link between acetylation and gene activation. We revisited this model, in the context of EGF-inducible gene expression and found that rather than a simple unifying model, there are two broad classes of genes; one in which high lysine acetylation activity is required for efficient gene activation, and a second group where the opposite occurs and high acetylation activity is inhibitory. We examined the latter class in more detail using EGR2 as a model gene and found that lysine acetylation levels are critical for several activation parameters, including the timing of expression onset, and overall amplitudes of the transcriptional response. In contrast, DUSP1 responds in the canonical manner and its transcriptional activity is promoted by acetylation. Single cell approaches demonstrate heterogenous activation kinetics of a given gene in response to EGF stimulation. Acetylation levels modify these heterogenous patterns and influence both allele activation frequencies and overall expression profile parameters. Our data therefore point to a complex interplay between acetylation equilibria and target gene induction where acetylation level thresholds are an important determinant of transcriptional induction dynamics that are sensed in a gene-specific manner.

Open chromatin profiling identifies AP1 as a transcriptional regulator in oesophageal adenocarcinoma.

Oesophageal adenocarcinoma (OAC) is one of the ten most prevalent forms of cancer and is showing a rapid increase in incidence and yet exhibits poor survival rates. Compared to many other common cancers, the molecular changes that occur in this disease are relatively poorly understood. However, genes encoding chromatin remodeling enzymes are frequently mutated in OAC. This is consistent with the emerging concept that cancer cells exhibit reprogramming of their chromatin environment which leads to subsequent changes in their transcriptional profile. Here, we have used ATAC-seq to interrogate the chromatin changes that occur in OAC using both cell lines and patient-derived material. We demonstrate that there are substantial changes in the regulatory chromatin environment in the cancer cells and using this data we have uncovered an important role for ETS and AP1 transcription factors in driving the changes in gene expression found in OAC cells.

Jun-Mediated Changes in Cell Adhesion Contribute to Mouse Embryonic Stem Cell Exit from Ground State Pluripotency.

Embryonic stem cells (ESC) are able to give rise to any somatic cell type. A lot is known about how ESC pluripotency is maintained, but comparatively less is known about how differentiation is promoted. Cell fate decisions are regulated by interactions between signaling and transcriptional networks. Recent studies have shown that the overexpression or downregulation of the transcription factor Jun can affect the ESC fate. Here we have focussed on the role of the Jun in the exit of mouse ESCs from ground state pluripotency and the onset of early differentiation. Transcriptomic analysis of differentiating ESCs reveals that Jun is required to upregulate a programme of genes associated with cell adhesion as ESCs exit the pluripotent ground state. Several of these Jun-regulated genes are shown to be required for efficient adhesion. Importantly this adhesion is required for the timely regulated exit of ESCs from ground state pluripotency and the onset of early differentiation events. Stem Cells 2016;34:1213-1224.

RUNX1 positively regulates a cell adhesion and migration program in murine hemogenic endothelium prior to blood emergence.

During ontogeny, the transcription factor RUNX1 governs the emergence of definitive hematopoietic cells from specialized endothelial cells called hemogenic endothelium (HE). The ultimate consequence of this endothelial-to-hematopoietic transition is the concomitant activation of the hematopoietic program and downregulation of the endothelial program. However, due to the rare and transient nature of the HE, little is known about the initial role of RUNX1 within this population. We, therefore, developed and implemented a highly sensitive DNA adenine methyltransferase identification-based methodology, including a novel data analysis pipeline, to map early RUNX1 transcriptional targets in HE cells. This novel transcription factor binding site identification protocol should be widely applicable to other low abundance cell types and factors. Integration of the RUNX1 binding profile with gene expression data revealed an unexpected early role for RUNX1 as a positive regulator of cell adhesion- and migration-associated genes within the HE. This suggests that RUNX1 orchestrates HE cell positioning and integration prior to the release of hematopoietic cells. Overall, our genome-wide analysis of the RUNX1 binding and transcriptional profile in the HE provides a novel comprehensive resource of target genes that will facilitate the precise dissection of the role of RUNX1 in early blood development.

Targeted genetic dependency screen facilitates identification of actionable mutations in FGFR4, MAP3K9, and PAK5 in lung cancer.

Approximately 70% of patients with non-small-cell lung cancer present with late-stage disease and have limited treatment options, so there is a pressing need to develop efficacious targeted therapies for these patients. This remains a major challenge as the underlying genetic causes of ~50% of non-small-cell lung cancers remain unknown. Here we demonstrate that a targeted genetic dependency screen is an efficient approach to identify somatic cancer alterations that are functionally important. By using this approach, we have identified three kinases with gain-of-function mutations in lung cancer, namely FGFR4, MAP3K9, and PAK5. Mutations in these kinases are activating toward the ERK pathway, and targeted depletion of the mutated kinases inhibits proliferation, suppresses constitutive activation of downstream signaling pathways, and results in specific killing of the lung cancer cells. Genomic profiling of patients with lung cancer is ushering in an era of personalized medicine; however, lack of actionable mutations presents a significant hurdle. Our study indicates that targeted genetic dependency screens will be an effective strategy to elucidate somatic variants that are essential for lung cancer cell viability.

Evaluation and validation of a robust single cell RNA-amplification protocol through transcriptional profiling of enriched lung cancer initiating cells.

BACKGROUND: Although profiling of RNA in single cells has broadened our understanding of development, cancer biology and mechanisms of disease dissemination, it requires the development of reliable and flexible methods. Here we demonstrate that the EpiStem RNA-Amp™ methodology reproducibly generates microgram amounts of cDNA suitable for RNA-Seq, RT-qPCR arrays and Microarray analysis. RESULTS: Initial experiments compared amplified cDNA generated by three commercial RNA-Amplification protocols (Miltenyi µMACS™ SuperAmp™, NuGEN Ovation® One-Direct System and EpiStem RNA-Amp™) applied to single cell equivalent levels of RNA (25-50 pg) using Affymetrix arrays. The EpiStem RNA-Amp™ kit exhibited the highest sensitivity and was therefore chosen for further testing. A comparison of Affymetrix array data from RNA-Amp™ cDNA generated from single MCF7 and MCF10A cells to reference controls of unamplified cDNA revealed a high degree of concordance. To assess the flexibility of the amplification system single cell RNA-Amp™ cDNA was also analysed using RNA-Seq and high-density qPCR, and showed strong cross-platform correlations. To exemplify the approach we used the system to analyse RNA profiles of small populations of rare cancer initiating cells (CICs) derived from a NSCLC patient-derived xenograft. RNA-Seq analysis was able to identify transcriptional differences in distinct subsets of CIC, with one group potentially enriched for metastasis formation. Pathway analysis revealed that the distinct transcriptional signatures demonstrated in the CIC subpopulations were significantly correlated with published stem-cell and epithelial-mesenchymal transition signatures. CONCLUSIONS: The combined results confirm the sensitivity and flexibility of the RNA-Amp™ method and demonstrate the suitability of the approach for identifying clinically relevant signatures in rare, biologically important cell populations.

The influence of mtDNA deletion on lung cancer cells under the conditions of hypoxia and irradiation.

PURPOSE: This study was to evaluate the influence of mtDNA deletion on the lung cancer cells under the conditions of hypoxia or irradiation. METHOD: The treatment conditions of lung cancer cell lines with (A549) and without mtDNA (ρ0A549: obtained by inducing from A549) included 2 h of hypoxia and 4 Gy irradiation (group 1: without treatment; group 2: 2 h of hypoxia; group 3: 4 Gy irradiation; group 4: 2 h of hypoxia plus 4 Gy irradiation). The Human OneArray™ microarray was used to hybridize with the Cy5-labeled aRNA in microarray sample preparation. Differentially expressed genes (DEGs) between the lung cancer cells with and without mtDNA were identified using NOISeq package in R. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using the online tool of DAVID. RESULT: In the KEGG pathway analysis of down-regulated DEGs, nineteen pathways were simultaneously enriched in the four groups, which were mainly metabolism- and biosynthesis-related pathways. Nine lung cancer-related pathways were enriched in group 4, and more cancer-associated DEGs, such as MYC, MAX, and E2F1 were found in group 4 than in the other groups. CONCLUSION: The mtDNA deletion could inhibit the biosynthesis and metabolism of lung cancer cells and promote the effect of hypoxia and radiation on lung cancer cells. MYC might be the key gene of the cooperation of hypoxia and radiation and MYC, MAX, and E2F1 might play roles in hypoxia- and radiation-induced cell death in lung cancer cells without mtDNA.

Nuclear localization of heparanase 2 (Hpa2) attenuates breast carcinoma growth and metastasis.

Unlike the intense research effort devoted to exploring the significance of heparanase in cancer, very little attention was given to Hpa2, a close homolog of heparanase. Here, we explored the role of Hpa2 in breast cancer. Unexpectedly, we found that patients endowed with high levels of Hpa2 exhibited a higher incidence of tumor metastasis and survived less than patients with low levels of Hpa2. Immunohistochemical examination revealed that in normal breast tissue, Hpa2 localizes primarily in the cell nucleus. In striking contrast, in breast carcinoma, Hpa2 expression is not only decreased but also loses its nuclear localization and appears diffuse in the cell cytoplasm. Importantly, breast cancer patients in which nuclear localization of Hpa2 is retained exhibited reduced lymph-node metastasis, suggesting that nuclear localization of Hpa2 plays a protective role in breast cancer progression. To examine this possibility, we engineered a gene construct that directs Hpa2 to the cell nucleus (Hpa2-Nuc). Notably, overexpression of Hpa2 in breast carcinoma cells resulted in bigger tumors, whereas targeting Hpa2 to the cell nucleus attenuated tumor growth and tumor metastasis. RNAseq analysis was performed to reveal differentially expressed genes (DEG) in Hpa2-Nuc tumors vs. control. The analysis revealed, among others, decreased expression of genes associated with the hallmark of Kras, beta-catenin, and TNF-alpha (via NFkB) signaling. Our results imply that nuclear localization of Hpa2 prominently regulates gene transcription, resulting in attenuation of breast tumorigenesis. Thus, nuclear Hpa2 may be used as a predictive parameter in personalized medicine for breast cancer patients.

ZMYM2 controls human transposable element transcription through distinct co-regulatory complexes.

ZMYM2 is a zinc finger transcriptional regulator that plays a key role in promoting and maintaining cell identity. It has been implicated in several diseases such as congenital anomalies of the kidney where its activity is diminished and cancer where it participates in oncogenic fusion protein events. ZMYM2 is thought to function through promoting transcriptional repression and here we provide more evidence to support this designation. Here we studied ZMYM2 function in human cells and demonstrate that ZMYM2 is part of distinct chromatin-bound complexes including the established LSD1-CoREST-HDAC1 corepressor complex. We also identify new functional and physical interactions with ADNP and TRIM28/KAP1. The ZMYM2-TRIM28 complex forms in a SUMO-dependent manner and is associated with repressive chromatin. ZMYM2 and TRIM28 show strong functional similarity and co-regulate a large number of genes. However, there are no strong links between ZMYM2-TRIM28 binding events and nearby individual gene regulation. Instead, ZMYM2-TRIM28 appears to regulate genes in a more regionally defined manner within TADs where it can directly regulate co-associated retrotransposon expression. We find that different types of ZMYM2 binding complex associate with and regulate distinct subclasses of retrotransposons, with ZMYM2-ADNP complexes at SINEs and ZMYM2-TRIM28 complexes at LTR elements. We propose a model whereby ZMYM2 acts directly through retrotransposon regulation, which may then potentially affect the local chromatin environment and associated coding gene expression.

eRNA profiling uncovers the enhancer landscape of oesophageal adenocarcinoma and reveals new deregulated pathways.

Cancer is driven by both genetic and epigenetic changes that impact on gene expression profiles and the resulting tumourigenic phenotype. Enhancers are transcriptional regulatory elements that are key to our understanding of how this rewiring of gene expression is achieved in cancer cells. Here, we have harnessed the power of RNA-seq data from hundreds of patients with oesophageal adenocarcinoma (OAC) or its precursor state Barrett's oesophagus coupled with open chromatin maps to identify potential enhancer RNAs and their associated enhancer regions in this cancer. We identify ~1000 OAC-specific enhancers and use these data to uncover new cellular pathways that are operational in OAC. Among these are enhancers for JUP, MYBL2, and CCNE1, and we show that their activity is required for cancer cell viability. We also demonstrate the clinical utility of our dataset for identifying disease stage and patient prognosis. Our data therefore identify an important set of regulatory elements that enhance our molecular understanding of OAC and point to potential new therapeutic directions.

Programmed fluctuations in sense/antisense transcript ratios drive sexual differentiation in S. pombe.

Strand-specific RNA sequencing of S. pombe revealed a highly structured programme of ncRNA expression at over 600 loci. Waves of antisense transcription accompanied sexual differentiation. A substantial proportion of ncRNA arose from mechanisms previously considered to be largely artefactual, including improper 3' termination and bidirectional transcription. Constitutive induction of the entire spk1+, spo4+, dis1+ and spo6+ antisense transcripts from an integrated, ectopic, locus disrupted their respective meiotic functions. This ability of antisense transcripts to disrupt gene function when expressed in trans suggests that cis production at native loci during sexual differentiation may also control gene function. Consistently, insertion of a marker gene adjacent to the dis1+ antisense start site mimicked ectopic antisense expression in reducing the levels of this microtubule regulator and abolishing the microtubule-dependent 'horsetail' stage of meiosis. Antisense production had no impact at any of these loci when the RNA interference (RNAi) machinery was removed. Thus, far from being simply 'genome chatter', this extensive ncRNA landscape constitutes a fundamental component in the controls that drive the complex programme of sexual differentiation in S. pombe.

A comparison of massively parallel nucleotide sequencing with oligonucleotide microarrays for global transcription profiling.

BACKGROUND: RNA-Seq exploits the rapid generation of gigabases of sequence data by Massively Parallel Nucleotide Sequencing, allowing for the mapping and digital quantification of whole transcriptomes. Whilst previous comparisons between RNA-Seq and microarrays have been performed at the level of gene expression, in this study we adopt a more fine-grained approach. Using RNA samples from a normal human breast epithelial cell line (MCF-10a) and a breast cancer cell line (MCF-7), we present a comprehensive comparison between RNA-Seq data generated on the Applied Biosystems SOLiD platform and data from Affymetrix Exon 1.0ST arrays. The use of Exon arrays makes it possible to assess the performance of RNA-Seq in two key areas: detection of expression at the granularity of individual exons, and discovery of transcription outside annotated loci. RESULTS: We found a high degree of correspondence between the two platforms in terms of exon-level fold changes and detection. For example, over 80% of exons detected as expressed in RNA-Seq were also detected on the Exon array, and 91% of exons flagged as changing from Absent to Present on at least one platform had fold-changes in the same direction. The greatest detection correspondence was seen when the read count threshold at which to flag exons Absent in the SOLiD data was set to t<1 suggesting that the background error rate is extremely low in RNA-Seq. We also found RNA-Seq more sensitive to detecting differentially expressed exons than the Exon array, reflecting the wider dynamic range achievable on the SOLiD platform. In addition, we find significant evidence of novel protein coding regions outside known exons, 93% of which map to Exon array probesets, and are able to infer the presence of thousands of novel transcripts through the detection of previously unreported exon-exon junctions. CONCLUSIONS: By focusing on exon-level expression, we present the most fine-grained comparison between RNA-Seq and microarrays to date. Overall, our study demonstrates that data from a SOLiD RNA-Seq experiment are sufficient to generate results comparable to those produced from Affymetrix Exon arrays, even using only a single replicate from each platform, and when presented with a large genome.

Cooperative behaviour and phenotype plasticity evolve during melanoma progression.

A major challenge for managing melanoma is its tumour heterogeneity based on individual co-existing melanoma cell phenotypes. These phenotypes display variable responses to standard therapies, and they drive individual steps of melanoma progression; hence, understanding their behaviour is imperative. Melanoma phenotypes are defined by distinct transcriptional states, which relate to different melanocyte lineage development phases, ranging from a mesenchymal, neural crest-like to a proliferative, melanocytic phenotype. It is thought that adaptive phenotype plasticity based on transcriptional reprogramming drives melanoma progression, but at which stage individual phenotypes dominate and moreover, how they interact is poorly understood. We monitored melanocytic and mesenchymal phenotypes throughout melanoma progression and detected transcriptional reprogramming at different stages, with a gain in mesenchymal traits in circulating melanoma cells (CTCs) and proliferative features in metastatic tumours. Intriguingly, we found that distinct phenotype populations interact in a cooperative manner, which generates tumours of greater "fitness," supports CTCs and expands organotropic cues in metastases. Fibronectin, expressed in mesenchymal cells, acts as key player in cooperativity and promotes survival of melanocytic cells. Our data reveal an important role for inter-phenotype communications at various stages of disease progression, suggesting these communications could act as therapeutic target.

Geno2proteo, a Tool for Batch Retrieval of DNA and Protein Sequences from Any Genomic or Protein Regions.

The interconversion of sequences that constitute the genome and the proteome is becoming increasingly important due to the generation of large amounts of DNA sequence data. Following mapping of DNA segments to the genome, one fundamentally important task is to find the amino acid sequences which are coded within a list of genomic sections. Conversely, given a series of protein segments, an important task is to find the genomic loci which code for a list of protein regions. To perform these tasks on a region by region basis is extremely laborious when a large number of regions are being studied. We have therefore implemented an R package geno2proteo which performs the two mapping tasks and subsequent sequence retrieval in a batch fashion. In order to make the tool more accessible to users, we have created a web interface of the R package which allows the users to perform the mapping tasks by going to the web page http://sharrocksresources.manchester.ac.uk/tofigaps and using the web service.

The influence of tumor heterogeneity on sensitivity of EGFR-mutant lung adenocarcinoma cells to EGFR-TKIs.

BACKGROUND: Studies have shown that extensive genetic or spatial heterogeneity is present within tumors. The present study explored the influence of heterogeneous cells in the tumor microenvironment (TME) on the sensitivity of EGFR-mutant lung adenocarcinoma cells to EGFR-TKIs and further investigated associated molecular mechanisms. METHODS: Tumor heterogeneity was simulated using transwell co-culture technique with H1975, A549 or MRC-5 cells grown in the upper chambers and PC-9 cells cultured in the bottom chamber. Each co-culture system was grouped based on different proportions of cells in the upper and lower chambers (from 1% to 300%), and each group was subdivided into erlotinib treatment group (erlotinib+) and erlotinib non-treatment group (erlotinib-). The viability of PC-9 cells was analyzed by CCK-8; HGF, IGF-1 and TGF-α were determined by ELISA; MET amplification was detected by qRT-PCR, and the relevant proteins were detected by Western blot. RESULTS: The viability of PC-9 cells increased with increased proportion of A549/PC-9 or MRC-5/PC-9 cells from 1% to 300%. HGF increased with increased proportion of H1975 cells from 1% to 300%. In all three co-culture systems, TGF-α production in the erlotinib+ subgroup was significantly lower than that in erlotinib- subgroup, and phosphorylated AKT protein showed an ascending tendency with the increased proportion of upper chamber cells relative to PC-9 cells from 1% to 300%. H1975 cells could induce MET amplification of PC-9 cells. MRC5 cells inhibited MET amplification. CONCLUSIONS: Tumor heterogeneity partially accounts for the resistance of EGFR-mutant lung adenocarcinoma cells to EGFR-TKIs. The possible mechanisms may be related to AKT signaling pathways activated by growth factors, which are secreted by heterogeneous cells in the TME under the pressure of EGFR-TKIs.

Enhancer Activation by Pharmacologic Displacement of LSD1 from GFI1 Induces Differentiation in Acute Myeloid Leukemia.

Pharmacologic inhibition of LSD1 promotes blast cell differentiation in acute myeloid leukemia (AML) with MLL translocations. The assumption has been that differentiation is induced through blockade of LSD1's histone demethylase activity. However, we observed that rapid, extensive, drug-induced changes in transcription occurred without genome-wide accumulation of the histone modifications targeted for demethylation by LSD1 at sites of LSD1 binding and that a demethylase-defective mutant rescued LSD1 knockdown AML cells as efficiently as wild-type protein. Rather, LSD1 inhibitors disrupt the interaction of LSD1 and RCOR1 with the SNAG-domain transcription repressor GFI1, which is bound to a discrete set of enhancers located close to transcription factor genes that regulate myeloid differentiation. Physical separation of LSD1/RCOR1 from GFI1 is required for drug-induced differentiation. The consequent inactivation of GFI1 leads to increased enhancer histone acetylation within hours, which directly correlates with the upregulation of nearby subordinate genes.

EINCR1 is an EGF inducible lincRNA overexpressed in lung adenocarcinomas.

Long non-coding RNAs are being increasingly recognised as important molecules involved in regulating a diverse array of biological functions. For example, many long non-coding RNAs have been associated with tumourigenesis and in this context their molecular functions often involves impacting on chromatin and transcriptional control processes. One important cellular control system that is often deregulated in cancer cells is the ERK MAP kinase pathway. Here we have investigated whether ERK pathway signaling in response to EGF stimulation, leads to changes in the production of long non-coding RNAs. We identify several different classes of EGF pathway-regulated lncRNAs. We focus on one of the inducible lincRNAs, EGF inducible long intergenic non-coding RNA 1 (EINCR1). EINCR1 is predominantly nuclear and shows delayed activation kinetics compared to other immediate-early EGF-inducible genes. In humans it is expressed in a tissue-specific manner and is mainly confined to the heart but it exhibits little evolutionary conservation. Importantly, in several cancers EINCR1 shows elevated expression levels which correlate with poor survival in lung adenocarcinoma patients. In the context of lung adenocarcinomas, EINCR1 expression is anti-correlated with the expression of several protein coding EGF-regulated genes. A potential functional connection is demonstrated as EINCR1 overexpression is shown to reduce the expression of EGF-regulated protein coding genes including FOS and FOSB.

Molecular analysis of circulating tumor cells identifies distinct copy-number profiles in patients with chemosensitive and chemorefractory small-cell lung cancer.

In most patients with small-cell lung cancer (SCLC)-a metastatic, aggressive disease-the condition is initially chemosensitive but then relapses with acquired chemoresistance. In a minority of patients, however, relapse occurs within 3 months of initial treatment; in these cases, disease is defined as chemorefractory. The molecular mechanisms that differentiate chemosensitive from chemorefractory disease are currently unknown. To identify genetic features that distinguish chemosensitive from chemorefractory disease, we examined copy-number aberrations (CNAs) in circulating tumor cells (CTCs) from pretreatment SCLC blood samples. After analysis of 88 CTCs isolated from 13 patients (training set), we generated a CNA-based classifier that we validated in 18 additional patients (testing set, 112 CTC samples) and in six SCLC patient-derived CTC explant tumors. The classifier correctly assigned 83.3% of the cases as chemorefractory or chemosensitive. Furthermore, a significant difference was observed in progression-free survival (PFS) (Kaplan-Meier P value = 0.0166) between patients designated as chemorefractory or chemosensitive by using the baseline CNA classifier. Notably, CTC CNA profiles obtained at relapse from five patients with initially chemosensitive disease did not switch to a chemorefractory CNA profile, which suggests that the genetic basis for initial chemoresistance differs from that underlying acquired chemoresistance.