Core transcriptional regulatory circuitry in human hepatocytes

Computational analysis of gene regulation in animal sleep deprivation

Sleep is an animal behavior shared by a wide range of species, suggesting that it must serve fundamental functions. However, the functions and molecular mechanisms underlying sleep are largely unknown. Through a meta-analysis of all available short-term sleep deprivation (SD) microarray data in mouse brain, we identified 91 key mouse SD-affected genes and two RBM3 isoforms showing opposite changes of expression during SD. Although most of the key SD-affected genes showed consistent changes of expression during SD across brain subregions despite their heterogeneous basal expression levels, we also identified the genes whose SD responses strongly depend upon the brain subregion. A gene regulatory network was also constructed for these genes showing that cAMP-responsive element (CRE) is one of the key cis-regulatory elements controlling SD-affected genes. We observed that SD during an animal's normal sleeping time could significantly disturb the circadian oscillation of clock genes. Surprisingly, synaptogenesis markers were significantly underexpressed in SD mice, differing from the previous findings in rat and fly. Comparing SD microarray data in mouse, rat, sparrow, and fly, we identified Egr and Nr4a gene families as conserved SD-affected genes, thus shedding new light on the origin of sleep in animals.

High-resolution profiling of the LEDGF/p75 chromatin interaction in the ENCODE region

Lens epithelium-derived growth factor/p75 (LEDGF/p75) is a transcriptional coactivator involved in stress response, autoimmune disease, cancer and HIV replication. A fusion between the nuclear pore protein NUP98 and LEDGF/p75 has been found in human acute and chronic myeloid leukemia and association of LEDGF/p75 with mixed-lineage leukemia (MLL)/menin is critical for leukemic transformation. During lentiviral replication, LEDGF/p75 tethers the pre-integration complex to the host chromatin resulting in a bias of integration into active transcription units (TUs). The consensus function of LEDGF/p75 is tethering of cargos to chromatin. In this regard, we determined the LEDGF/p75 chromatin binding profile. To this purpose, we used DamID technology and focused on the highly annotated ENCODE (Encyclopedia of DNA Elements) regions. LEDGF/p75 primarily binds downstream of the transcription start site of active TUs in agreement with the enrichment of HIV-1 integration sites at these locations. We show that LEDGF/p75 binding is not restricted to stress response elements in the genome, and correlation analysis with more than 200 genomic features revealed an association with active chromatin markers, such as H3 and H4 acetylation, H3K4 monomethylation and RNA polymerase II binding. Interestingly, some associations did not correlate with HIV-1 integration indicating that not all LEDGF/p75 complexes on the chromosome are amenable to HIV-1 integration.

Foxa2 is essential for mouse endometrial gland development and fertility

During embryonic development, Foxa2 is required for the formation of the node and notochord, and ablation of this gene results in defects in gastrulation, neural tube patterning, and gut morphogenesis. Foxa2 has been shown to be expressed specifically in the glandular epithelium of the murine uterus. To study the uterine function of Foxa2, this gene was conditionally ablated in the mouse uterus by crossing mice with floxed Foxa2 alleles, Foxa2(loxP/loxP), with the Pgr(cre) mouse model. Pgr(cre/+) Foxa2(loxP/loxP) mice showed significantly reduced fertility. Analysis of the uterus on Day 5.5 of pregnancy showed disrupted blastocyst implantation. Pgr(cre/+) Foxa2(loxP/loxP) mice also showed a severe impairment of the uterus to respond to the artificial induction of the decidual response. Morphological examination of the uteri of these mice showed a severe reduction in the number of endometrial glands. The loss of endometrial glands resulted in the reduction of leukemia inhibitory factor (Lif) expression. The lack of a decidual response could be partially rescued by an intrauterine injection of LIF before the initiation of the decidual response. This analysis demonstrates that Foxa2 regulates endometrial gland development and that mice with a loss of endometrial glands cannot support implantation in part due to the loss of LIF, which is a requisite for fertility in the mouse.

Networks inferred from biochemical data reveal profound differences in toll-like receptor and inflammatory signaling between normal and transformed hepatocytes

Systematic study of cell signaling networks increasingly involves high throughput proteomics, transcriptional profiling, and automated literature mining with the aim of assembling large scale interaction networks. In contrast, functional analysis of cell signaling usually focuses on a much smaller sets of proteins and eschews computation but focuses directly on cellular responses to environment and perturbation. We sought to combine these two traditions by collecting cell response measures on a reasonably large scale and then attempting to infer differences in network topology between two cell types. Human hepatocytes and hepatocellular carcinoma cell lines were exposed to inducers of inflammation, innate immunity, and proliferation in the presence and absence of small molecule drugs, and multiplex biochemical measurement was then performed on intra- and extracellular signaling molecules. We uncovered major differences between primary and transformed hepatocytes with respect to the engagement of toll-like receptor and NF-κB-dependent secretion of chemokines and cytokines that prime and attract immune cells. Overall, our results serve as a proof of principle for an approach to network analysis that is systematic, comparative, and biochemically focused. More specifically, our data support the hypothesis that hepatocellular carcinoma cells down-regulate normal inflammatory and immune responses to avoid immune editing.

Genomic analyses of hormone signaling and gene regulation

Many cellular signaling pathways ultimately control specific patterns of gene expression in the nucleus through a variety of signal-regulated transcription factors (TFs), including nuclear hormone receptors (NRs). The advent of genomic technologies for examining signal-regulated transcriptional responses and TF binding on a genomic scale has dramatically increased our understanding of the cellular programs that control hormonal signaling and gene regulation. Studies of TFs, especially NRs, using genomic approaches have revealed novel and unexpected features of hormone-regulated transcription, and a global view is beginning to emerge. In this review, we discuss the genomic methodologies that have been applied to the study of hormone-regulated gene expression, the results that have been obtained from using them, and the future prospects for these approaches. Given the wealth of information about hormone-dependent gene regulation by NRs, we have focused this review on the knowledge gained from genomic studies of their function.

A quantitative model of transcriptional regulation reveals the influence of binding location on expression

Understanding the mechanistic basis of transcriptional regulation has been a central focus of molecular biology since its inception. New high-throughput chromatin immunoprecipitation experiments have revealed that most regulatory proteins bind thousands of sites in mammalian genomes. However, the functional significance of these binding sites remains unclear. We present a quantitative model of transcriptional regulation that suggests the contribution of each binding site to tissue-specific gene expression depends strongly on its position relative to the transcription start site. For three cell types, we show that, by considering binding position, it is possible to predict relative expression levels between cell types with an accuracy approaching the level of agreement between different experimental platforms. Our model suggests that, for the transcription factors profiled in these cell types, a regulatory site's influence on expression falls off almost linearly with distance from the transcription start site in a 10 kilobase range. Binding to both evolutionarily conserved and non-conserved sequences contributes significantly to transcriptional regulation. Our approach also reveals the quantitative, tissue-specific role of individual proteins in activating or repressing transcription. These results suggest that regulator binding position plays a previously unappreciated role in influencing expression and blurs the classical distinction between proximal promoter and distal binding events.

Gene expression is controlled, in large part, by regulatory proteins called transcription factors that bind specific sites in the genome. A major focus of molecular biology has been understanding how these transcription factors interact with the cell's transcriptional machinery, the genome, and with each other to turn genes' expression on and off in various physiological contexts. Previous approaches for modeling transcriptional regulation have focused on the complex combinatorial interactions between groups of transcription factors at regulatory sites, or on the specific activating or repressive functions of individual proteins. In this work, we present a new modeling framework and demonstrate that an equally important, and previously overlooked, consideration in predicting the effect that a regulatory site has on gene expression is simply its location relative to the transcription start site of nearby genes. Our results show that, in general, the closer a binding event is to a gene's transcription start site, the more it influences expression. We also show that considering the particular proteins bound at a regulatory site helps predict the expression of nearby genes. However, considering the sequence conservation level of these sites does not lead to more accurate predictions.

Kruppel-like factor 4 (Klf4) prevents embryonic stem (ES) cell differentiation by regulating Nanog gene expression

Transcription factor Kruppel-like factor 4 (Klf4) is essential for somatic cell reprogramming. In addition, Klf4 seems to play a redundant role along with other Klf family proteins in embryonic stem (ES) cell self-renewal. However, how Klf4 regulates ES cell self-renewal and somatic cell reprogramming is still poorly understood. Here we report that Klf4 is required for both ES cell self-renewal and maintenance of pluripotency and that the expression of Klf4 prevents ES cell differentiation in response to withdrawal of leukemia inhibitory factor (LIF) or bone morphogenetic protein 4 (BMP4). In addition, Klf4 directly binds to the promoter region of Nanog and regulates its expression. Expression of Nanog prevents ES cell differentiation even when Klf4 gene expression is knocked down. On the other hand, knockdown of Nanog expression induces differentiation of ES cells that overexpress Klf4. Taken together, these results demonstrate that Klf4 functions upstream of Nanog in ES cell self-renewal and in preventing ES cell differentiation.

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.

Integrated approach for the identification of human hepatocyte nuclear factor 4alpha target genes using protein binding microarrays

Hepatocyte nuclear factor 4 alpha (HNF4alpha), a member of the nuclear receptor superfamily, is essential for liver function and is linked to several diseases including diabetes, hemophilia, atherosclerosis, and hepatitis. Although many DNA response elements and target genes have been identified for HNF4alpha, the complete repertoire of binding sites and target genes in the human genome is unknown. Here, we adapt protein binding microarrays (PBMs) to examine the DNA-binding characteristics of two HNF4alpha species (rat and human) and isoforms (HNF4alpha2 and HNF4alpha8) in a high-throughput fashion. We identified approximately 1400 new binding sequences and used this dataset to successfully train a Support Vector Machine (SVM) model that predicts an additional approximately 10,000 unique HNF4alpha-binding sequences; we also identify new rules for HNF4alpha DNA binding. We performed expression profiling of an HNF4alpha RNA interference knockdown in HepG2 cells and compared the results to a search of the promoters of all human genes with the PBM and SVM models, as well as published genome-wide location analysis. Using this integrated approach, we identified approximately 240 new direct HNF4alpha human target genes, including new functional categories of genes not typically associated with HNF4alpha, such as cell cycle, immune function, apoptosis, stress response, and other cancer-related genes.

CONCLUSION

We report the first use of PBMs with a full-length liver-enriched transcription factor and greatly expand the repertoire of HNF4alpha-binding sequences and target genes, thereby identifying new functions for HNF4alpha. We also establish a web-based tool, HNF4 Motif Finder, that can be used to identify potential HNF4alpha-binding sites in any sequence.

Identification of interacting transcription factors regulating tissue gene expression in human

BACKGROUND

Tissue gene expression is generally regulated by multiple transcription factors (TFs). A major first step toward understanding how tissues achieve their specificity is to identify, at the genome scale, interacting TFs regulating gene expression in different tissues. Despite previous discoveries, the mechanisms that control tissue gene expression are not fully understood.

RESULTS

We have integrated a function conservation approach, which is based on evolutionary conservation of biological function, and genes with highest expression level in human tissues to predict TF pairs controlling tissue gene expression. To this end, we have identified 2549 TF pairs associated with a certain tissue. To find interacting TFs controlling tissue gene expression in a broad spatial and temporal manner, we looked for TF pairs common to the same type of tissues and identified 379 such TF pairs, based on which TF-TF interaction networks were further built. We also found that tissue-specific TFs may play an important role in recruiting non-tissue-specific TFs to the TF-TF interaction network, offering the potential for coordinating and controlling tissue gene expression across a variety of conditions.

CONCLUSION

The findings from this study indicate that tissue gene expression is regulated by large sets of interacting TFs either on the same promoter of a gene or through TF-TF interaction networks.

Differential binding and co-binding pattern of FOXA1 and FOXA3 and their relation to H3K4me3 in HepG2 cells revealed by ChIP-seq

FOXA1 and FOXA3 binding patterns in HepG2 cells, together with their possible molecular interactions with FOXA2 and each other, are revealed by ChIP-seq.

Background

The forkhead box/winged helix family members FOXA1, FOXA2, and FOXA3 are of high importance in development and specification of the hepatic linage and the continued expression of liver-specific genes.

Results

Here, we present a genome-wide location analysis of FOXA1 and FOXA3 binding sites in HepG2 cells through chromatin immunoprecipitation with detection by sequencing (ChIP-seq) studies and compare these with our previous results on FOXA2. We found that these factors often bind close to each other in different combinations and consecutive immunoprecipitation of chromatin for one and then a second factor (ChIP-reChIP) shows that this occurs in the same cell and on the same DNA molecule, suggestive of molecular interactions. Using co-immunoprecipitation, we further show that FOXA2 interacts with both FOXA1 and FOXA3 in vivo, while FOXA1 and FOXA3 do not appear to interact. Additionally, we detected diverse patterns of trimethylation of lysine 4 on histone H3 (H3K4me3) at transcriptional start sites and directionality of this modification at FOXA binding sites. Using the sequence reads at polymorphic positions, we were able to predict allele specific binding for FOXA1, FOXA3, and H3K4me3. Finally, several SNPs associated with diseases and quantitative traits were located in the enriched regions.

Conclusions

We find that ChIP-seq can be used not only to create gene regulatory maps but also to predict molecular interactions and to inform on the mechanisms for common quantitative variation.

Evolution of transcriptional control in mammals

Changes in gene expression directed by transcriptional regulators can give rise to new phenotypes. While gene expression profiles can be maintained across large evolutionary distances, transcription factor-DNA interactions diverge rapidly. The application of new genome-wide methodologies has begun refining our global understanding of when and where mammalian transcription factors interact with DNA, thereby providing new insight into the mechanisms of transcriptional evolution. The interplay between cis and trans regulation of gene expression is an increasingly active area of investigation, and recent studies suggest that mutations in cis-regulatory DNA can explain many inter-species differences in gene expression.

Novel immortalized human fetal liver cell line, cBAL111, has the potential to differentiate into functional hepatocytes

Background

A clonal cell line that combines both stable hepatic function and proliferation capacity is desirable for in vitro applications that depend on hepatic function, such as pharmacological or toxicological assays and bioartificial liver systems. Here we describe the generation and characterization of a clonal human cell line for in vitro hepatocyte applications.

Results

Cell clones derived from human fetal liver cells were immortalized by over-expression of telomerase reverse transcriptase. The resulting cell line, cBAL111, displayed hepatic functionality similar to the parental cells prior to immortalization, and did not grow in soft agar. Cell line cBAL111 expressed markers of immature hepatocytes, like glutathione S transferase and cytokeratin 19, as well as progenitor cell marker CD146 and was negative for lidocaine elimination. On the other hand, the cBAL111 cells produced urea, albumin and cytokeratin 18 and eliminated galactose. In contrast to hepatic cell lines NKNT-3 and HepG2, all hepatic functions were expressed in cBAL111, although there was considerable variation in their levels compared with primary mature hepatocytes. When transplanted in the spleen of immunodeficient mice, cBAL111 engrafted into the liver and partly differentiated into hepatocytes showing expression of human albumin and carbamoylphosphate synthetase without signs of cell fusion.

Conclusion

This novel liver cell line has the potential to differentiate into mature hepatocytes to be used for in vitro hepatocyte applications.

The transcriptional foundation of pluripotency

A fundamental goal in biology is to understand the molecular basis of cell identity. Pluripotent embryonic stem (ES) cell identity is governed by a set of transcription factors centred on the triumvirate of Oct4, Sox2 and Nanog. These proteins often bind to closely localised genomic sites. Recent studies have identified additional transcriptional modulators that bind to chromatin near sites occupied by Oct4, Sox2 and Nanog. This suggests that the combinatorial control of gene transcription might be fundamental to the ES cell state. Here we discuss how these observations advance our understanding of the transcription factor network that controls pluripotent identity and highlight unresolved issues that arise from these studies.

A systems genetics approach implicates USF1, FADS3, and other causal candidate genes for familial combined hyperlipidemia

We hypothesized that a common SNP in the 3' untranslated region of the upstream transcription factor 1 (USF1), rs3737787, may affect lipid traits by influencing gene expression levels, and we investigated this possibility utilizing the Mexican population, which has a high predisposition to dyslipidemia. We first associated rs3737787 genotypes in Mexican Familial Combined Hyperlipidemia (FCHL) case/control fat biopsies, with global expression patterns. To identify sets of co-expressed genes co-regulated by similar factors such as transcription factors, genetic variants, or environmental effects, we utilized weighted gene co-expression network analysis (WGCNA). Through WGCNA in the Mexican FCHL fat biopsies we identified two significant Triglyceride (TG)-associated co-expression modules. One of these modules was also associated with FCHL, the other FCHL component traits, and rs3737787 genotypes. This USF1-regulated FCHL-associated (URFA) module was enriched for genes involved in lipid metabolic processes. Using systems genetics procedures we identified 18 causal candidate genes in the URFA module. The FCHL causal candidate gene fatty acid desaturase 3 (FADS3) was associated with TGs in a recent Caucasian genome-wide significant association study and we replicated this association in Mexican FCHL families. Based on a USF1-regulated FCHL-associated co-expression module and SNP rs3737787, we identify a set of causal candidate genes for FCHL-related traits. We then provide evidence from two independent datasets supporting FADS3 as a causal gene for FCHL and elevated TGs in Mexicans.

By integrating a genetic polymorphism with genome-wide gene expression levels, we were able to attribute function to a genetic polymorphism in the USF1 gene. The USF1 gene has previously been associated with a common dyslipidemia, FCHL. FCHL is characterized by elevated levels of total cholesterol, triglycerides, or both. We demonstrate that this genetic polymorphism in USF1 contributes to FCHL disease risk by modulating the expression of a group of genes functionally related to lipid metabolism, and that this modulation is mediated by USF1. One of the genes whose expression is modulated by USF1 is FADS3, which was also implicated in a recent genome-wide association study for lipid traits. We demonstrated that a genetic polymorphism from the FADS3 region, which was associated with triglycerides in a GWAS study of Caucasians, was also associated with triglycerides in Mexican FCHL families. Our analysis provides novel insight into the gene expression profile contributing to FCHL disease risk, and identifies FADS3 as a new gene for FCHL in Mexicans.

Mechanisms of liver development: concepts for understanding liver disorders and design of novel therapies

The study of liver development has significantly contributed to developmental concepts about morphogenesis and differentiation of other organs. Knowledge of the mechanisms that regulate hepatic epithelial cell differentiation has been essential in creating efficient cell culture protocols for programmed differentiation of stem cells to hepatocytes as well as developing cell transplantation therapies. Such knowledge also provides a basis for the understanding of human congenital diseases. Importantly, much of our understanding of organ development has arisen from analyses of patients with liver deficiencies. We review how the liver develops in the embryo and discuss the concepts that operate during this process. We focus on the mechanisms that control the differentiation and organization of the hepatocytes and cholangiocytes and refer to other reviews for the development of nonepithelial tissue in the liver. Much progress in the characterization of liver development has been the result of genetic studies of human diseases; gaining a better understanding of these mechanisms could lead to new therapeutic approaches for patients with liver disorders.

Characterization of glucocorticoid receptor and hepatocyte nuclear factor 4alpha (HNF4alpha) binding to the hnf4alpha gene in the liver

Hepatocyte nuclear factor 4alpha (HNF4alpha) plays a crucial role in hepatocyte differentiation, liver organogenesis and regulation of liver functions. In mouse liver, HNF4alpha is expressed from two promoters, P1 and P2, the latter being very weakly active and only in the embryo. Previously, using transfection assays we identified an enhancer upstream of P1 that mediates both HNF4alpha transactivation and glucocorticoid induction and showed that HNF4alpha1, originated from P1, represses activity of the P2 promoter, possibly through its indirect recruitment to the promoter. However, glucocorticoid receptor (GR) binding to the enhancer was not shown and HNF4alpha binding to P2, first reported in isolated human hepatocytes, was not confirmed in mouse liver. Here, to analyse glucocorticoid inducibility and auto-regulation of the hnf4alpha gene in the liver, we accurately mapped and quantitatively assessed GR and HNF4alpha binding to enhancer and HNF4alpha recruitment to the P2 promoter using chromatin immunoprecipitation (ChIP) and real-time PCR. We proved that GR binds to enhancer from embryonic day (E) 17.5 onward and HNF4alpha even earlier. We showed that HNF4alpha binds to P2 independently of the activation function (AF) 1 domain in adult liver. We mapped the binding region between -400 and -200 bp upstream of the transcription start site. Although Sp1 binds within this region in vitro, we did not find evidence of a role of this factor in HNF4alpha recruitment. Our results suggest that, in the liver, HNF4alpha expression may be induced by glucocorticoids around birth and positive auto-regulation of the gene may take place early in development. They support a model of P2 repression involving HNF4alpha recruitment to promoter, possibly through interaction with several promoter-bound factors.

Identification of an endogenous ligand bound to a native orphan nuclear receptor

Orphan nuclear receptors have been instrumental in identifying novel signaling pathways and therapeutic targets. However, identification of ligands for these receptors has often been based on random compound screens or other biased approaches. As a result, it remains unclear in many cases if the reported ligands are the true endogenous ligands, – i.e., the ligand that is bound to the receptor in an unperturbed in vivo setting. Technical limitations have limited our ability to identify ligands based on this rigorous definition. The orphan receptor hepatocyte nuclear factor 4 α (HNF4α) is a key regulator of many metabolic pathways and linked to several diseases including diabetes, atherosclerosis, hemophilia and cancer. Here we utilize an affinity isolation/mass-spectrometry (AIMS) approach to demonstrate that HNF4α is selectively occupied by linoleic acid (LA, C18:2ω6) in mammalian cells and in the liver of fed mice. Receptor occupancy is dramatically reduced in the fasted state and in a receptor carrying a mutation derived from patients with Maturity Onset Diabetes of the Young 1 (MODY1). Interestingly, however, ligand occupancy does not appear to have a significant effect on HNF4α transcriptional activity, as evidenced by genome-wide expression profiling in cells derived from human colon. We also use AIMS to show that LA binding is reversible in intact cells, indicating that HNF4α could be a viable drug target. This study establishes a general method to identify true endogenous ligands for nuclear receptors (and other lipid binding proteins), independent of transcriptional function, and to track in vivo receptor occupancy under physiologically relevant conditions.

Foxa2-dependent hepatic gene regulatory networks depend on physiological state

Bile acids are powerful detergents produced by the liver to aid in the absorption of dietary lipids. We recently reported a novel role for Foxa2 in bile acid metabolism. The winged helix transcription factor Foxa2 is required to prevent intrahepatic cholestasis and liver injury in mice fed a cholic acid-enriched diet. Here, we use functional genomics to study how Foxa2 regulates its targets in a cholic acid-dependent manner. We found that multiple signaling pathways essential for the hepatic response to acute liver injury are impaired in livers of Foxa2-deficient mice, suggesting that the deletion of Foxa2 in the hepatocyte affects the liver on a large scale. We also discovered distinct feed-forward regulatory loops controlling Foxa2-dependent targets in a cholic acid-dependent or -independent manner. We show that Foxa2 interacts with different transcription factors to achieve gene expression responses appropriate for each physiologic state.

Let-7 and miR-200 microRNAs: guardians against pluripotency and cancer progression

Micro (mi)RNAs are emerging as important regulators of cellular differentiation, their importance underscored by the fact that they are often dysregulated during carcinogenesis. Two evolutionary conserved families, let-7 and miR-200, regulate key differentiation processes during development. Loss of let-7 in cancer results in reverse embryogenesis and dedifferentiation, and miR-200 has been identified as a powerful regulator of epithelial-to-mesenchymal transition (EMT). Recent findings have connected let-7 with stem cell maintenance and point at a connection between EMT and stem cell formation. A part of tumor progression can be viewed as a continuum of progressive dedifferentiation (EMT) with a cell at the endpoint that has stem cell-like properties. I propose that steps of this process are driven by specific changes in the expression of let-7 and miR-200 family members.

Gene set-based module discovery in the breast cancer transcriptome

Background

Although microarray-based studies have revealed global view of gene expression in cancer cells, we still have little knowledge about regulatory mechanisms underlying the transcriptome. Several computational methods applied to yeast data have recently succeeded in identifying expression modules, which is defined as co-expressed gene sets under common regulatory mechanisms. However, such module discovery methods are not applied cancer transcriptome data.

Results

In order to decode oncogenic regulatory programs in cancer cells, we developed a novel module discovery method termed EEM by extending a previously reported module discovery method, and applied it to breast cancer expression data. Starting from seed gene sets prepared based on cis-regulatory elements, ChIP-chip data, and gene locus information, EEM identified 10 principal expression modules in breast cancer based on their expression coherence. Moreover, EEM depicted their activity profiles, which predict regulatory programs in each subtypes of breast tumors. For example, our analysis revealed that the expression module regulated by the Polycomb repressive complex 2 (PRC2) is downregulated in triple negative breast cancers, suggesting similarity of transcriptional programs between stem cells and aggressive breast cancer cells. We also found that the activity of the PRC2 expression module is negatively correlated to the expression of EZH2, a component of PRC2 which belongs to the E2F expression module. E2F-driven EZH2 overexpression may be responsible for the repression of the PRC2 expression modules in triple negative tumors. Furthermore, our network analysis predicts regulatory circuits in breast cancer cells.

Conclusion

These results demonstrate that the gene set-based module discovery approach is a powerful tool to decode regulatory programs in cancer cells.

Down-regulation of hepatic HNF4alpha gene expression during hyperinsulinemia via SREBPs

Mutations in the coding region of hepatocyte nuclear factor 4alpha (HNF4alpha), and its upstream promoter (P2) that drives expression in the pancreas, are known to lead to maturity-onset diabetes of the young 1 (MODY1). HNF4alpha also controls gluconeogenesis and lipid metabolism in the liver, where the proximal promoter (P1) predominates. However, very little is known about the role of hepatic HNF4alpha in diabetes. Here, we examine the expression of hepatic HNF4alpha in two diabetic mouse models, db/db mice (type 2, insulin resistant) and streptozotocin-treated mice (type 1, insulin deficient). We found that the level of HNF4alpha protein and mRNA was decreased in the liver of db/db mice but increased in streptozotocin-treated mice. Because insulin increases the activity of sterol regulatory element-binding proteins (SREBP)-1c and -2, we also examined the effect of SREBPs on hepatic HNF4alpha gene expression and found that, like insulin, ectopic expression of SREBPs decreases the level of hepatic HNF4alpha protein and mRNA both in vitro in primary hepatocytes and in vivo in the liver of C57BL/6 mice. Finally, we use gel shift, chromatin immunoprecipitation, small interfering RNA, and reporter gene analysis to show that SREBP2 binds the human HNF4alpha P1 promoter and negatively regulates its expression. These data indicate that hyperinsulinemia down-regulates HNF4alpha in the liver through the up-regulation of SREBPs, thereby establishing a link between these two critical transcription factor pathways that regulate lipid and glucose metabolism in the liver. These findings also provide new insights into diabetes-associated complications such as fatty liver disease.

Computational modeling of the hematopoietic erythroid-myeloid switch reveals insights into cooperativity, priming, and irreversibility

Hematopoietic stem cell lineage choices are decided by genetic networks that are turned ON/OFF in a switch-like manner. However, prior to lineage commitment, genes are primed at low expression levels. Understanding the underlying molecular circuitry in terms of how it governs both a primed state and, at the other extreme, a committed state is of relevance not only to hematopoiesis but also to developmental systems in general. We develop a computational model for the hematopoietic erythroid-myeloid lineage decision, which is determined by a genetic switch involving the genes PU.1 and GATA-1. Dynamical models based upon known interactions between these master genes, such as mutual antagonism and autoregulation, fail to make the system bistable, a desired feature for robust lineage determination. We therefore suggest a new mechanism involving a cofactor that is regulated as well as recruited by one of the master genes to bind to the antagonistic partner that is necessary for bistability and hence switch-like behavior. An interesting fallout from this architecture is that suppression of the cofactor through external means can lead to a loss of cooperativity, and hence to a primed state for PU.1 and GATA-1. The PU.1–GATA-1 switch also interacts with another mutually antagonistic pair, –FOG-1. The latter pair inherits the state of its upstream master genes and further reinforces the decision due to several feedback loops, thereby leading to irreversible commitment. The genetic switch, which handles the erythroid-myeloid lineage decision, is an example of a network that implements both a primed and a committed state by regulating cooperativity through recruitment of cofactors. Perturbing the feedback between the master regulators and downstream targets suggests potential reprogramming strategies. The approach points to a framework for lineage commitment studies in general and could aid the search for lineage-determining genes.

An important question in developmental biology is how different lineage choices are regulated at the genetic level. Robust lineage decisions are implemented by genetic switches, whereby one set of master genes are ON and another set are OFF, leading to a specific expression pattern of genes for a particular lineage. We develop a computational model to illustrate these principles as applied to the hematopoietic erythroid-myeloid lineage choice, where two master regulator genes, PU.1 and GATA-1, function as a genetic switch. The model, which is based upon known interactions, suggests missing interactions between the master genes, which we hypothesize, so as to reproduce the desired dynamics. Furthermore, there exist feedback interactions between the master genes and their downstream targets. When these are included in the model, the dynamics imply that the feedback is responsible for irreversible commitment. Our results suggest the search for missing interactions between the master genes in terms of a coregulated cofactor. The second important result of the model is that reprogramming irreversible cell fate may be possible by perturbing feedback regulation between the master genes and their downstream targets. Hence, dynamical modeling provides prediction of novel mechanisms and also strategies for reprogramming the fates of cells.

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.

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

Mammalian tissue- and/or time-specific transcription is primarily regulated in a combinatorial fashion through interactions between a specific set of transcriptional regulatory factors (TRFs) and their cognate cis-regulatory elements located in the regulatory regions. In exploring the DNA regions and TRFs involved in combinatorial transcriptional regulation, we noted that individual knockdown of a set of human liver-enriched TRFs such as HNF1A, HNF3A, HNF3B, HNF3G and HNF4A resulted in perturbation of the expression of several single TRF genes, such as HNF1A, HNF3G and CEBPA genes. We thus searched the potential binding sites for these five TRFs in the highly conserved genomic regions around these three TRF genes and found several putative combinatorial regulatory regions. Chromatin immunoprecipitation analysis revealed that almost all of the putative regulatory DNA regions were bound by the TRFs as well as two coactivators (CBP and p300). The strong transcription-enhancing activity of the putative combinatorial regulatory region located downstream of the CEBPA gene was confirmed. EMSA demonstrated specific bindings of these HNFs to the target DNA region. Finally, co-transfection reporter assays with various combinations of expression vectors for these HNF genes demonstrated the transcriptional activation of the CEBPA gene in a combinatorial manner by these TRFs.

Multistable decision switches for flexible control of epigenetic differentiation

It is now recognized that molecular circuits with positive feedback can induce two different gene expression states (bistability) under the very same cellular conditions. Whether, and how, cells make use of the coexistence of a larger number of stable states (multistability) is however largely unknown. Here, we first examine how autoregulation, a common attribute of genetic master regulators, facilitates multistability in two-component circuits. A systematic exploration of these modules' parameter space reveals two classes of molecular switches, involving transitions in bistable (progression switches) or multistable (decision switches) regimes. We demonstrate the potential of decision switches for multifaceted stimulus processing, including strength, duration, and flexible discrimination. These tasks enhance response specificity, help to store short-term memories of recent signaling events, stabilize transient gene expression, and enable stochastic fate commitment. The relevance of these circuits is further supported by biological data, because we find them in numerous developmental scenarios. Indeed, many of the presented information-processing features of decision switches could ultimately demonstrate a more flexible control of epigenetic differentiation.

Genome-scale validation of deep-sequencing libraries

Chromatin immunoprecipitation followed by high-throughput (HTP) sequencing (ChIP-seq) is a powerful tool to establish protein-DNA interactions genome-wide. The primary limitation of its broad application at present is the often-limited access to sequencers. Here we report a protocol, Mab-seq, that generates genome-scale quality evaluations for nucleic acid libraries intended for deep-sequencing. We show how commercially available genomic microarrays can be used to maximize the efficiency of library creation and quickly generate reliable preliminary data on a chromosomal scale in advance of deep sequencing. We also exploit this technique to compare enriched regions identified using microarrays with those identified by sequencing, demonstrating that they agree on a core set of clearly identified enriched regions, while characterizing the additional enriched regions identifiable using HTP sequencing.

Mapping key features of transcriptional regulatory circuitry in embryonic stem cells

The process by which a single fertilized egg develops into a human being with more than 200 cell types--each with a distinct gene expression pattern controlling its cellular state--is poorly understood. Knowledge of the transcriptional regulatory circuitry that establishes and maintains gene expression programs in mammalian cells is fundamental to understanding development and should provide the foundation for improved diagnosis and treatment of disease. Although it is not yet feasible to map the entirety of this circuitry in vertebrate cells, recent work in embryonic stem (ES) cells has demonstrated that core features of the circuitry can be discovered through studies involving selected regulators. Here, we highlight the fundamental insights that have emerged from studies that examined the role of transcription factors, chromatin regulators, signaling pathways, and noncoding RNAs in the regulatory circuitry of ES cells. Maps of regulatory circuitry and the insights that have emerged from these studies have improved our understanding of global gene expression and are facilitating efforts to reprogram cells for disease therapeutics and regenerative medicine.

Microarray analysis of Foxa2 mutant mouse embryos reveals novel gene expression and inductive roles for the gastrula organizer and its derivatives

Background

The Spemann/Mangold organizer is a transient tissue critical for patterning the gastrula stage vertebrate embryo and formation of the three germ layers. Despite its important role during development, there are still relatively few genes with specific expression in the organizer and its derivatives. Foxa2 is a forkhead transcription factor that is absolutely required for formation of the mammalian equivalent of the organizer, the node, the axial mesoderm and the definitive endoderm (DE). However, the targets of Foxa2 during embryogenesis, and the molecular impact of organizer loss on the gastrula embryo, have not been well defined.

Results

To identify genes specific to the Spemann/Mangold organizer, we performed a microarray-based screen that compared wild-type and Foxa2 mutant embryos at late gastrulation stage (E7.5). We could detect genes that were consistently down-regulated in replicate pools of mutant embryos versus wild-type, and these included a number of known node and DE markers. We selected 314 genes without previously published data at E7.5 and screened for expression by whole mount in situ hybridization. We identified 10 novel expression patterns in the node and 5 in the definitive endoderm. We also found significant reduction of markers expressed in secondary tissues that require interaction with the organizer and its derivatives, such as cardiac mesoderm, vasculature, primitive streak, and anterior neuroectoderm.

Conclusion

The genes identified in this screen represent novel Spemann/Mangold organizer genes as well as potential Foxa2 targets. Further investigation will be needed to define these genes as novel developmental regulatory factors involved in organizer formation and function. We have placed these genes in a Foxa2-dependent genetic regulatory network and we hypothesize how Foxa2 may regulate a molecular program of Spemann/Mangold organizer development. We have also shown how early loss of the organizer and its inductive properties in an otherwise normal embryo, impacts on the molecular profile of surrounding tissues.

Positional distribution of human transcription factor binding sites

We developed a method for estimating the positional distribution of transcription factor (TF) binding sites using ChIP-chip data, and applied it to recently published experiments on binding sites of nine TFs: OCT4, SOX2, NANOG, HNF1A, HNF4A, HNF6, FOXA2, USF1 and CREB1. The data were obtained from a genome-wide coverage of promoter regions from 8-kb upstream of the transcription start site (TSS) to 2-kb downstream. The number of target genes of each TF ranges from few hundred to several thousand. We found that for each of the nine TFs the estimated binding site distribution is closely approximated by a mixture of two components: a narrow peak, localized within 300-bp upstream of the TSS, and a distribution of almost uniform density within the tested region. Using Gene Ontology (GO) and Enrichment analysis, we were able to associate (for each of the TFs studied) the target genes of both types of binding with known biological processes. Most GO terms were enriched either among the proximal targets or among those with a uniform distribution of binding sites. For example, the three stemness-related TFs have several hundred target genes that belong to ‘development’ and ‘morphogenesis’ whose binding sites belong to the uniform distribution.

Analysis of regulatory network topology reveals functionally distinct classes of microRNAs

MicroRNAs (miRNAs) negatively regulate the expression of target genes at the post-transcriptional level. Little is known about the crosstalk between miRNAs and transcription factors (TFs). Here we provide data suggesting that the interaction patterns between TFs and miRNAs can influence the biological functions of miRNAs. From this global survey, we find that a regulated feedback loop, in which two TFs regulate each other and one miRNA regulates both of the factors, is the most significantly overrepresented network motif. Mathematical modeling shows that the miRNA in this motif stabilizes the feedback loop to resist environmental perturbation, providing one mechanism to explain the robustness of developmental programs that is contributed by miRNAs. Furthermore, on the basis of a network motif profile analysis, we demonstrate the existence of two classes of miRNAs with distinct network topological properties. The first class of miRNAs is regulated by a large number of TFs, whereas the second is regulated by only a few TFs. The differential expression level of the two classes of miRNAs in embryonic developmental stages versus adult tissues suggests that the two classes may have fundamentally different biological functions. Our results demonstrate that the TFs and miRNAs extensively interact with each other and the biological functions of miRNAs may be wired in the regulatory network topology.

Toward a systems-level understanding of developmental regulatory networks

Developmental regulatory networks constitute all the interconnections among molecular components that guide embryonic development. Developmental transcriptional regulatory networks (TRNs) are circuits of transcription factors and cis-acting DNA elements that control expression of downstream regulatory and effector genes. Developmental networks comprise functional subnetworks that are deployed sequentially in requisite spatiotemporal patterns. Here, we discuss integrative genomics approaches for elucidating TRNs, with an emphasis on those involved in Drosophila mesoderm development and mammalian embryonic stem cell maintenance and differentiation. As examples of regulatory subnetworks, we consider the transcriptional and signaling regulation of genes that interact to control cell morphology and migration. Finally, we describe integrative experimental and computational strategies for defining the entirety of molecular interactions underlying developmental regulatory networks.

Genome-wide occupancy of SREBP1 and its partners NFY and SP1 reveals novel functional roles and combinatorial regulation of distinct classes of genes

The sterol regulatory element-binding protein (SREBP) family member SREBP1 is a critical transcriptional regulator of cholesterol and fatty acid metabolism and has been implicated in insulin resistance, diabetes, and other diet-related diseases. We globally identified the promoters occupied by SREBP1 and its binding partners NFY and SP1 in a human hepatocyte cell line using chromatin immunoprecipitation combined with genome tiling arrays (ChIP-chip). We find that SREBP1 occupies the promoters of 1,141 target genes involved in diverse biological pathways, including novel targets with roles in lipid metabolism and insulin signaling. We also identify a conserved SREBP1 DNA-binding motif in SREBP1 target promoters, and we demonstrate that many SREBP1 target genes are transcriptionally activated by treatment with insulin and glucose using gene expression microarrays. Finally, we show that SREBP1 cooperates extensively with NFY and SP1 throughout the genome and that unique combinations of these factors target distinct functional pathways. Our results provide insight into the regulatory circuitry in which SREBP1 and its network partners coordinate a complex transcriptional response in the liver with cues from the diet.

Transcription factors (TFs) are DNA-binding proteins that regulate the transcription of their target genes. TFs typically bind in proximity to the transcription start sites of their target genes in a region called the promoter. SREBP1 is a TF that increases the transcription of numerous genes involved in cholesterol and fat metabolism and has been linked to diet-related diseases such as insulin resistance and type 2 diabetes. Using microarray technology, we identified all of the promoters in the human genome that are bound by SREBP1 and two associated TFs called NFY and SP1 in a human liver cell line. Our findings greatly expand the number of genes and biological pathways that may be regulated by SREBP1 and reveal that different combinations of SREBP1 and its partners preferentially target genes involved in different pathways. Thus, in contrast to traditional studies that focus on individual genes, we have used a genomics approach to provide a novel global view of the regulatory circuitry in which SREBP1 and its partners coordinate a transcriptional response in the liver with cues from the diet.

Variation in homeodomain DNA binding revealed by high-resolution analysis of sequence preferences

Most homeodomains are unique within a genome, yet many are highly conserved across vast evolutionary distances, implying strong selection on their precise DNA-binding specificities. We determined the binding preferences of the majority (168) of mouse homeodomains to all possible 8-base sequences, revealing rich and complex patterns of sequence specificity and showing that there are at least 65 distinct homeodomain DNA-binding activities. We developed a computational system that successfully predicts binding sites for homeodomain proteins as distant from mouse as Drosophila and C. elegans, and we infer full 8-mer binding profiles for the majority of known animal homeodomains. Our results provide an unprecedented level of resolution in the analysis of this simple domain structure and suggest that variation in sequence recognition may be a factor in its functional diversity and evolutionary success.

Global analysis of in vivo Foxa2-binding sites in mouse adult liver using massively parallel sequencing

Foxa2 (HNF3β) is a one of three, closely related transcription factors that are critical to the development and function of the mouse liver. We have used chromatin immunoprecipitation and massively parallel Illumina 1G sequencing (ChIP–Seq) to create a genome-wide profile of in vivo Foxa2-binding sites in the adult liver. More than 65% of the ∼11.5 k genomic sites associated with Foxa2 binding, mapped to extended gene regions of annotated genes, while more than 30% of intragenic sites were located within first introns. 20.5% of all sites were further than 50 kb from any annotated gene, suggesting an association with novel gene regions. QPCR analysis demonstrated a strong positive correlation between peak height and fold enrichment for Foxa2-binding sites. We measured the relationship between Foxa2 and liver gene expression by overlapping Foxa2-binding sites with a SAGE transcriptome profile, and found that 43.5% of genes expressed in the liver were also associated with Foxa2 binding. We also identified potential Foxa2-interacting transcription factors whose motifs were enriched near Foxa2-binding sites. Our comprehensive results for in vivo Foxa2-binding sites in the mouse liver will contribute to resolving transcriptional regulatory networks that are important for adult liver function.

Meta-analysis approach identifies candidate genes and associated molecular networks for type-2 diabetes mellitus

Background

Multiple functional genomics data for complex human diseases have been published and made available by researchers worldwide. The main goal of these studies is the detailed analysis of a particular aspect of the disease. Complementary, meta-analysis approaches try to extract supersets of disease genes and interaction networks by integrating and combining these individual studies using statistical approaches.

Results

Here we report on a meta-analysis approach that integrates data of heterogeneous origin in the domain of type-2 diabetes mellitus (T2DM). Different data sources such as DNA microarrays and, complementing, qualitative data covering several human and mouse tissues are integrated and analyzed with a Bootstrap scoring approach in order to extract disease relevance of the genes. The purpose of the meta-analysis is two-fold: on the one hand it identifies a group of genes with overall disease relevance indicating common, tissue-independent processes related to the disease; on the other hand it identifies genes showing specific alterations with respect to a single study. Using a random sampling approach we computed a core set of 213 T2DM genes across multiple tissues in human and mouse, including well-known genes such as Pdk4, Adipoq, Scd, Pik3r1, Socs2 that monitor important hallmarks of T2DM, for example the strong relationship between obesity and insulin resistance, as well as a large fraction (128) of yet barely characterized novel candidate genes. Furthermore, we explored functional information and identified cellular networks associated with this core set of genes such as pathway information, protein-protein interactions and gene regulatory networks. Additionally, we set up a web interface in order to allow users to screen T2DM relevance for any – yet non-associated – gene.

Conclusion

In our paper we have identified a core set of 213 T2DM candidate genes by a meta-analysis of existing data sources. We have explored the relation of these genes to disease relevant information and – using enrichment analysis – we have identified biological networks on different layers of cellular information such as signaling and metabolic pathways, gene regulatory networks and protein-protein interactions. The web interface is accessible via .

Cis-regulatory modules in the mammalian liver: composition depends on strength of Foxa2 consensus site

Foxa2 is a critical transcription factor that controls liver development and plays an important role in hepatic gluconeogensis in adult mice. Here, we use genome-wide location analysis for Foxa2 to identify its targets in the adult liver. We then show by computational analyses that Foxa2 containing cis-regulatory modules are not constructed from a random assortment of binding sites for other transcription factors expressed in the liver, but rather that their composition depends on the strength of the Foxa2 consensus site present. Genes containing a cis-regulatory module with a medium or weak Foxa2 consensus site are much more liver-specific than the genes with a strong consensus site. We not only provide a better understanding of the mechanisms of Foxa2 regulation but also introduce a novel method for identification of different cis-regulatory modules involving a single factor.

Associating transcription factor-binding site motifs with target GO terms and target genes

The roles and target genes of many transcription factors (TFs) are still unknown. To predict the roles of TFs, we present a computational method for associating Gene Ontology (GO) terms with TF-binding motifs. The method works by ranking all genes as potential targets of the TF, and reporting GO terms that are significantly associated with highly ranked genes. We also present an approach, whereby these predicted GO terms can be used to improve predictions of TF target genes. This uses a novel gene-scoring function that reflects the insight that genes annotated with GO terms predicted to be associated with the TF are more likely to be its targets. We construct validation sets of GO terms highly associated with known targets of various yeast and human TF. On the yeast reference sets, our prediction method identifies at least one correct GO term for 73% of the TF, 49% of the correct GO terms are predicted and almost one-third of the predicted GO terms are correct. Results on human reference sets are similarly encouraging. Validation of our target gene prediction method shows that its accuracy exceeds that of simple motif scanning.

The diabetic phenotype in HNF4A mutation carriers is moderated by the expression of HNF4A isoforms from the P1 promoter during fetal development

OBJECTIVE

Mutations in the alternatively spliced HNF4A gene cause maturity-onset diabetes of the young (MODY). We characterized the spatial and developmental expression patterns of HNF4A transcripts in human tissues and investigated their role as potential moderators of the MODY phenotype.

RESEARCH DESIGN AND METHODS

We measured the expression of HNF4A isoforms in human adult tissues and gestationally staged fetal pancreas by isoform-specific real-time PCR. The correlation between mutation position and age of diagnosis or age-related penetrance was assessed in a cohort of 190 patients with HNF4A mutations.

RESULTS

HNF4A was expressed exclusively from the P2 promoter in adult pancreas, but from 9 weeks until at least 26 weeks after conception, up to 23% of expression in fetal pancreas was of P1 origin. HNF4A4-6 transcripts were not detected in any tissue. In whole pancreas, HNF4A9 expression was greater than in islets isolated from the endocrine pancreas (relative level 22 vs. 7%). Patients with mutations in exons 9 and 10 (absent from HNF4A3, HNF4A6, and HNF4A9 isoforms) developed diabetes later than those with mutations in exons 2-8, where all isoforms were affected (40 vs. 24 years; P = 0.029). Exon 9/10 mutations were also associated with a reduced age-related penetrance (53 vs. 10% without diabetes at age 55 years; P < 0.00001).

CONCLUSIONS

We conclude that isoforms derived from the HNF4A P1 promoter are expressed in human fetal, but not adult, pancreas, and that their presence during pancreatic development may moderate the diabetic phenotype in individuals with mutations in the HNF4A gene.

Liver-specific hepatocyte nuclear factor-4alpha deficiency: greater impact on gene expression in male than in female mouse liver

Hepatocyte nuclear factor (HNF)-4alpha is a liver-enriched transcription factor that regulates numerous liver-expressed genes including several sex-specific cytochrome P450 genes. Presently, a liver-specific HNF4alpha-deficient mouse model was used to characterize the impact of liver HNF4alpha deficiency on a global scale using 41,174 feature microarrays. A total of 4994 HNF4alpha-dependent genes were identified, of which about 1000 fewer genes responded to the loss of HNF4alpha in female liver as compared with male liver. Sex differences in the impact of liver HNF4alpha deficiency were even more dramatic when genes showing sex-specific expression were examined. Thus, 372 of the 646 sex-specific genes characterized by a dependence on HNF4alpha responded to the loss of HNF4alpha in males only, as compared with only 61 genes that responded in females only. Moreover, in male liver, 78% of 508 male-specific genes were down-regulated and 42% of 356 female-specific genes were up-regulated in response to the loss of HNF4alpha, with sex specificity lost for 90% of sex-specific genes. This response to HNF4alpha deficiency is similar to the response of male mice deficient in the GH-activated transcription factor signal transducer and activator of transcription 5b (STAT5b), where 90% of male-specific genes were down-regulated and 61% of female-specific genes were up-regulated, suggesting these two factors cooperatively regulate liver sex specificity by mechanisms that are primarily active in males. Finally, 203 of 648 genes previously shown to bind HNF4alpha near the transcription start site in mouse hepatocytes were affected by HNF4alpha deficiency in mouse liver, with the HNF4alpha-bound gene set showing a 5-fold enrichment for genes positively regulated by HNF4alpha. Thus, a substantial fraction of the HNF4alpha-dependent genes reported here are likely to be direct targets of HNF4alpha.

Tcf3 is an integral component of the core regulatory circuitry of embryonic stem cells

Embryonic stem (ES) cells have a unique regulatory circuitry, largely controlled by the transcription factors Oct4, Sox2, and Nanog, which generates a gene expression program necessary for pluripotency and self-renewal. How external signals connect to this regulatory circuitry to influence ES cell fate is not known. We report here that a terminal component of the canonical Wnt pathway in ES cells, the transcription factor T-cell factor-3 (Tcf3), co-occupies promoters throughout the genome in association with the pluripotency regulators Oct4 and Nanog. Thus, Tcf3 is an integral component of the core regulatory circuitry of ES cells, which includes an autoregulatory loop involving the pluripotency regulators. Both Tcf3 depletion and Wnt pathway activation cause increased expression of Oct4, Nanog, and other pluripotency factors and produce ES cells that are refractory to differentiation. Our results suggest that the Wnt pathway, through Tcf3, brings developmental signals directly to the core regulatory circuitry of ES cells to influence the balance between pluripotency and differentiation.

Stem cells, the molecular circuitry of pluripotency and nuclear reprogramming

Reprogramming of somatic cells to a pluripotent embryonic stem cell-like state has been achieved by nuclear transplantation of a somatic nucleus into an enucleated egg and most recently by introducing defined transcription factors into somatic cells. Nuclear reprogramming is of great medical interest, as it has the potential to generate a source of patient-specific cells. Here, we review strategies to reprogram somatic cells to a pluripotent embryonic state and discuss our understanding of the molecular mechanisms of reprogramming based on recent insights into the regulatory circuitry of the pluripotent state.

Genome wide prediction of HNF4alpha functional binding sites by the use of local and global sequence context

An application of machine learning algorithms enables prediction of the functional context of transcription factor binding sites in the human genome.

We report an application of machine learning algorithms that enables prediction of the functional context of transcription factor binding sites in the human genome. We demonstrate that our method allowed de novo identification of hepatic nuclear factor (HNF)4α binding sites and significantly improved an overall recognition of faithful HNF4α targets. When applied to published findings, an unprecedented high number of false positives were identified. The technique can be applied to any transcription factor.

Cell cycle genes are the evolutionarily conserved targets of the E2F4 transcription factor

Maintaining quiescent cells in G0 phase is achieved in part through the multiprotein subunit complex known as DREAM, and in human cell lines the transcription factor E2F4 directs this complex to its cell cycle targets. We found that E2F4 binds a highly overlapping set of human genes among three diverse primary tissues and an asynchronous cell line, which suggests that tissue-specific binding partners and chromatin structure have minimal influence on E2F4 targeting. To investigate the conservation of these transcription factor binding events, we identified the mouse genes bound by E2f4 in seven primary mouse tissues and a cell line. E2f4 bound a set of mouse genes that was common among mouse tissues, but largely distinct from the genes bound in human. The evolutionarily conserved set of E2F4 bound genes is highly enriched for functionally relevant regulatory interactions important for maintaining cellular quiescence. In contrast, we found minimal mRNA expression perturbations in this core set of E2f4 bound genes in the liver, kidney, and testes of E2f4 null mice. Thus, the regulatory mechanisms maintaining quiescence are robust even to complete loss of conserved transcription factor binding events.

OSCAR: one-class SVM for accurate recognition of cis-elements

MOTIVATION

Traditional methods to identify potential binding sites of known transcription factors still suffer from large number of false predictions. They mostly use sequence information in a position-specific manner and neglect other types of information hidden in the proximal promoter regions. Recent biological and computational researches, however, suggest that there exist not only locational preferences of binding, but also correlations between transcription factors.

RESULTS

In this article, we propose a novel approach, OSCAR, which utilizes one-class SVM algorithms, and incorporates multiple factors to aid the recognition of transcription factor binding sites. Using both synthetic and real data, we find that our method outperforms existing algorithms, especially in the high sensitivity region. The performance of our method can be further improved by taking into account locational preference of binding events. By testing on experimentally-verified binding sites of GATA and HNF transcription factor families, we show that our algorithm can infer the true co-occurring motif pairs accurately, and by considering the co-occurrences of correlated motifs, we not only filter out false predictions, but also increase the sensitivity.

AVAILABILITY

An online server based on OSCAR is available at http://bioinfo.au.tsinghua.edu.cn/oscar.

Nuclear receptor hepatocyte nuclear factor 4alpha1 competes with oncoprotein c-Myc for control of the p21/WAF1 promoter

The dichotomy between cellular differentiation and proliferation is a fundamental aspect of both normal development and tumor progression; however, the molecular basis of this opposition is not well understood. To address this issue, we investigated the mechanism by which the nuclear receptor hepatocyte nuclear factor 4alpha1 (HNF4alpha1) regulates the expression of the human cyclin-dependent kinase inhibitor gene p21/WAF1 (CDKN1A). We found that HNF4alpha1, a transcription factor that plays a central role in differentiation in the liver, pancreas, and intestine, activates the expression of p21 primarily by interacting with promoter-bound Sp1 at both the proximal promoter region and at newly identified sites in a distal region (-2.4 kb). Although HNF4alpha1 also binds two additional regions containing putative HNF4alpha binding sites, HNF4alpha1 mutants deficient in DNA binding activate the p21 promoter to the same extent as wild-type HNF4alpha1, indicating that direct DNA binding by HNF4alpha1 is not necessary for p21 activation. We also observed an in vitro and in vivo interaction between HNF4alpha1 and c-Myc as well as a competition between these two transcription factors for interaction with promoter-bound Sp1 and regulation of p21. Finally, we show that c-Myc competes with HNF4alpha1 for control of apolipoprotein C3 (APOC3), a gene associated with the differentiated hepatic phenotype. These results suggest a general model by which a differentiation factor (HNF4alpha1) and a proliferation factor (c-Myc) may compete for control of genes involved in cell proliferation and differentiation.