Assessing genome-wide dynamic changes in enhancer activity during early mESC differentiation by FAIRE-STARR-seq

Embryonic stem cells (ESCs) can differentiate into any given cell type and therefore represent a versatile model to study the link between gene regulation and differentiation. To quantitatively assess the dynamics of enhancer activity during the early stages of murine ESC differentiation, we analyzed accessible genomic regions using STARR-seq, a massively parallel reporter assay. This resulted in a genome-wide quantitative map of active mESC enhancers, in pluripotency and during the early stages of differentiation. We find that only a minority of accessible regions is active and that such regions are enriched near promoters, characterized by specific chromatin marks, enriched for distinct sequence motifs, and modeling shows that active regions can be predicted from sequence alone. Regions that change their activity upon retinoic acid-induced differentiation are more prevalent at distal intergenic regions when compared to constitutively active enhancers. Further, analysis of differentially active enhancers verified the contribution of individual TF motifs toward activity and inducibility as well as their role in regulating endogenous genes. Notably, the activity of retinoic acid receptor alpha (RARα) occupied regions can either increase or decrease upon the addition of its ligand, retinoic acid, with the direction of the change correlating with spacing and orientation of the RARα consensus motif and the co-occurrence of additional sequence motifs. Together, our genome-wide enhancer activity map elucidates features associated with enhancer activity levels, identifies regulatory regions disregarded by computational prediction tools, and provides a resource for future studies into regulatory elements in mESCs.     

Long ncRNA expression associates with tissue-specific enhancers

Long non-coding RNAs (ncRNA) have recently been demonstrated to be expressed from a subset of enhancers and to be required for the distant regulation of gene expression. Several approaches to predict enhancers have been developed based on various chromatin marks and occupancy of enhancer-binding proteins. Despite the rapid advances in the field, no consensus how to define tissue specific enhancers yet exists. Here, we identify 2,695 long ncRNAs annotated by ENCODE (corresponding to 28% of all ENCODE annotated long ncRNAs) that overlap tissue-specific enhancers. We use a recently developed algorithm to predict tissue-specific enhancers, PreSTIGE, that is based on the H3K4me1 mark and tissue specific expression of mRNAs. The expression of the long ncRNAs overlapping enhancers is significantly higher when the enhancer is predicted as active in a specific cell line, suggesting a general interdependency of active enhancers and expression of long ncRNAs. This dependency is not identified using previous enhancer prediction algorithms that do not account for expression of their downstream targets. The predicted enhancers that overlap annotated long ncRNAs generally have a lower ratio of H3K4me1 to H3K4me3, suggesting that enhancers expressing long ncRNAs might be associated with specific epigenetic marks. In conclusion, we demonstrate the tissue-specific predictive power of PreSTIGE and provide evidence for thousands of long ncRNAs that are expressed from active tissue-specific enhancers, suggesting a particularly important functional relationship between long ncRNAs and enhancer activity in determining tissue-specific gene expression.     

Long ncRNA expression associates with tissue-specific enhancers

Long non-coding RNAs (ncRNA) have recently been demonstrated to be expressed from a subset of enhancers and to be required for the distant regulation of gene expression. Several approaches to predict enhancers have been developed based on various chromatin marks and occupancy of enhancer-binding proteins. Despite the rapid advances in the field, no consensus how to define tissue specific enhancers yet exists. Here, we identify 2,695 long ncRNAs annotated by ENCODE (corresponding to 28% of all ENCODE annotated long ncRNAs) that overlap tissue-specific enhancers. We use a recently developed algorithm to predict tissue-specific enhancers, PreSTIGE, that is based on the H3K4me1 mark and tissue specific expression of mRNAs. The expression of the long ncRNAs overlapping enhancers is significantly higher when the enhancer is predicted as active in a specific cell line, suggesting a general interdependency of active enhancers and expression of long ncRNAs. This dependency is not identified using previous enhancer prediction algorithms that do not account for expression of their downstream targets. The predicted enhancers that overlap annotated long ncRNAs generally have a lower ratio of H3K4me1 to H3K4me3, suggesting that enhancers expressing long ncRNAs might be associated with specific epigenetic marks. In conclusion, we demonstrate the tissue-specific predictive power of PreSTIGE and provide evidence for thousands of long ncRNAs that are expressed from active tissue-specific enhancers, suggesting a particularly important functional relationship between long ncRNAs and enhancer activity in determining tissue-specific gene expression.     

Nuclear receptor response elements mediate induction of intestinal MDR1 by rifampin

Intestinal P-glycoprotein, which is encoded by the MDR1 gene, plays an important role in the absorption and presystemic elimination of many xenobiotics. Hence, an understanding of the factors regulating its expression and function is of substantial interest. In addition to genetic factors, exposure to drugs such as rifampin can profoundly affect its expression. So far, the mechanisms by which rifampin induces MDR1 expression are poorly understood. Recent studies demonstrate that the nuclear receptor PXR (pregnane X receptor) is involved in xenobiotic induction of CYP3A4. Because CYP3A4 and MDR1 are often co-induced, we investigated whether a similar mechanism is also involved in MDR1 induction. The human colon carcinoma cell line LS174T was used as an intestinal model to study induction because in these cells the endogenous MDR1 gene is highly inducible by rifampin. The 5'-upstream region of human MDR1 was examined for the presence of potential PXR response elements. Several binding sites were identified that form a complex regulatory cluster at about -8 kilobase pairs. Only one DR4 motif within this cluster is necessary for induction by rifampin. We conclude that induction of MDR1 is mediated by a DR4 motif in the upstream enhancer at about -8 kilobase pairs, to which PXR binds.     

通过系统分析揭示组蛋白修饰模式与转录因子结合之间的关联

转录因子对顺势调控元件的选择性结合,在哺乳动物细胞类型特异的基因表达中扮演重要的角色.这个过程受到染色质表观遗传状态的潜在调控.近期,染色质免疫共沉淀结合测序的研究提供了大量泛基因组水平的数据,阐述转录因子结合以及组蛋白修饰的位点,这为系统地研究转录因子和表观遗传标记之间的空间及调控关系提供了基础.该研究对公共数据库中的染色质免疫共沉淀结合测序数据进行整合分析,涉及5个细胞系中的85种转录因子、9种组蛋白修饰,目的是研究转录因子结合位点与组蛋白修饰模式以及基因表达在泛基因组水平上的关联.作者发现,不同转录因子与组蛋白修饰的共定位模式高度一致,并且组蛋白修饰在距离转录因子结合位点约500碱基对的位置富集.作者还发现,转录因子结合位点的占有率与活性组蛋白修饰的水平和双峰模式正相关,并且启动子区域组蛋白修饰的双峰和共定位模式和基因的高转录水平相一致.组蛋白修饰模式、转录因子结合位点的占有率与基因转录之间的关联暗示了细胞可能利用的基因表达调控机制.     

Regulation of CYP3A4 by pregnane X receptor: The role of nuclear receptors competing for response element binding

Induction of the major drug metabolizing enzyme CYP3A4 by xenobiotics contributes to the pronounced interindividual variability of its expression and often results in clinically relevant drug-drug interactions. It is mainly mediated by PXR, which regulates CYP3A4 expression by binding to several specific elements in the 5′ upstream regulatory region of the gene. Induction itself shows a marked interindividual variability, whose underlying determinants are only partly understood. In this study, we investigated the role of nuclear receptor binding to PXR response elements in CYP3A4, as a potential non-genetic mechanism contributing to interindividual variability of induction. By in vitro DNA binding experiments, we showed that several nuclear receptors bind efficiently to the proximal promoter ER6 and distal xenobiotic-responsive enhancer module DR3 motifs. TRα1, TRβ1, COUP-TFI, and COUP-TFII further demonstrated dose-dependent repression of PXR-mediated CYP3A4 enhancer/promoter reporter activity in transient transfection in the presence and absence of the PXR inducer rifampin, while VDR showed this effect only in the absence of treatment. By combining functional in vitro characterization with hepatic expression analysis, we predict that TRα1, TRβ1, COUP-TFI, and COUP-TFII show a strong potential for the repression of PXR-mediated activation of CYP3A4 in vivo. In summary, our results demonstrate that nuclear receptor binding to PXR response elements interferes with PXR-mediated expression and induction of CYP3A4 and thereby contributes to the interindividual variability of induction.     

Retinoic acid-mediated down-regulation of Oct3/4 coincides with the loss of promoter occupancy in vivo.

Oct3/4, a hallmark of the earliest stages of embryogenesis, is expressed in undifferentiated embryonal carcinoma (EC) and embryonic stem (ES) cells. Oct3/4 gene expression is dependent on the promoter region, the proximal enhancer and the newly identified distal enhancer. We have analysed in vivo occupancy of these elements. In undifferentiated EC and ES cells, strong footprints were detected at specific sites of all three regulatory elements. These were promptly lost upon RA treatment in ES cells and in P19 EC cells, in parallel with sharply reduced Oct3/4 mRNA levels. Thus, the occupancy of regulatory elements is coupled with Oct3/4 expression, and RA treatment causes coordinated factor displacement, leading to extinction of gene activity. In F9 EC cells, footprint was first abolished at the proximal enhancer. However, this loss of binding site occupancy did not result in a decrease in Oct3/4 mRNA levels. The partial factor displacement seen in F9 EC cells, combined with the observation that EC and ES cells utilize the proximal and distal enhancers in differential manner, indicate the complex pattern of Oct3/4 gene regulation, which could reflect a cell type- and lineage-specific expression of the gene in vivo.