Immediate early genes are considered to play important roles in dynamic gene regulatory networks following exposure to appropriate stimuli. One of the immediate early genes, early growth response gene 1 (EGR-1), has been implicated in differentiation of human monoblastoma cells along the monocytic commitment following treatment with phorbol ester. EGR-1 has been thought to work as a modifier of monopoiesis, but the precise function of EGR-1 in monocytic differentiation has not been fully elucidated.

Gene regulatory networks in living cells are controlled by the interaction of multiple cell type-specific transcription regulators with DNA binding sites in target genes. Interferon regulatory factor 8 (IRF8), also known as interferon consensus sequence binding protein (ICSBP), is a transcription factor expressed predominantly in myeloid and lymphoid cell lineages. To find the functional direct target genes of IRF8, the gene expression profiles of siRNA knockdown samples and genome-wide binding locations by ChIP-chip were analyzed in THP-1 myelomonocytic leukemia cells. Consequently, 84 genes were identified as functional direct targets. The ETS family transcription factor PU.1, also known as SPI1, binds to IRF8 and regulates basal transcription in macrophages. Using the same approach, we identified 53 direct target genes of PU.1; these overlapped with 19 IRF8 targets. These 19 genes included key molecules of IFN signaling such as OAS1 and IRF9, but excluded other IFN-related genes amongst the IRF8 functional direct target genes. We suggest that IRF8 and PU.1 can have both combined, and independent actions on different promoters in myeloid cells.

Transcriptional Regulatory Networks (TRNs) coordinate multiple transcription factors (TFs) in concert to maintain tissue homeostasis and cellular function. The re-establishment of target cell TRNs has been previously implicated in direct trans-differentiation studies where the newly introduced TFs switch on a set of key regulatory factors to induce de novo expression and function. However, the extent to which TRNs in starting cell types, such as dermal fibroblasts, protect cells from undergoing cellular reprogramming remains largely unexplored. In order to identify TFs specific to maintaining the fibroblast state, we performed systematic knockdown of 18 fibroblast-enriched TFs and analyzed differential mRNA expression against the same 18 genes, building a Matrix-RNAi. The resulting expression matrix revealed seven highly interconnected TFs. Interestingly, suppressing four out of seven TFs generated lipid droplets and induced PPARG and CEBPA expression in the presence of adipocyte-inducing medium only, while negative control knockdown cells maintained fibroblastic character in the same induction regime. Global gene expression analyses further revealed that the knockdown-induced adipocytes expressed genes associated with lipid metabolism and significantly suppressed fibroblast genes. Overall, this study reveals the critical role of the TRN in protecting cells against aberrant reprogramming, and demonstrates the vulnerability of donor cell's TRNs, offering a novel strategy to induce transgene-free trans-differentiations.

Hydroxymethylcytosines (hmC), one of several reported cytosine modifications, was recently found to be enriched in embryonic stem cells and neuronal cells, and thought to play an important role in regulating gene expression and cell specification. However, unlike methylcytosines (mC), the fate of hmC beyond DNA replication is not well understood. Here, to monitor the status of hmC during DNA replication, we prepared a stable episomal vector-based monitoring system called MoCEV in 293T cells. The MoCEV system containing fully hydroxymethylated-cytosine fragments revealed a significant modification towards mC after several rounds of DNA replication. Strikingly this modification was specifically observed at the CpG sites (71.9% of cytosines), whereas only 1.1% of modified cytosines were detected at the non-CpG sites. Since the unmodified MoCEV did not undergo any DNA methylation during cell division, the results strongly suggest that somatic cells undergo hmC to mC specifically at the CpG sites during cell division.

Histone modifications play an important role in gene regulation. Acetylation of histone 3 lysine 9 (H3K9ac) is generally associated with transcription initiation and unfolded chromatin, thereby positively influencing gene expression. Deep sequencing of the 5' ends of gene transcripts using DeepCAGE delivers detailed information about the architecture and expression level of gene promoters. The combination of H3K9ac ChIP-chip and DeepCAGE in a myeloid leukemia cell line (THP-1) allowed us to study the spatial distribution of H3K9ac around promoters using a novel clustering approach. The promoter classes were analyzed for association with relevant genomic sequence features.

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.

Mammals are composed of hundreds of different cell types with specialized functions. Each of these cellular phenotypes are controlled by different combinations of transcription factors. Using a human non islet cell insulinoma cell line (TC-YIK) which expresses insulin and the majority of known pancreatic beta cell specific genes as an example, we describe a general approach to identify key cell-type-specific transcription factors (TFs) and their direct and indirect targets. By ranking all human TFs by their level of enriched expression in TC-YIK relative to a broad collection of samples (FANTOM5), we confirmed known key regulators of pancreatic function and development. Systematic siRNA mediated perturbation of these TFs followed by qRT-PCR revealed their interconnections with NEUROD1 at the top of the regulation hierarchy and its depletion drastically reducing insulin levels. For 15 of the TF knock-downs (KD), we then used Cap Analysis of Gene Expression (CAGE) to identify thousands of their targets genome-wide (KD-CAGE). The data confirm NEUROD1 as a key positive regulator in the transcriptional regulatory network (TRN), and ISL1, and PROX1 as antagonists. As a complimentary approach we used ChIP-seq on four of these factors to identify NEUROD1, LMX1A, PAX6, and RFX6 binding sites in the human genome. Examining the overlap between genes perturbed in the KD-CAGE experiments and genes with a ChIP-seq peak within 50 kb of their promoter, we identified direct transcriptional targets of these TFs. Integration of KD-CAGE and ChIP-seq data shows that both NEUROD1 and LMX1A work as the main transcriptional activators. In the core TRN (i.e., TF-TF only), NEUROD1 directly transcriptionally activates the pancreatic TFs HSF4, INSM1, MLXIPL, MYT1, NKX6-3, ONECUT2, PAX4, PROX1, RFX6, ST18, DACH1, and SHOX2, while LMX1A directly transcriptionally activates DACH1, SHOX2, PAX6, and PDX1. Analysis of these complementary datasets suggests the need for caution in interpreting ChIP-seq datasets. (1) A large fraction of binding sites are at distal enhancer sites and cannot be directly associated to their targets, without chromatin conformation data. (2) Many peaks may be non-functional: even when there is a peak at a promoter, the expression of the gene may not be affected in the matching perturbation experiment.

Transcriptional regulatory networks (TRN) control the underlying mechanisms behind cellular functions and they are defined by a set of core transcription factors regulating cascades of peripheral genes. Here we report SPI1, CEBPA, MNDA and IRF8 as core transcription factors of monocyte TRN and demonstrate functional inductions of phagocytosis, inflammatory response and chemotaxis activities in human dermal fibroblasts. The Gene Ontology and KEGG pathway analyses also revealed notable representation of genes involved in immune response and endocytosis in fibroblasts. Moreover, monocyte TRN-inducers triggered multiple monocyte-specific genes based on the transcription factor motif response analysis and suggest that complex cellular TRNs are uniquely amenable to elicit cell-specific functions in unrelated cell types.

Combinatorial interactions of transcription modulators are critical to regulate cell-specific expression and to drive direct cell reprogramming (e.g. trans-differentiation). However, the identification of key transcription modulators from myriad of candidate genes is laborious and time consuming. To rapidly identify key regulatory factors involved in direct cell reprogramming, we established a multiplex single-cell screening system using a fibroblast-to-monocyte transition model. The system implements a single-cell 'shotgun-transduction' strategy followed by nested-single-cell-polymerase chain reaction (Nesc-PCR) gene expression analysis. To demonstrate this, we simultaneously transduced 18 monocyte-enriched transcription modulators in fibroblasts followed by selection of single cells expressing monocyte-specific CD14 and HLA-DR cell-surface markers from a heterogeneous population. Highly multiplex Nesc-PCR expression analysis revealed a variety of gene combinations with a significant enrichment of SPI1 (86/86) and a novel transcriptional modulator, HCLS1 (76/86), in the CD14(+)/HLA-DR(+) single cells. We could further demonstrate the synergistic role of HCLS1 in regulating monocyte-specific gene expressions and phagocytosis in dermal fibroblasts in the presence of SPI1. This study establishes a platform for a multiplex single-cell screening of combinatorial transcription modulators to drive any direct cell reprogramming.

Fungal community analyses in homes have been attracting attention because fungi are now generally considered to be allergens. Currently, these analyses are generally conducted using the culture method, although fungal communities in households often contain species that are difficult to culture. In contrast, next-generation sequencing (NGS) represents a comprehensive, labor- and time-saving approach that can facilitate species identification. However, the reliability of the NGS method has not been compared to that of the culture method. In this study, in an attempt to demonstrate the reliability of this application, we used the NGS method to target the internal transcribed spacer 1 (ITS1) in the fungal genome, conducted fungal community analyses for 18 house-dust samples and analyzed fungal community structures. The NGS method positively correlated with the culture method regarding the relative abundance of Aspergillus, Penicillium, Cladosporium and yeasts, which represent the major fungal components found in houses. Furthermore, several genera, such as Malassezia, could be sensitively detected. Our results imply that the reliability of the NGS method is comparable to that of the culture method and indicates that easily available databases may require modifications, including the removal of registrations that have not been sufficiently classified at the genus level.

Combinatorial interactions among transcription factors are critical to directing tissue-specific gene expression. To build a global atlas of these combinations, we have screened for physical interactions among the majority of human and mouse DNA-binding transcription factors (TFs). The complete networks contain 762 human and 877 mouse interactions. Analysis of the networks reveals that highly connected TFs are broadly expressed across tissues, and that roughly half of the measured interactions are conserved between mouse and human. The data highlight the importance of TF combinations for determining cell fate, and they lead to the identification of a SMAD3/FLI1 complex expressed during development of immunity. The availability of large TF combinatorial networks in both human and mouse will provide many opportunities to study gene regulation, tissue differentiation, and mammalian evolution.

Human DICER1 protein cleaves double-stranded RNA into small sizes, a crucial step in production of single-stranded RNAs which are mediating factors of cytoplasmic RNA interference. Here, we clearly demonstrate that human DICER1 protein localizes not only to the cytoplasm but also to the nucleoplasm. We also find that human DICER1 protein associates with the NUP153 protein, one component of the nuclear pore complex. This association is detected predominantly in the cytoplasm but is also clearly distinguishable at the nuclear periphery. Additional characterization of the NUP153-DICER1 association suggests NUP153 plays a crucial role in the nuclear localization of the DICER1 protein.

In the FANTOM5 project, transcription initiation events across the human and mouse genomes were mapped at a single base-pair resolution and their frequencies were monitored by CAGE (Cap Analysis of Gene Expression) coupled with single-molecule sequencing. Approximately three thousands of samples, consisting of a variety of primary cells, tissues, cell lines, and time series samples during cell activation and development, were subjected to a uniform pipeline of CAGE data production. The analysis pipeline started by measuring RNA extracts to assess their quality, and continued to CAGE library production by using a robotic or a manual workflow, single molecule sequencing, and computational processing to generate frequencies of transcription initiation. Resulting data represents the consequence of transcriptional regulation in each analyzed state of mammalian cells. Non-overlapping peaks over the CAGE profiles, approximately 200,000 and 150,000 peaks for the human and mouse genomes, were identified and annotated to provide precise location of known promoters as well as novel ones, and to quantify their activities.