Long noncoding RNA programs active chromatin domain to coordinate homeotic gene activation
Kevin C. Wang,Yul W. Yang,Bo Liu,Amartya Sanyal,Ryan Corces-Zimmerman,Yong Chen,Bryan R. Lajoie,Angeline Protacio,Ryan A. Flynn,Rajnish A. Gupta,Joanna Wysocka,Ming Lei,Job Dekker,Jill A. Helms,Howard Y. Chang
2013-01-01
Abstract:The genome is extensively transcribed into long intergenic noncoding RNAs (lincRNAs), many of which are implicated in gene silencing1,2. Potential roles of lincRNAs in gene activation are much less understood3,4,5. Development and homeostasis require coordinate regulation of neighboring genes through a process termed locus control6. Some locus control elements and enhancers transcribe lincRNAs7,8,9,10, hinting at possible roles in long range control. In vertebrates, 39 Hox genes, encoding homeodomain transcription factors critical for positional identity, are clustered in four chromosomal loci; the Hox genes are expressed in nested anterior-posterior and proximaldistal patterns co-linear with their genomic position from 3′ to 5′of the cluster11. Here we identify HOTTIP, a lincRNA transcribed from the 5′ tip of the HOXA locus that coordinates the activation of multiple 5′ HOXA genes in vivo. Chromosomal looping brings HOTTIP into close proximity to its target genes. HOTTIP directly binds the adaptor protein WDR5 and targets WDR5/MLL complexes across HOXA, driving histone H3 lysine 4 trimethylation and gene transcription. Induced proximity is necessary and sufficient for HOTTIP activation of its target genes. Thus, by serving as key intermediates that transmit information from higher order chromosomal looping †To whom correspondence should be addressed. howchang@stanford.edu. *These authors contributed equally to this work. Supplementary Information is linked to the online version of the paper at www.nature.com/nature. Author Contributions K.C.W., R.A.G., and H.Y.C. initiated the project; K.C.W. and H.Y.C. designed the experiments; K.C.W., Y.W.Y., B.L., A.S., R.C.Z., B.L., A.P., R.A.F., J.D., and J.A.H. conducted the experiments and analyzed the data; Y.C. and M.L. purified the recombinant proteins; J.W. provided antibodies and cell lines; K.C.W. and H.Y.C. prepared the manuscript with inputs from all co-authors. Author Information: Sequence for human HOTTIP has been deposited with GenBank under the accession number GU724873. Reprints and permissions information is available at www.nature.com/nature/reprints. The authors declare no competing financial interests. NIH Public Access Author Manuscript Nature. Author manuscript; available in PMC 2013 June 03. Published in final edited form as: Nature. 2011 April 7; 472(7341): 120–124. doi:10.1038/nature09819. N IH PA Athor M anscript N IH PA Athor M anscript N IH PA Athor M anscript into chromatin modifications, lincRNAs may organize chromatin domains to coordinate longrange gene activation. We examined chromosome structure and histone modifications in primary fibroblasts derived from several anatomic sites12, and found distinctive differences in the HOXA locus. High throughput chromosome conformation capture (5C)13 across HOXA revealed that its higher order structure is dependent on positional identity. In anatomically distal cells (e.g. foreskin and foot fibroblasts), we detected abundant chromatin interactions within the transcriptionally active 5′ HOXA locus (with reference to the directions of transcription of constituent Hox genes), pointing to a compact and looped conformation. In contrast, no long-range chromatin interactions are detected within the transcriptionally silent 3′ HOXA which appears largely linear (Fig. 1A). Strikingly, anatomically proximal cells (e.g. lung fibroblasts) have the diametrically opposite pattern. The ON and OFF states of Hox and other key developmental genes are maintained by the MLL/Trithorax (Trx) and Polycomb group (PcG) proteins, which mediate trimethylation of histone H3 lysine 4 (H3K4me3) to activate genes or lysine 27 (H3K27me3) to repress genes14. The portions of HOXA in tight physical interaction are marked by broad domains of H3K4me3, whereas H3K27me3 marks the physically extended and transcriptional silent regions (Fig. 1A). On the very 5′ and 3′ edges of the two respective interaction clusters are two lincRNA loci that exhibit distinct chromatin modifications. The 3′element has been previously identified as the myelopoiesis-associated lincRNA HOTAIRM115. The 5′ element, for which we suggest the name HOTTIP for “HOXA transcript at the distal tip”, exhibits bivalent H3K4me3 and H3K27me3, a histone modification pattern associated with poised regulatory sequences16. (Henceforth we refer to the DNA element with italics and the RNA without italics). Comparison with RNA polymerase II occupancy and RNA expression showed that the bivalent H3K4me3 and H3K27me3 modifications on HOTTIP do not require HOTTIP transcription, but transcription of HOTTIP is associated with increased H3K4me3 and decreased H3K27me3 at HOTTIP itself (Fig. 1A, left). ChIP analysis confirmed that HOTTIP is occupied by both PRC2 and MLL complexes, consistent with the bivalent histone marks (Fig. S1A). HOTTIP transcription yields a 3764 nucleotide, spliced, and polyadenylated lincRNA that initiates ~330 bases upstream of HOXA13. Only the strand antisense to HOXA genes is transcribed (Fig. S1B). Genes near the 5′ end of each HOX cluster tend to be expressed in more posterior and/or distal anatomical locations. Consistent with its genomic location 5′ to HOXA13, HOTTIP is expressed in distal and/or posterior anatomic sites, (Fig. 1B). In situ hybridization of developing mouse and chick embryos confirmed that HOTTIP is expressed in posterior and distal sites in vivo, suggesting a conserved expression pattern from development to adulthood (Fig. 1C and Fig. S1C). Even in distal cells where HOTTIP is expressed, its RNA level is very low and estimated to be ~0.3 copies per cell (Fig. S2). We employed small interfering RNAs (siRNAs) to knock down HOTTIP RNA in fibroblasts from a distal anatomic site (foreskin), and examined expression of 5′ HOXA genes by quantitative reverse transcription PCR. Notably, HOTTIP knockdown abrogated expression of distal HOXA genes across 40 kilobases with a trend dependent on the distance to HOTTIP. The strongest blockade was observed for HOXA13 and HOXA11, with progressively less severe effects on HOXA10, HOXA9, and HOXA7 (Fig. 2A). The effect on gene transcription appeared to be unidirectional, as there were no appreciable changes in the levels of EVX1, located ~40 kilobases 5′ of the HOXA cluster (data not shown). HOTTIP knockdown did not affect expression of the highly homologous HOXD genes, other control genes, nor induce antisense transcription at its own locus (Fig. 2B, Fig. S3A). Wang et al. Page 2 Nature. Author manuscript; available in PMC 2013 June 03. N IH PA Athor M anscript N IH PA Athor M anscript N IH PA Athor M anscript Multiple independent siRNAs targeting HOTTIP yielded similar results (Fig. S3B). These results suggest that HOTTIP RNA is necessary to coordinate activation of 5′ HOXA genes. We next addressed the function of HOTTIP RNA in vivo in the developing chick limb bud (Fig. 2C). While prior genetic studies of ncRNAs involved deletion or insertion into the gene locus17, we wished to distinguish the functions of HOTTIP RNA from its corresponding DNA element. HOTTIP can nucleate H3K4 and H3K27 methylation independent of transcription (Fig. 1A), and the precise genomic distance between upstream enhancer elements and Hox genes is critical for their proper co-linear activation17. Therefore, we employed RNAi in chick embryos, where replication-competent retroviruses can deliver short-hairpin RNAs (shRNAs) with high penetrance and precise spatiotemporal control18 (Fig. S4). In the limb bud, 5′ HoxA genes are transcribed in a nested pattern along the proximal-distal axis19. In this tissue, HoxA function is highly redundant with that of the HoxD locus, which allowed us to assess altered HoxA expression patterns without major changes in anatomic landmarks20. We injected retroviruses carrying shRNAs against chick HOTTIP into upper limb buds of stage 13 chicks; RT-PCR and in situ hybridization were performed on both control and knockdown samples after 2-4 days. Knockdown of HOTTIP by two independent shRNAs in limb buds decreased expression of HoxA13, HoxA11, and HoxA10—again with a graded impact depending on genomic proximity to HOTTIP. Vector control or an shRNA that fails to deplete HOTTIP had little effect on Hox gene expression (Fig. 2D). In situ hybridization on whole embryos (Fig. 2E) and sections (Fig. S5) revealed that HOTTIP most strongly affects HoxA gene expression at the distal edge of the developing limb bud, where the 5′ HoxA genes are most strongly expressed. By stage 36, limbs depleted of HOTTIP showed notable shortening and bending of distal bony elements, including the radius, ulna, and third digit (~20% length reduction for each compared to contralateral and stage-matched limbs treated with control virus, p < 0.05, Student′s t-test, Fig. 2F). This phenotype resembled some of the defects in mice lacking HoxA11 and HoxA1321,22,23. Together, these data indicate that HOTTIP RNA controls activation of distal Hox genes in vivo. The broad impact of HOTTIP on gene activation across HOXA locus is reminiscent of the broad domains of chromatin modifications demarcating active and silent chromosomal domains12. 5C analysis of control and HOTTIP-depleted cells showed little change in higher order chromosomal structure, suggesting that the chromosomal looping is pre-configured and upstream of gene expression (Fig. S6A). In contrast, HOTTIP knockdown led to broad loss of H3K4me3 and H3K4me2 across the HOXA locus, mostly prominently over 5′ HOXA and HOTTIP itself (Fig. 3A, Figs. S6B and S7). HOTTIP knockdown also increased H3K27me3 focally over HOTTIP, but had little impact on H3K27me3 across HOXA. These results indicate that HOTTIP RNA is required for maintenance of H3K4me3 across the HOXA. These findings also imply that loss of 5′ HOXA gene transcription upon HOTTIP knockdown is likely due to loss of H3K4me3 (or other changes) rather than ectopic spread of H3K27me3. H3K4 methylation of the HOX loci is carried out by the ML