Dan E. Webster

Control of somatic tissue differentiation by the long non-coding RNA TINCR

Several of the thousands of human long non-coding RNAs (lncRNAs) have been functionally characterized; however, potential roles for lncRNAs in somatic tissue differentiation remain poorly understood. Here we show that a 3.7-kilobase lncRNA, terminal differentiation-induced ncRNA (TINCR), controls human epidermal differentiation by a post-transcriptional mechanism. TINCR is required for high messenger RNA abundance of key differentiation genes, many of which are mutated in human skin diseases, including FLG, LOR, ALOXE3, ALOX12B, ABCA12, CASP14 and ELOVL3. TINCR-deficient epidermis lacked terminal differentiation ultrastructure, including keratohyalin granules and intact lamellar bodies. Genome-scale RNA interactome analysis revealed that TINCR interacts with a range of differentiation mRNAs. TINCR–mRNA interaction occurs through a 25-nucleotide ‘TINCR box’ motif that is strongly enriched in interacting mRNAs and required for TINCR binding. A high-throughput screen to analyse TINCR binding capacity to approximately 9,400 human recombinant proteins revealed direct binding of TINCR RNA to the staufen1 (STAU1) protein. STAU1-deficient tissue recapitulated the impaired differentiation seen with TINCR depletion. Loss of UPF1 and UPF2, both of which are required for STAU1-mediated RNA decay, however, did not have differentiation effects. Instead, the TINCR–STAU1 complex seems to mediate stabilization of differentiation mRNAs, such as KRT80. These data identify TINCR as a key lncRNA required for somatic tissue differentiation, which occurs through lncRNA binding to differentiation mRNAs to ensure their expression.

Suppression of progenitor differentiation requires the long noncoding RNA ANCR

Long noncoding RNAs (lncRNAs) regulate diverse processes, yet a potential role for lncRNAs in maintaining the undifferentiated state in somatic tissue progenitor cells remains uncharacterized. We used transcriptome sequencing and tiling arrays to compare lncRNA expression in epidermal progenitor populations versus differentiating cells. We identified ANCR (anti-differentiation ncRNA) as an 855-base-pair lncRNA down-regulated during differentiation. Depleting ANCR in progenitor-containing populations, without any other stimuli, led to rapid differentiation gene induction. In epidermis, ANCR loss abolished the normal exclusion of differentiation from the progenitor-containing compartment. The ANCR lncRNA is thus required to enforce the undifferentiated cell state within epidermis.

BRAFV600E remodels the melanocyte transcriptome and induces BANCR to regulate melanoma cell migration

Aberrations of protein-coding genes are a focus of cancer genomics; however, the impact of oncogenes on expression of the ~50% of transcripts without protein-coding potential, including long noncoding RNAs (lncRNAs), has been largely uncharacterized. Activating mutations in the BRAF oncogene are present in >70% of melanomas, 90% of which produce active mutant BRAF(V600E) protein. To define the impacts of oncogenic BRAF on the melanocyte transcriptome, massively parallel cDNA sequencing (RNA-seq) was performed on genetically matched normal human melanocytes with and without BRAF(V600E) expression. To enhance potential disease relevance by verifying expression of altered genes in BRAF-driven cancer tissue, parallel RNA-seq was also undertaken of two BRAF(V600E)-mutant human melanomas. BRAF(V600E) regulated expression of 1027 protein-coding transcripts and 39 annotated lncRNAs, as well as 70 unannotated, potentially novel, intergenic transcripts. These transcripts display both tissue-specific and multi-tissue expression profiles and harbor distinctive regulatory chromatin marks and transcription factor binding sites indicative of active transcription. Coding potential analysis of the 70 unannotated transcripts suggested that most may represent newly identified lncRNAs. BRAF-regulated lncRNA 1 (BANCR) was identified as a recurrently overexpressed, previously unannotated 693-bp transcript on chromosome 9 with a potential functional role in melanoma cell migration. BANCR knockdown reduced melanoma cell migration, and this could be rescued by the chemokine CXCL11. Combining RNA-seq of oncogene-expressing normal cells with RNA-seq of their corresponding human cancers may represent a useful approach to discover new oncogene-regulated RNA transcripts of potential clinical relevance in cancer.

A LncRNA-MAF:MAFB Transcription Factor Network Regulates Epidermal Differentiation

Progenitor differentiation requires remodeling of genomic expression; however, in many tissues, such as epidermis, the spectrum of remodeled genes and the transcription factors (TFs) that control them are not fully defined. We performed kinetic transcriptome analysis during regeneration of differentiated epidermis and identified gene sets enriched in progenitors (594 genes), in early (159 genes), and in late differentiation (387 genes). Module mapping of 1,046 TFs identified MAF and MAFB as necessary and sufficient for progenitor differentiation. MAF:MAFB regulated 393 genes altered in this setting. Integrative analysis identified ANCR and TINCR lncRNAs as essential upstream MAF:MAFB regulators. ChIP-seq analysis demonstrated MAF:MAFB binding to known epidermal differentiation TF genes whose expression they controlled, including GRHL3, ZNF750, KLF4, and PRDM1. Each of these TFs rescued expression of specific MAF:MAFB target gene subsets in the setting of MAF:MAFB loss, indicating they act downstream of MAF:MAFB. A lncRNA-TF network is thus essential for epidermal differentiation.

The noncoding RNAs SNORD50A and SNORD50B bind K-Ras and are recurrently deleted in human cancer

Small nucleolar RNAs (snoRNAs) are conserved noncoding RNAs best studied as ribonucleoprotein (RNP) guides in RNA modification. To explore their role in cancer, we compared 5,473 tumor-normal genome pairs to identify snoRNAs with frequent copy number loss. The SNORD50A-SNORD50B snoRNA locus was deleted in 10–40% of 12 common cancers, where its loss was associated with reduced survival. A human protein microarray screen identified direct SNORD50A and SNORD50B RNA binding to K-Ras. Loss of SNORD50A and SNORD50B increased the amount of GTP-bound, active K-Ras and hyperactivated Ras-ERK1/ERK2 signaling. Loss of these snoRNAs also increased binding by farnesyltransferase to K-Ras and increased K-Ras prenylation, suggesting that KRAS mutation might synergize with SNORD50A and SNORD50B loss in cancer. In agreement with this hypothesis, CRISPR-mediated deletion of SNORD50A and SNORD50B in KRAS-mutant tumor cells enhanced tumorigenesis, and SNORD50A and SNORD50B deletion and oncogenic KRAS mutation co-occurred significantly in multiple human tumor types. SNORD50A and SNORD50B snoRNAs thus directly bind and inhibit K-Ras and are recurrently deleted in human cancer.

Dicer-microRNA-Myc circuit promotes transcription of hundreds of long noncoding RNAs

Long noncoding RNAs (lncRNAs) are important regulators of cell fate, yet little is known about mechanisms controlling lncRNA expression. Here we show that transcription is quantitatively different for lncRNAs and mRNAs—as revealed by deficiency of Dicer (Dcr), a key RNase that generates microRNAs (miRNAs). Dcr loss in mouse embryonic stem cells led unexpectedly to decreased levels of hundreds of lncRNAs. The canonical Dgcr8-Dcr-miRNA pathway is required for robust lncRNA transcriptional initiation and elongation. Computational and genetic epistasis analyses demonstrated that Dcr activation of the oncogenic transcription factor cMyc is partly responsible for lncRNA expression. A quantitative metric of mRNA-lncRNA decoupling revealed that Dcr and cMyc differentially regulate lncRNAs versus mRNAs in diverse cell types and in vivo. Thus, numerous lncRNAs may be modulated as a class in development and disease, notably where Dcr and cMyc act.

Genomic Profiling of a Human Organotypic Model of AEC Syndrome Reveals ZNF750 as an Essential Downstream Target of Mutant TP63

The basis for impaired differentiation in TP63 mutant ankyloblepharon-ectodermal dysplasia-clefting (AEC) syndrome is unknown. Human epidermis harboring AEC TP63 mutants recapitulated this impairment, along with downregulation of differentiation activators, including HOPX, GRHL3, KLF4, PRDM1, and ZNF750. Gene-set enrichment analysis indicated that disrupted expression of epidermal differentiation programs under the control of ZNF750 and KLF4 accounted for the majority of disrupted epidermal differentiation resulting from AEC mutant TP63. Chromatin immunoprecipitation (ChIP) analysis and ChIP-sequencing of TP63 binding in differentiated keratinocytes revealed ZNF750 as a direct target of wild-type and AEC mutant TP63. Restoring ZNF750 to AEC model tissue rescued activator expression and differentiation, indicating that AEC TP63-mediated ZNF750 inhibition contributes to differentiation defects in AEC. Incorporating disease-causing mutants into regenerated human tissue can thus dissect pathomechanisms and identify targets that reverse disease features.

Enhancer-targeted genome editing selectively blocks innate resistance to oncokinase inhibition

Thousands of putative enhancers are characterized in the human genome, yet few have been shown to have a functional role in cancer progression. Inhibiting oncokinases, such as EGFR, ALK, ERBB2, and BRAF, is a mainstay of current cancer therapy but is hindered by innate drug resistance mediated by up-regulation of the HGF receptor, MET. The mechanisms mediating such genomic responses to targeted therapy are unknown. Here, we identify lineage-specific enhancers at the MET locus for multiple common tumor types, including a melanoma lineage-specific enhancer 63 kb downstream from the MET TSS. This enhancer displays inducible chromatin looping with the MET promoter to up-regulate MET expression upon BRAF inhibition. Epigenomic analysis demonstrated that the melanocyte-specific transcription factor, MITF, mediates this enhancer function. Targeted genomic deletion (<7 bp) of the MITF motif within the MET enhancer suppressed inducible chromatin looping and innate drug resistance, while maintaining MITF-dependent, inhibitor-induced melanoma cell differentiation. Epigenomic analysis can thus guide functional disruption of regulatory DNA to decouple pro- and anti-oncogenic functions of a dominant transcription factor and block innate resistance to oncokinase therapy.

Identification of proteins binding coding and non-coding human RNAs using protein microarrays

BackgroundThe regulation and function of mammalian RNAs has been increasingly appreciated to operate via RNA-protein interactions. With the recent discovery of thousands of novel human RNA molecules by high-throughput RNA sequencing, efficient methods to uncover RNA-protein interactions are urgently required. Existing methods to study proteins associated with a given RNA are laborious and require substantial amounts of cell-derived starting material. To overcome these limitations, we have developed a rapid and large-scale approach to characterize binding of in vitro transcribed labeled RNA to ~9,400 human recombinant proteins spotted on protein microarrays.ResultsWe have optimized methodology to probe human protein microarrays with full-length RNA molecules and have identified 137 RNA-protein interactions specific for 10 coding and non-coding RNAs. Those proteins showed strong enrichment for common human RNA binding domains such as RRM, RBD, as well as K homology and CCCH type zinc finger motifs. Previously unknown RNA-protein interactions were discovered using this technique, and these interactions were biochemically verified between TP53 mRNA and Staufen1 protein as well as between HRAS mRNA and CNBP protein. Functional characterization of the interaction between Staufen 1 protein and TP53 mRNA revealed a novel role for Staufen 1 in preserving TP53 RNA stability.ConclusionsOur approach demonstrates a scalable methodology, allowing rapid and efficient identification of novel human RNA-protein interactions using RNA hybridization to human protein microarrays. Biochemical validation of newly identified interactions between TP53-Stau1 and HRAS-CNBP using reciprocal pull-down experiments, both in vitro and in vivo, demonstrates the utility of this approach to study uncharacterized RNA-protein interactions.