Identification of DNA regions and a set of transcriptional regulatory factors involved in transcriptional regulation of several human liver-enriched transcription factor genes

Transcriptomic Analysis of Hepatitis B Infected Liver for Prediction of Hepatocellular Carcinoma

Hepatocellular cancer (HCC) is a leading cause of cancer-related mortality worldwide, and chronic hepatitis B virus infection (CHB) has been a major risk factor for HCC development. The pathogenesis of HBV-related HCC has been a major focus revealing the interplay of a multitude of intracellular signaling pathways, yet the precise mechanisms and their implementations to clinical practice remain to be elucidated. This study utilizes publicly available transcriptomic data from the livers of CHB patients in order to identify a population with a higher risk of malignant transformation. We report the identification of a novel list of genes (PCM1) which can generate clear transcriptomic sub-groups among HBV-infected livers. PCM1 includes genes related to cell cycle activity and liver cancer development. In addition, markers of inflammation, M1 macrophages and gamma delta T cell infiltration are present within the signature. Genes within PCM1 are also able to differentiate HCC from normal liver, and some genes within the signature are associated with poor prognosis of HCC at the mRNA level. The analysis of the immunohistochemical stainings validated that proteins coded by a group of PCM1 genes were overexpressed in liver cancer, while minimal or no expression was detected in normal liver. Altogether, our findings suggest that PCM1 can be developed into a clinically applicable method to identify CHB patients with a higher risk of HCC development.

SNAIL1-mediated downregulation of FOXA proteins facilitates the inactivation of transcriptional enhancer elements at key epithelial genes in colorectal cancer cells

Phenotypic conversion of tumor cells through epithelial-mesenchymal transition (EMT) requires massive gene expression changes. How these are brought about is not clear. Here we examined the impact of the EMT master regulator SNAIL1 on the FOXA family of transcription factors which are distinguished by their particular competence to induce chromatin reorganization for the activation of transcriptional enhancer elements. We show that the expression of SNAIL1 and FOXA genes is anticorrelated in transcriptomes of colorectal tumors and cell lines. In cellular EMT models, ectopically expressed Snail1 directly represses FOXA1 and triggers downregulation of all FOXA family members, suggesting that loss of FOXA expression promotes EMT. Indeed, cells with CRISPR/Cas9-induced FOXA-deficiency acquire mesenchymal characteristics. Furthermore, ChIP-seq data analysis of FOXA chromosomal distribution in relation to chromatin structural features which characterize distinct states of transcriptional activity, revealed preferential localization of FOXA factors to transcriptional enhancers at signature genes that distinguish epithelial from mesenchymal colon tumors. To validate the significance of this association, we investigated the impact of FOXA factors on structure and function of enhancers at the CDH1, CDX2 and EPHB3 genes. FOXA-deficiency and expression of dominant negative FOXA2 led to chromatin condensation at these enhancer elements. Site-directed mutagenesis of FOXA binding sites in reporter gene constructs and by genome-editing in situ impaired enhancer activity and completely abolished the active chromatin state of the EPHB3 enhancer. Conversely, expression of FOXA factors in cells with inactive CDX2 and EPHB3 enhancers led to chromatin opening and de novo deposition of the H3K4me1 and H3K27ac marks. These findings establish the pioneer function of FOXA factors at enhancer regions of epithelial genes and demonstrate their essential role in maintaining enhancer structure and function. Thus, by repressing FOXA family members, SNAIL1 targets transcription factors at strategically important positions in gene-regulatory hierarchies, which may facilitate transcriptional reprogramming during EMT.

Cancer patient mortality is overwhelmingly due to distant organ metastases. Epithelial-mesenchymal transition is a process thought to facilitate local invasion and dissemination of cancer cells, thereby promoting metastasis. The conversion of epithelial cells into mesenchymal, fibroblast-like cells requires profound gene expression changes. A few transcription factors like SNAIL1 can initiate these changes, but are unlikely to be solely responsible for all of them. In our study we asked, whether destabilization of epithelial gene expression programs could involve FOXA transcription factors. FOXA factors represent a special subgroup of regulatory proteins, so-called pioneer factors, with unique roles in the activation of transcriptional enhancers which are key regulatory DNA elements that orchestrate spatio-temporal gene expression. In a model of colorectal cancer we found that SNAIL1 represses FOXA factors, and demonstrate that FOXA factors are associated with enhancer elements at epithelial signature genes. Indeed, FOXA factors are sufficient to initiate enhancer activation and necessary to maintain their activity. Our findings indicate that SNAIL1 induces pervasive repression of epithelial genes through a hierarchical scheme of alterations in transcription factor expression which may be applicable to other instances of cell fate changes and transcriptional reprogramming.

Epigenetic modification of histone 3 lysine 27: mediator subunit MED25 is required for the dissociation of polycomb repressive complex 2 from the promoter of cytochrome P450 2C9

The Mediator complex is vital for the transcriptional regulation of eukaryotic genes. Mediator binds to nuclear receptors at target response elements and recruits chromatin-modifying enzymes and RNA polymerase II. Here, we examine the involvement of Mediator subunit MED25 in the epigenetic regulation of human cytochrome P450 2C9 (CYP2C9). MED25 is recruited to the CYP2C9 promoter through association with liver-enriched HNF4α, and we show that MED25 influences the H3K27 status of the HNF4α binding region. This region was enriched for the activating marker H3K27ac and histone acetyltransferase CREBBP after MED25 overexpression but was trimethylated when MED25 expression was silenced. The epigenetic regulator Polycomb repressive complex (PRC2), which represses expression by methylating H3K27, plays an important role in target gene regulation. Silencing MED25 correlated with increased association of PRC2 not only with the promoter region chromatin but with HNF4α itself. We confirmed the involvement of MED25 for fully functional preinitiation complex recruitment and transcriptional output in vitro. Formaldehyde-assisted isolation of regulatory elements (FAIRE) revealed chromatin conformation changes that were reliant on MED25, indicating that MED25 induced a permissive chromatin state that reflected increases in CYP2C9 mRNA. For the first time, we showed evidence that a functionally relevant human gene is transcriptionally regulated by HNF4α via MED25 and PRC2. CYP2C9 is important for the metabolism of many exogenous chemicals including pharmaceutical drugs as well as endogenous substrates. Thus, MED25 is important for regulating the epigenetic landscape resulting in transcriptional activation of a highly inducible gene, CYP2C9.

Clustering of tissue-specific sub-TADs accompanies the regulation of HoxA genes in developing limbs

HoxA genes exhibit central roles during development and causal mutations have been found in several human syndromes including limb malformation. Despite their importance, information on how these genes are regulated is lacking. Here, we report on the first identification of bona fide transcriptional enhancers controlling HoxA genes in developing limbs and show that these enhancers are grouped into distinct topological domains at the sub-megabase scale (sub-TADs). We provide evidence that target genes and regulatory elements physically interact with each other through contacts between sub-TADs rather than by the formation of discreet “DNA loops”. Interestingly, there is no obvious relationship between the functional domains of the enhancers within the limb and how they are partitioned among the topological domains, suggesting that sub-TAD formation does not rely on enhancer activity. Moreover, we show that suppressing the transcriptional activity of enhancers does not abrogate their contacts with HoxA genes. Based on these data, we propose a model whereby chromatin architecture defines the functional landscapes of enhancers. From an evolutionary standpoint, our data points to the convergent evolution of HoxA and HoxD regulation in the fin-to-limb transition, one of the major morphological innovations in vertebrates.

Hox genes encode transcription factors with crucial roles during development. These genes are grouped in four different clusters names HoxA, B, C, and D. Mutations in genes of the HoxA and D clusters have been found in several human syndromes, affecting in some cases limb development. Despite their essential role and contrary to the genes of the HoxD cluster, little is known about how the HoxA genes are regulated. Here, we identified a large set of regulatory elements controlling HoxA genes during limb development. By studying spatial chromatin organization at the HoxA region, we found that the regulatory elements are spatially clustered regardless of their activity. Clustering of enhancers define tissue-specific chromatin domains that interact specifically with each other and with active genes in the limb. Our findings give support to the emerging concept that chromatin architecture defines the functional properties of genomes. Additionally, our study suggests a common constraint of the chromatin topology in the evolution of HoxA and HoxD regulation in the emergence of the hand/foot, which is one of the major morphological innovations in vertebrates.

Statistical significance of combinatorial regulations

More than three transcription factors often work together to enable cells to respond to various signals. The detection of combinatorial regulation by multiple transcription factors, however, is not only computationally nontrivial but also extremely unlikely because of multiple testing correction. The exponential growth in the number of tests forces us to set a strict limit on the maximum arity. Here, we propose an efficient branch-and-bound algorithm called the "limitless arity multiple-testing procedure" (LAMP) to count the exact number of testable combinations and calibrate the Bonferroni factor to the smallest possible value. LAMP lists significant combinations without any limit, whereas the family-wise error rate is rigorously controlled under the threshold. In the human breast cancer transcriptome, LAMP discovered statistically significant combinations of as many as eight binding motifs. This method may contribute to uncover pathways regulated in a coordinated fashion and find hidden associations in heterogeneous data.

DNA methylation status of nuclear-encoded mitochondrial genes underlies the tissue-dependent mitochondrial functions

Background

Mitochondria are semi-autonomous, semi-self-replicating organelles harboring their own DNA (mitochondrial DNA, mtDNA), and their dysregulation is involved in the development of various diseases. While mtDNA does not generally undergo epigenetic modifications, almost all mitochondrial proteins are encoded by nuclear DNA. However, the epigenetic regulation of nuclear-encoded mitochondrial genes (nuclear mt genes) has not been comprehensively analyzed.

Results

We analyzed the DNA methylation status of 899 nuclear mt genes in the liver, brain, and heart tissues of mouse, and identified 636 nuclear mt genes carrying tissue-dependent and differentially methylated regions (T-DMRs). These nuclar mt genes are involved in various mitochondrial functions and they also include genes related to human diseases. T-DMRs regulate the expression of nuclear mt genes. Nuclear mt genes with tissue-specific hypomethylated T-DMRs were characterized by enrichment of the target genes of specific transcription factors such as FOXA2 in the liver, and CEBPA and STAT1 in the brain.

Conclusions

A substantial proportion of nuclear mt genes contained T-DMRs, and the DNA methylation status of numerous T-DMRs should underlie tissue-dependent mitochondrial functions.

Molecular interactions between HNF4a, FOXA2 and GABP identified at regulatory DNA elements through ChIP-sequencing

Gene expression is regulated by combinations of transcription factors, which can be mapped to regulatory elements on a genome-wide scale using ChIP experiments. In a previous ChIP-chip study of USF1 and USF2 we found evidence also of binding of GABP, FOXA2 and HNF4a within the enriched regions. Here, we have applied ChIP-seq for these transcription factors and identified 3064 peaks of enrichment for GABP, 7266 for FOXA2 and 18783 for HNF4a. Distal elements with USF2 signal was frequently bound also by HNF4a and FOXA2. GABP peaks were found at transcription start sites, whereas 94% of FOXA2 and 90% of HNF4a peaks were located at other positions. We developed a method to accurately define TFBS within peaks, and found the predicted sites to have an elevated conservation level compared to peak centers; however the majority of bindings were not evolutionary conserved. An interaction between HNF4a and GABP was seen at TSS, with one-third of the HNF4a positive promoters being bound also by GABP, and this interaction was verified by co-immunoprecipitations.

Regulatory interdependence of myeloid transcription factors revealed by Matrix RNAi analysis

The knockdown of 78 transcription factors in differentiating human THP-1 cells using matrix RNAi reveals their interdependence

Background

With the move towards systems biology, we need sensitive and reliable ways to determine the relationships between transcription factors and their target genes. In this paper we analyze the regulatory relationships between 78 myeloid transcription factors and their coding genes by using the matrix RNAi system in which a set of transcription factor genes are individually knocked down and the resultant expression perturbation is quantified.

Results

Using small interfering RNAs we knocked down the 78 transcription factor genes in monocytic THP-1 cells and monitored the perturbation of the expression of the same 78 transcription factors and 13 other transcription factor genes as well as 5 non-transcription factor genes by quantitative real-time RT-PCR, thereby building a 78 × 96 matrix of perturbation and measurement. This approach identified 876 cases where knockdown of one transcription factor significantly affected the expression of another (from a potential 7,488 combinations). Our study also revealed cell-type-specific transcriptional regulatory networks in two different cell types.

Conclusions

By considering whether the targets of a given transcription factor are naturally up- or downregulated during phorbol 12-myristate 13-acetate-induced differentiation, we could classify these edges as pro-differentiative (229), anti-differentiative (76) or neither (571) using expression profiling data obtained in the FANTOM4 study. This classification analysis suggested that several factors could be involved in monocytic differentiation, while others such as MYB and the leukemogenic fusion MLL-MLLT3 could help to maintain the initial undifferentiated state by repressing the expression of pro-differentiative factors or maintaining expression of anti-differentiative factors.

Identification of an inter-transcription factor regulatory network in human hepatoma cells by Matrix RNAi

Transcriptional regulation by transcriptional regulatory factors (TRFs) of their target TRF genes is central to the control of gene expression. To study a static multi-tiered inter-TRF regulatory network in the human hepatoma cells, we have applied a Matrix RNAi approach in which siRNA knockdown and quantitative RT-PCR are used in combination on the same set of TRFs to determine their interdependencies. This approach focusing on several liver-enriched TRF families, each of which consists of structurally homologous members, revealed many significant regulatory relationships. These include the cross-talks between hepatocyte nuclear factors (HNFs) and the other TRF groups such as CCAAT/enhancer-binding proteins (CEBPs), retinoic acid receptors (RARs), retinoid receptors (RXRs) and RAR-related orphan receptors (RORs), which play key regulatory functions in human hepatocytes and liver. In addition, various multi-component regulatory motifs, which make up the complex inter-TRF regulatory network, were identified. A large part of the regulatory edges identified by the Matrix RNAi approach could be confirmed by chromatin immunoprecipitation. The resultant significant edges enabled us to depict the inter-TRF TRN forming an apparent regulatory hierarchy of (FOXA1, RXRA) → TCF1 → (HNF4A, ONECUT1) → (RORC, CEBPA) as the main streamline.