Genome-wide analysis of RNA-chromatin interactions in lizards as a mean for functional lncRNA identification.

BACKGROUND: Long non-coding RNAs (lncRNAs) are defined as transcribed molecules longer than 200 nucleotides with little to no protein-coding potential. LncRNAs can regulate gene expression of nearby genes (cis-acting) or genes located on other chromosomes (trans-acting). Several methodologies have been developed to capture lncRNAs associated with chromatin at a genome-wide level. Analysis of RNA-DNA contacts can be combined with epigenetic and RNA-seq data to define potential lncRNAs involved in the regulation of gene expression. RESULTS: We performed Chromatin Associated RNA sequencing (ChAR-seq) in Anolis carolinensis to obtain the genome-wide map of the associations that RNA molecules have with chromatin. We analyzed the frequency of DNA contacts for different classes of RNAs and were able to define cis- and trans-acting lncRNAs. We integrated the ChAR-seq map of RNA-DNA contacts with epigenetic data for the acetylation of lysine 16 on histone H4 (H4K16ac), a mark connected to actively transcribed chromatin in lizards. We successfully identified three trans-acting lncRNAs significantly associated with the H4K16ac signal, which are likely involved in the regulation of gene expression in A. carolinensis. CONCLUSIONS: We show that the ChAR-seq method is a powerful tool to explore the RNA-DNA map of interactions. Moreover, in combination with epigenetic data, ChAR-seq can be applied in non-model species to establish potential roles for predicted lncRNAs that lack functional annotations.

ELAV links paused Pol II to alternative polyadenylation in the Drosophila nervous system.

Alternative polyadenylation (APA) has been implicated in a variety of developmental and disease processes. A particularly dramatic form of APA occurs in the developing nervous system of flies and mammals, whereby various developmental genes undergo coordinate 3' UTR extension. In Drosophila, the RNA-binding protein ELAV inhibits RNA processing at proximal polyadenylation sites, thereby fostering the formation of exceptionally long 3' UTRs. Here, we present evidence that paused Pol II promotes recruitment of ELAV to extended genes. Replacing promoters of extended genes with heterologous promoters blocks normal 3' extension in the nervous system, while extension-associated promoters can induce 3' extension in ectopic tissues expressing ELAV. Computational analyses suggest that promoter regions of extended genes tend to contain paused Pol II and associated cis-regulatory elements such as GAGA. ChIP-seq assays identify ELAV in the promoter regions of extended genes. Our study provides evidence for a regulatory link between promoter-proximal pausing and APA.

Histone H2A deubiquitinase activity of the Polycomb repressive complex PR-DUB.

Polycomb group (PcG) proteins are transcriptional repressors that control processes ranging from the maintenance of cell fate decisions and stem cell pluripotency in animals to the control of flowering time in plants. In Drosophila, genetic studies identified more than 15 different PcG proteins that are required to repress homeotic (HOX) and other developmental regulator genes in cells where they must stay inactive. Biochemical analyses established that these PcG proteins exist in distinct multiprotein complexes that bind to and modify chromatin of target genes. Among those, Polycomb repressive complex 1 (PRC1) and the related dRing-associated factors (dRAF) complex contain an E3 ligase activity for monoubiquitination of histone H2A (refs 1-4). Here we show that the uncharacterized Drosophila PcG gene calypso encodes the ubiquitin carboxy-terminal hydrolase BAP1. Biochemically purified Calypso exists in a complex with the PcG protein ASX, and this complex, named Polycomb repressive deubiquitinase (PR-DUB), is bound at PcG target genes in Drosophila. Reconstituted recombinant Drosophila and human PR-DUB complexes remove monoubiquitin from H2A but not from H2B in nucleosomes. Drosophila mutants lacking PR-DUB show a strong increase in the levels of monoubiquitinated H2A. A mutation that disrupts the catalytic activity of Calypso, or absence of the ASX subunit abolishes H2A deubiquitination in vitro and HOX gene repression in vivo. Polycomb gene silencing may thus entail a dynamic balance between H2A ubiquitination by PRC1 and dRAF, and H2A deubiquitination by PR-DUB.

The role of the histone H2A ubiquitinase Sce in Polycomb repression.

Polycomb group (PcG) proteins exist in multiprotein complexes that modify chromatin to repress transcription. Drosophila PcG proteins Sex combs extra (Sce; dRing) and Posterior sex combs (Psc) are core subunits of PRC1-type complexes. The Sce:Psc module acts as an E3 ligase for monoubiquitylation of histone H2A, an activity thought to be crucial for repression by PRC1-type complexes. Here, we created an Sce knockout allele and show that depletion of Sce results in loss of H2A monoubiquitylation in developing Drosophila. Genome-wide profiling identified a set of target genes co-bound by Sce and all other PRC1 subunits. Analyses in mutants lacking individual PRC1 subunits reveals that these target genes comprise two distinct classes. Class I genes are misexpressed in mutants lacking any of the PRC1 subunits. Class II genes are only misexpressed in animals lacking the Psc-Su(z)2 and Polyhomeotic (Ph) subunits but remain stably repressed in the absence of the Sce and Polycomb (Pc) subunits. Repression of class II target genes therefore does not require Sce and H2A monoubiquitylation but might rely on the ability of Psc-Su(z)2 and Ph to inhibit nucleosome remodeling or to compact chromatin. Similarly, Sce does not provide tumor suppressor activity in larval tissues under conditions in which Psc-Su(z)2, Ph and Pc show such activity. Sce and H2A monoubiquitylation are therefore only crucial for repression of a subset of genes and processes regulated by PRC1-type complexes. Sce synergizes with the Polycomb repressive deubiquitinase (PR-DUB) complex to repress transcription at class I genes, suggesting that H2A monoubiquitylation must be appropriately balanced for their transcriptional repression.

Isolation of Boechera stricta Developing Embryos for Hi-C.

The possibility of analyzing chromatin topology in developing plant embryos is hampered by inaccessibility of the embryo sac, deeply embedded in the maternal seed tissue, following double fertilization. Here we describe a protocol to isolate, purify, and prepare developing Boechera stricta embryos for chromosome conformation capture-based methods as in situ Hi-C experiments. Early globular embryos can be isolated by air-pressure microaspiration, and subsequently washed to eliminate residual cells from the endosperm and maternal seed coat, allowing for pure sampling of selected stages of embryogenesis. This protocol allows for the possibility of comparing genome topology during plant embryonic differentiation since early until late embryo development stages.

Plant In Situ Hi-C Experimental Protocol and Bioinformatic Analysis.

Hi-C enables the characterization of the 0conformation of the genome in the three-dimensional nuclear space. This technique has revolutionized our ability to detect interactions between linearly distant genomic sites on a genome-wide scale. Here, we detail a protocol to carry out in situ Hi-C in plants and describe a straightforward bioinformatics pipeline for the analysis of such data, in particular for comparing samples from different organs or conditions.

Evolution of Genome-Organizing Long Non-coding RNAs in Metazoans.

Long non-coding RNAs (lncRNAs) have important regulatory functions across eukarya. It is now clear that many of these functions are related to gene expression regulation through their capacity to recruit epigenetic modifiers and establish chromatin interactions. Several lncRNAs have been recently shown to participate in modulating chromatin within the spatial organization of the genome in the three-dimensional space of the nucleus. The identification of lncRNA candidates is challenging, as it is their functional characterization. Conservation signatures of lncRNAs are different from those of protein-coding genes, making identifying lncRNAs under selection a difficult task, and the homology between lncRNAs may not be readily apparent. Here, we review the evidence for these higher-order genome organization functions of lncRNAs in animals and the evolutionary signatures they display.

Chromatin Loop Formation Induced by a Subtelomeric Protosilencer Represses EPA Genes in Candida glabrata.

Adherence, an important virulence factor, is mediated by the EPA (Epithelial Adhesin) genes in the opportunistic pathogen Candida glabrata Expression of adhesin-encoding genes requires tight regulation to respond to harsh environmental conditions within the host. The majority of EPA genes are localized in subtelomeric regions regulated by subtelomeric silencing, which depends mainly on Rap1 and the Sir proteins. In vitro adhesion to epithelial cells is primarily mediated by Epa1. EPA1 forms a cluster with EPA2 and EPA3 in the right telomere of chromosome E (E-R). This telomere contains a cis-acting regulatory element, the protosilencer Sil2126 between EPA3 and the telomere. Interestingly, Sil2126 is only active in the context of its native telomere. Replacement of the intergenic regions between EPA genes in E-R revealed that cis-acting elements between EPA2 and EPA3 are required for Sil2126 activity when placed 32 kb away from the telomere (Sil@-32kb). Sil2126 contains several putative binding sites for Rap1 and Abf1, and its activity depends on these proteins. Indeed, Sil2126 binds Rap1 and Abf1 at its native position and also when inserted at -32 kb, a silencing-free environment in the parental strain. In addition, we found that Sil@-32kb and Sil2126 at its native position can physically interact with the intergenic regions between EPA1-EPA2 and EPA2-EPA3 respectively, by chromosome conformation capture assays. We speculate that Rap1 and Abf1 bound to Sil2126 can recruit the Silent Information Regulator complex, and together mediate silencing in this region, probably through the formation of a chromatin loop.

The tumour suppressor L(3)mbt inhibits neuroepithelial proliferation and acts on insulator elements.

In Drosophila, defects in asymmetric cell division often result in the formation of stem-cell-derived tumours. Here, we show that very similar terminal brain tumour phenotypes arise through a fundamentally different mechanism. We demonstrate that brain tumours in l(3)mbt mutants originate from overproliferation of neuroepithelial cells in the optic lobes caused by derepression of target genes in the Salvador-Warts-Hippo (SWH) pathway. We use ChIP-sequencing to identify L(3)mbt binding sites and show that L(3)mbt binds to chromatin insulator elements. Mutating l(3)mbt or inhibiting expression of the insulator protein gene mod(mdg4) results in upregulation of SWH pathway reporters. As l(3)mbt tumours are rescued by mutations in bantam or yorkie or by overexpression of Expanded, the deregulation of SWH pathway target genes is an essential step in brain tumour formation. Therefore, very different primary defects result in the formation of brain tumours, which behave quite similarly in their advanced stages.

Essential role of the glycosyltransferase sxc/Ogt in polycomb repression.

Polycomb group proteins are conserved transcriptional repressors that control animal and plant development. Here, we found that the Drosophila Polycomb group gene super sex combs (sxc) encodes Ogt, the highly conserved glycosyltransferase that catalyzes the addition of N-acetylglucosamine (GlcNAc) to proteins in animals and plants. Genome-wide profiling in Drosophila revealed that GlcNAc-modified proteins are highly enriched at Polycomb response elements. Among different Polycomb group proteins, Polyhomeotic is glycosylated by Sxc/Ogt in vivo. sxc/Ogt-null mutants lacked O-linked GlcNAcylation and failed to maintain Polycomb transcriptional repression even though Polycomb group protein complexes were bound at their target sites. Polycomb repression appears to be a critical function of Sxc/Ogt in Drosophila and may be mediated by the glycosylation of Polyhomeotic.

Dynamic regulation by polycomb group protein complexes controls pattern formation and the cell cycle in Drosophila.

Polycomb group (PcG) proteins form conserved regulatory complexes that modify chromatin to repress transcription. Here, we report genome-wide binding profiles of PhoRC, the Drosophila PcG protein complex containing the DNA-binding factor Pho/dYY1 and dSfmbt. PhoRC constitutively occupies short Polycomb response elements (PREs) of a large set of developmental regulator genes in both embryos and larvae. The majority of these PREs are co-occupied by the PcG complexes PRC1 and PRC2. Analysis of PcG mutants shows that the PcG system represses genes required for anteroposterior, dorsoventral, and proximodistal patterning of imaginal discs and that it also represses cell cycle regulator genes. Many of these genes are regulated in a dynamic manner, and our results suggest that the PcG system restricts signaling-mediated activation of target genes to appropriate cells. Analysis of cell cycle regulators indicates that the PcG system also dynamically modulates the expression levels of certain genes, providing a possible explanation for the tumor phenotype of PcG mutants.

Human papillomavirus E6/E7 oncogenes promote mouse ear regeneration by increasing the rate of wound re-epithelization and epidermal growth.

Mammals have limited regeneration capacity. We report here that, in transgenic mice (Tg(bK6-E6/E7)), the expression of the E6/E7 oncogenes of human papilloma virus type 16 (HPV16) under the control of the bovine keratin 6 promoter markedly improves the mouse's capacity to repair portions of the ear after being wounded. Increased repair capacity correlates with an increased number of epidermal proliferating cells. In concordance with the expected effects of the E6 and E7 oncogenes, levels of p53 decreased and those of p16 in epidermal cells increased. In addition, we observed that wound re-epithelization proceeded faster in transgenic than in wild-type animals. After the initial re-epithelization, epidermal cell migration from the intact surrounding tissue appears to be a major contributor to the growing epidermis, especially in the repairing tissue of transgenic mice. We also found that there is a significantly higher number of putative epidermal stem cells in Tg(bK6-E6/E7) than in wild-type mice. Remarkably, hair follicles and cartilage regenerated within the repaired ear tissue, without evidence of tumor formation. We propose that the ability to regenerate ear portions is limited by the capacity of the epidermis to repair itself and grow.

Pcl-PRC2 is needed to generate high levels of H3-K27 trimethylation at Polycomb target genes.

PRC2 is thought to be the histone methyltransferase (HMTase) responsible for H3-K27 trimethylation at Polycomb target genes. Here we report the biochemical purification and characterization of a distinct form of Drosophila PRC2 that contains the Polycomb group protein polycomblike (Pcl). Like PRC2, Pcl-PRC2 is an H3-K27-specific HMTase that mono-, di- and trimethylates H3-K27 in nucleosomes in vitro. Analysis of Drosophila mutants that lack Pcl unexpectedly reveals that Pcl-PRC2 is required to generate high levels of H3-K27 trimethylation at Polycomb target genes but is dispensable for the genome-wide H3-K27 mono- and dimethylation that is generated by PRC2. In Pcl mutants, Polycomb target genes become derepressed even though H3-K27 trimethylation at these genes is only reduced and not abolished, and even though targeting of the Polycomb protein complexes PhoRC and PRC1 to Polycomb response elements is not affected. Pcl-PRC2 is thus the HMTase that generates the high levels of H3-K27 trimethylation in Polycomb target genes that are needed to maintain a Polycomb-repressed chromatin state.