Reconstruction of monocyte transcriptional regulatory network accompanies monocytic functions in human fibroblasts

MNDA, a PYHIN factor involved in transcriptional regulation and apoptosis control in leukocytes

Inflammation control is critical during the innate immune response. Such response is triggered by the detection of molecules originating from pathogens or damaged host cells by pattern-recognition receptors (PRRs). PRRs subsequently initiate intra-cellular signalling through different pathways, resulting in i) the production of inflammatory cytokines, including type I interferon (IFN), and ii) the initiation of a cascade of events that promote both immediate host responses as well as adaptive immune responses. All human PYRIN and HIN-200 domains (PYHIN) protein family members were initially proposed to be PRRs, although this view has been challenged by reports that revealed their impact on other cellular mechanisms. Of relevance here, the human PYHIN factor myeloid nuclear differentiation antigen (MNDA) has recently been shown to directly control the transcription of genes encoding factors that regulate programmed cell death and inflammation. While MNDA is mainly found in the nucleus of leukocytes of both myeloid (neutrophils and monocytes) and lymphoid (B-cell) origin, its subcellular localization has been shown to be modulated in response to genotoxic agents that induce apoptosis and by bacterial constituents, mediators of inflammation. Prior studies have noted the importance of MNDA as a marker for certain forms of lymphoma, and as a clinical prognostic factor for hematopoietic diseases characterized by defective regulation of apoptosis. Abnormal expression of MNDA has also been associated with altered levels of cytokines and other inflammatory mediators. Refining our comprehension of the regulatory mechanisms governing the expression of MNDA and other PYHIN proteins, as well as enhancing our definition of their molecular functions, could significantly influence the management and treatment strategies of numerous human diseases. Here, we review the current state of knowledge regarding PYHIN proteins and their role in innate and adaptive immune responses. Emphasis will be placed on the regulation, function, and relevance of MNDA expression in the control of gene transcription and RNA stability during cell death and inflammation.

Structural mechanism of dsDNA recognition by the hMNDA HIN domain: New insights into the DNA-binding model of a PYHIN protein

The hematopoietic interferon-inducible nuclear (HIN) domain of the PYHIN family of proteins recognizes double-stranded DNA (dsDNA) through different dsDNA-binding modes. These modes apparently confer different roles upon these proteins in the regulation of innate immune responses, gene transcription, and apoptosis. Myeloid cell nuclear differentiation antigen (MNDA), a member of the human PYHIN family, binds DNA and regulates gene transcription in monocytes. However, the mechanism of DNA recognition and DNA-binding modes of human MNDA (hMNDA) remain unclear. Here, we determined the crystal structure of the hMNDA-HIN domain in complex with dsDNA at 2.4 Å resolution, and reveal that hMNDA-HIN binds to dsDNA in a sequence-independent manner. Structure and mutation studies indicated that hMNDA-HIN binds to dsDNA through a unique mode, involving two dsDNA-binding interfaces. Interface I exhibits an AIM2-like dsDNA-binding mode, and interface II has a previously unreported mode of dsDNA-binding. These results provide new insights into the DNA-binding modes of this PYHIN protein.

Quantitative proteomic analysis of human serum using tandem mass tags to predict cardiovascular risks in patients with psoriasis

Although biomarker candidates associated with psoriasis have been suggested, those for predicting the risk of cardiovascular disease (CVD) early in patients with psoriasis are lacking. We aimed to identify candidate biomarkers that can predict the occurrence of CVD in psoriasis patients. We pursued quantitative proteomic analysis of serum samples composed of three groups: psoriasis patients with and those without CVD risk factors, and healthy controls. Age/Sex-matched serum samples were selected and labeled with 16-plex tandem mass tag (TMT) and analyzed using liquid chromatography-mass spectrometry and subsequent verification with ELISA. Of the 184 proteins that showed statistical significance (P-value < 0.05) among the three groups according to TMT-based quantitative analysis, 98 proteins showed significant differences (> 2.0-fold) between the psoriasis groups with and without CVD risk factors. Verification by ELISA revealed that caldesmon (CALD1), myeloid cell nuclear differentiation antigen (MNDA), and zyxin (ZYX) levels were significantly increased in the psoriasis group with CVD risk factors. Further network analysis identified pathways including integrin signaling, which could be related to platelet aggregation, and actin cytoskeleton signaling. Three novel candidates (MNDA, ZYX, and CALD1) could be potential biomarkers for predicting CVD risks in psoriasis patients. We expect these biomarker candidates can be used to predict CVD risk in psoriasis patients in clinical settings although further studies including large validation are needed.

Prediction of transcription factors associated with DNA demethylation during human cellular development

DNA methylation of CpG dinucleotides is an important epigenetic modification involved in the regulation of mammalian gene expression, with each type of cell developing a specific methylation profile during its differentiation. Recently, it has been shown that a small subgroup of transcription factors (TFs) might promote DNA demethylation at their binding sites. We developed a bioinformatics pipeline to predict from genome-wide DNA methylation data TFs that promote DNA demethylation at their binding site. We applied the pipeline to International Human Epigenome Consortium methylome data and selected 393 candidate transcription factor binding motifs and associated 383 TFs that are likely associated with DNA demethylation. Validation of a subset of the candidate TFs using an in vitro assay suggested that 28 of 49 TFs from various TF families had DNA-demethylation-promoting activity; TF families, such as bHLH and ETS, contained both TFs with and without the activity. The identified TFs showed large demethylated/methylated CpG ratios and their demethylated CpGs showed significant bias toward hypermethylation in original cells. Furthermore, the identified TFs promoted demethylation of distinct sets of CpGs, with slight overlap of the targeted CpGs among TF family members, which was consistent with the results of a gene ontology (GO) term analysis of the identified TFs. Gene expression analysis of the identified TFs revealed that multiple TFs from various families are specifically expressed in human cells and tissues. Together, our results suggest that a large number of TFs from various TF families are associated with cell-type-specific DNA demethylation during human cellular development.

Development of p53 knockout U87MG cell line for unbiased drug delivery testing system using CRISPR-Cas9 and transcriptomic analysis

Antibody-drug conjugates offers many advantages as a drug delivery platform that allows for highly specific targeting of cell types and genes. Ideally, testing the efficacy of these systems requires two cell types to be different only in the gene targeted by the drug, with the rest of the cellular machinery unchanged, in order to minimize other potential differences from obscuring the effects of the drug. In this study, we created multiple variants of U87MG cells with targeted mutation in the TP53 gene using the CRISPR-Cas9 system, and determined that their major transcriptional differences stem from the loss of p53 function. Using the transcriptome data, we predicted which mutant clones would have less divergent phenotypes from the wild type and thereby serve as the best candidates to be used as drug delivery testing platforms. Further in vitro and in vivo assays of cell morphology, proliferation rate and target antigen-mediated uptake supported our predictions. Based on the combined analysis results, we successfully selected the best qualifying mutant clone. This study serves as proof-of-principle of the approach and paves the way for extending to additional cell types and target genes.

Emerging Role of PYHIN Proteins as Antiviral Restriction Factors

Innate immune sensors and restriction factors are cellular proteins that synergize to build an effective first line of defense against viral infections. Innate sensors are usually constitutively expressed and capable of detecting pathogen-associated molecular patterns (PAMPs) via specific pattern recognition receptors (PRRs) to stimulate the immune response. Restriction factors are frequently upregulated by interferons (IFNs) and may inhibit viral pathogens at essentially any stage of their replication cycle. Members of the Pyrin and hematopoietic interferon-inducible nuclear (HIN) domain (PYHIN) family have initially been recognized as important sensors of foreign nucleic acids and activators of the inflammasome and the IFN response. Accumulating evidence shows, however, that at least three of the four members of the human PYHIN family restrict viral pathogens independently of viral sensing and innate immune activation. In this review, we provide an overview on the role of human PYHIN proteins in the innate antiviral immune defense and on viral countermeasures.

Transcriptome analysis reveals inadequate spermatogenesis and immediate radical immune reactions during organ culture in vitro spermatogenesis

Cultivation of neonatal mouse testis tissue can induce spermatogenesis and produce fertile sperms. However, in vitro spermatogenesis mediated by the current organ culture method comes short in fully mimicking the in vivo counterpart, partly due to a lack of knowledge underlying molecular phenotypes of in vitro spermatogenesis. In this study, we investigated transcriptome of cultured testis tissues using microarray method. Principle component analysis of the transcriptome data revealed delay and/or arrest of spermatogenesis and immediate radical immune reactions in the cultured testis tissues. The delay/arrest of spermatogenesis occurred before and during early meiotic phase, resulting in inefficient progression of meiosis. The immune reaction, on the other hand, was drastic and overwhelming, in which TLR4-NF-kB signaling was speculated to be involved. Notably, treatment with TAK242, an inhibitor of TLR4-NF-kB signaling pathway, ameliorated the macrophage activation which otherwise would exacerbate the inflammation. Thus, the present study revealed for the first time at molecular level that the deficiency of germ cell differentiation and the immense immune reaction are major abnormalities in the cultured testis tissues.

Direct cell-fate conversion of somatic cells: Toward regenerative medicine and industries

Cells of multicellular organisms have diverse characteristics despite having the same genetic identity. The distinctive phenotype of each cell is determined by molecular mechanisms such as epigenetic changes that occur throughout the lifetime of an individual. Recently, technologies that enable modification of the fate of somatic cells have been developed, and the number of studies using these technologies has increased drastically in the last decade. Various cell types, including neuronal cells, cardiomyocytes, and hepatocytes, have been generated using these technologies. Although most direct reprogramming methods employ forced transduction of a defined sets of transcription factors to reprogram cells in a manner similar to induced pluripotent cell technology, many other strategies, such as methods utilizing chemical compounds and microRNAs to change the fate of somatic cells, have also been developed. In this review, we summarize transcription factor-based reprogramming and various other reprogramming methods. Additionally, we describe the various industrial applications of direct reprogramming technologies.

Interferon-Responsive Genes Are Targeted during the Establishment of Human Cytomegalovirus Latency

Human cytomegalovirus (HCMV) latency is an active process which remodels the latently infected cell to optimize latent carriage and reactivation. This is achieved, in part, through the expression of viral genes, including the G-protein-coupled receptor US28. Here, we use an unbiased proteomic screen to assess changes in host proteins induced by US28, revealing that interferon-inducible genes are downregulated by US28. We validate that major histocompatibility complex (MHC) class II and two pyrin and HIN domain (PYHIN) proteins, myeloid cell nuclear differentiation antigen (MNDA) and IFI16, are downregulated during experimental latency in primary human CD14+ monocytes. We find that IFI16 is targeted rapidly during the establishment of latency in a US28-dependent manner but only in undifferentiated myeloid cells, a natural site of latent carriage. Finally, by overexpressing IFI16, we show that IFI16 can activate the viral major immediate early promoter and immediate early gene expression during latency via NF-κB, a function which explains why downregulation of IFI16 during latency is advantageous for the virus.IMPORTANCE Human cytomegalovirus (HCMV) is a ubiquitous herpesvirus which infects 50 to 100% of humans worldwide. HCMV causes a lifelong subclinical infection in immunocompetent individuals but is a serious cause of mortality and morbidity in the immunocompromised and neonates. In particular, reactivation of HCMV in the transplant setting is a major cause of transplant failure and related disease. Therefore, a molecular understanding of HCMV latency and reactivation could provide insights into potential ways to target the latent viral reservoir in at-risk patient populations.

Preparation of Induced Pluripotent Stem Cells Using Human Peripheral Blood Monocytes

Since induced pluripotent stem (iPS) cells have been established, in recent years, clinical transplantation of cells differentiated from iPS cells derived from human skin fibroblasts is been in progress. On the contrary, monocytes have complete genome information without damage and gene recombination, they are contained in the peripheral blood by ∼3%–8% and differentiate into dendritic cells that are the type of control tower for immune cells. However, generation of monocyte-derived iPS cells has only been successful when special persistent Sendai virus vectors have been used. Therefore, in this study, as a preculture method for monocytes, a culture method for maintaining activity without using any cytokine was established, and using a commercially available vector without genetic toxicity without damaging the chromosome of the cell, iPS cells derived from monocytes were successfully produced. This cell has the ability to differentiate into three germ layers, and when compared with commercially available iPS cells, there was no significant difference between self-renewal and gene expression in the three germ layers. In future, we will compare the differentiation induction of monocyte-derived iPS cells with dendritic cells and investigate the production of dendritic cells that can cope with various antigens.

Mitochondrial biogenesis is transcriptionally repressed in lysosomal lipid storage diseases

Perturbations in mitochondrial function and homeostasis are pervasive in lysosomal storage diseases, but the underlying mechanisms remain unknown. Here, we report a transcriptional program that represses mitochondrial biogenesis and function in lysosomal storage diseases Niemann-Pick type C (NPC) and acid sphingomyelinase deficiency (ASM), in patient cells and mouse tissues. This mechanism is mediated by the transcription factors KLF2 and ETV1, which are both induced in NPC and ASM patient cells. Mitochondrial biogenesis and function defects in these cells are rescued by the silencing of KLF2 or ETV1. Increased ETV1 expression is regulated by KLF2, while the increase of KLF2 protein levels in NPC and ASM stems from impaired signaling downstream sphingosine-1-phosphate receptor 1 (S1PR1), which normally represses KLF2. In patient cells, S1PR1 is barely detectable at the plasma membrane and thus unable to repress KLF2. This manuscript provides a mechanistic pathway for the prevalent mitochondrial defects in lysosomal storage diseases.

Editorial note

This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

A screening system to identify transcription factors that induce binding site-directed DNA demethylation

Background

DNA methylation is a fundamental epigenetic modification that is involved in many biological systems such as differentiation and disease. We and others recently showed that some transcription factors (TFs) are involved in the site-specific determination of DNA demethylation in a binding site-directed manner, although the reports of such TFs are limited.

Results

Here, we develop a screening system to identify TFs that induce binding site-directed DNA methylation changes. The system involves the ectopic expression of target TFs in model cells followed by DNA methylome analysis and overrepresentation analysis of the corresponding TF binding motif at differentially methylated regions. It successfully identified binding site-directed demethylation of SPI1, which is known to promote DNA demethylation in a binding site-directed manner. We extended our screening system to 15 master TFs involved in cellular differentiation and identified eight novel binding site-directed DNA demethylation-inducing TFs (RUNX3, GATA2, CEBPB, MAFB, NR4A2, MYOD1, CEBPA, and TBX5). Gene ontology and tissue enrichment analysis revealed that these TFs demethylate genomic regions associated with corresponding biological roles. We also describe the characteristics of binding site-directed DNA demethylation induced by these TFs, including the targeting of highly methylated CpGs, local DNA demethylation, and the overlap of demethylated regions between TFs of the same family.

Conclusions

Our results show the usefulness of the developed screening system for the identification of TFs that induce DNA demethylation in a site-directed manner.

Electronic supplementary material

The online version of this article (10.1186/s13072-017-0169-6) contains supplementary material, which is available to authorized users.

Identifying novel transcription factors involved in the inflammatory response by using binding site motif scanning in genomic regions defined by histone acetylation

The innate immune response to pathogenic challenge is a complex, multi-staged process involving thousands of genes. While numerous transcription factors that act as master regulators of this response have been identified, the temporal complexity of gene expression changes in response to pathogen-associated molecular pattern receptor stimulation strongly suggest that additional layers of regulation remain to be uncovered. The evolved pathogen response program in mammalian innate immune cells is understood to reflect a compromise between the probability of clearing the infection and the extent of tissue damage and inflammatory sequelae it causes. Because of that, a key challenge to delineating the regulators that control the temporal inflammatory response is that an innate immune regulator that may confer a selective advantage in the wild may be dispensable in the lab setting. In order to better understand the complete transcriptional response of primary macrophages to the bacterial endotoxin lipopolysaccharide (LPS), we designed a method that integrates temporally resolved gene expression and chromatin-accessibility measurements from mouse macrophages. By correlating changes in transcription factor binding site motif enrichment scores, calculated within regions of accessible chromatin, with the average temporal expression profile of a gene cluster, we screened for transcriptional factors that regulate the cluster. We have validated our predictions of LPS-stimulated transcriptional regulators using ChIP-seq data for three transcription factors with experimentally confirmed functions in innate immunity. In addition, we predict a role in the macrophage LPS response for several novel transcription factors that have not previously been implicated in immune responses. This method is applicable to any experimental situation where temporal gene expression and chromatin-accessibility data are available.

RUNX1 regulates site specificity of DNA demethylation by recruitment of DNA demethylation machineries in hematopoietic cells

RUNX1 is an essential master transcription factor in hematopoietic development and plays important roles in immune functions. Although the gene regulatory mechanism of RUNX1 has been characterized extensively, the epigenetic role of RUNX1 remains unclear. Here, we demonstrate that RUNX1 contributes DNA demethylation in a binding site-directed manner in human hematopoietic cells. Overexpression analysis of RUNX1 showed the RUNX1-binding site-directed DNA demethylation. The RUNX1-mediated DNA demethylation was also observed in DNA replication-arrested cells, suggesting an involvement of active demethylation mechanism. Coimmunoprecipitation in hematopoietic cells showed physical interactions between RUNX1 and DNA demethylation machinery enzymes TET2, TET3, TDG, and GADD45. Further chromatin immunoprecipitation sequencing revealed colocalization of RUNX1 and TET2 in the same genomic regions, indicating recruitment of DNA demethylation machinery by RUNX1. Finally, methylome analysis revealed significant overrepresentation of RUNX1-binding sites at demethylated regions during hematopoietic development. Collectively, the present data provide evidence that RUNX1 contributes site specificity of DNA demethylation by recruitment of TET and other demethylation-related enzymes to its binding sites in hematopoietic cells.

Comprehensive transcriptome mining of the direct conversion of mesodermal cells

The direct reprogramming of somatic cells is a promising approach for regenerative medicine, especially in the production of mesoderm layer-derived cells. Meta-analysis studies provide precise insight into the undergoing processes and help increase the efficiency of reprogramming. Here, using 27 high-throughput expression data sets, we analyzed the direct reprogramming of mesodermal cells in humans and mice. Fibroblast-derived cells showed a common expression pattern of up- and down-regulated genes that were mainly involved in the suppression of the fibroblast-specific gene expression program, and may be used as markers of the initiation of reprogramming. Furthermore, we found a specific gene expression profile for each fibroblast-derived cell studied, and each gene set appeared to play specific functional roles in its cell type, suggesting their use as markers for their mature state. Furthermore, using data from protein-DNA interactions, we identified the main transcription factors (TFs) involved in the conversion process and ranked them based on their importance in their gene regulatory networks. In summary, our meta-analysis approach provides new insights on the direct conversion of mesodermal somatic cells, introduces a list of genes as markers for initiation and maturation, and identifies TFs for which manipulating their expression may increase the efficiency of direct conversion.

RUNX1 induces DNA replication independent active DNA demethylation at SPI1 regulatory regions

Background

SPI1 is an essential transcription factor (TF) for the hematopoietic lineage, in which its expression is tightly controlled through a −17-kb upstream regulatory region and a promoter region. Both regulatory regions are demethylated during hematopoietic development, although how the change of DNA methylation status is performed is still unknown.

Results

We found that the ectopic overexpression of RUNX1 (another key TF in hematopoiesis) in HEK-293T cells induces almost complete DNA demethylation at the −17-kb upstream regulatory region and partial but significant DNA demethylation at the proximal promoter region. This DNA demethylation occurred in mitomycin-C-treated nonproliferating cells at both regulatory regions, suggesting active DNA demethylation. Furthermore, ectopic RUNX1 expression induced significant endogenous SPI1 expression, although its expression level was much lower than that of natively SPI1-expressing monocyte cells.

Conclusions

These results suggest the novel role of RUNX1 as an inducer of DNA demethylation at the SPI1 regulatory regions, although the mechanism of RUNX1-induced DNA demethylation remains to be explored.

Loss of variation of state detected in soybean metabolic and human myelomonocytic leukaemia cell transcriptional networks under external stimuli

Soybean (Glycine max) is sensitive to flooding stress, and flood damage at the seedling stage is a barrier to growth. We constructed two mathematical models of the soybean metabolic network, a control model and a flooded model, from metabolic profiles in soybean plants. We simulated the metabolic profiles with perturbations before and after the flooding stimulus using the two models. We measured the variation of state that the system could maintain from a state-space description of the simulated profiles. The results showed a loss of variation of state during the flooding response in the soybean plants. Loss of variation of state was also observed in a human myelomonocytic leukaemia cell transcriptional network in response to a phorbol-ester stimulus. Thus, we detected a loss of variation of state under external stimuli in two biological systems, regardless of the regulation and stimulus types. Our results suggest that a loss of robustness may occur concurrently with the loss of variation of state in biological systems. We describe the possible applications of the quantity of variation of state in plant genetic engineering and cell biology. Finally, we present a hypothetical "external stimulus-induced information loss" model of biological systems.

Asymmetric Regulation of Peripheral Genes by Two Transcriptional Regulatory Networks

Transcriptional regulatory network (TRN) reconstitution and deconstruction occur simultaneously during reprogramming; however, it remains unclear how the starting and targeting TRNs regulate the induction and suppression of peripheral genes. Here we analyzed the regulation using direct cell reprogramming from human dermal fibroblasts to monocytes as the platform. We simultaneously deconstructed fibroblastic TRN and reconstituted monocytic TRN; monocytic and fibroblastic gene expression were analyzed in comparison with that of fibroblastic TRN deconstruction only or monocytic TRN reconstitution only. Global gene expression analysis showed cross-regulation of TRNs. Detailed analysis revealed that knocking down fibroblastic TRN positively affected half of the upregulated monocytic genes, indicating that intrinsic fibroblastic TRN interfered with the expression of induced genes. In contrast, reconstitution of monocytic TRN showed neutral effects on the majority of fibroblastic gene downregulation. This study provides an explicit example that demonstrates how two networks together regulate gene expression during cell reprogramming processes and contributes to the elaborate exploration of TRNs.

Distinct Patterns of Expression of Transcription Factors in Response to Interferonβ and Interferonλ1

After viral infection, type I and III interferons (IFNs) are coexpressed by respiratory epithelial cells (RECs) and activate the ISGF3 transcription factor (TF) complex to induce expression of a cell-specific set of interferon-stimulated genes (ISGs). Type I and III IFNs share a canonical signaling pathway, suggesting that they are redundant. Animal and in vitro models, however, have shown that they are not redundant. Because TFs dictate cellular phenotype and function, we hypothesized that focusing on TF-ISG will reveal critical combinatorial and nonredundant functions of type I or III IFN. We treated BEAS-2B human RECs with increasing doses of IFNβ or IFNλ1 and measured expression of TF-ISG. ISGs were expressed in a dose-dependent manner with a nonlinear jump at intermediate doses. At subsaturating combinations of IFNβ and IFNλ1, many ISGs were expressed in a pattern that we modeled with a cubic equation that mathematically defines this threshold effect. Uniquely, IFNβ alone induced early and transient IRF1 transcript and protein expression, while IFNλ1 alone induced IRF1 protein expression at low levels that were sustained through 24 h. In combination, saturating doses of these 2 IFNs together enhanced and sustained IRF1 expression. We conclude that the cubic model quantitates combinatorial effects of IFNβ and IFNλ1 and that IRF1 may mediate nonredundancy of type I or III IFN in RECs.

Paradigm shifts in genomics through the FANTOM projects

Big leaps in science happen when scientists from different backgrounds interact. In the past 15 years, the FANTOM Consortium has brought together scientists from different fields to analyze and interpret genomic data produced with novel technologies, including mouse full-length cDNAs and, more recently, expression profiling at single-nucleotide resolution by cap-analysis gene expression. The FANTOM Consortium has provided the most comprehensive mouse cDNA collection for functional studies and extensive maps of the human and mouse transcriptome comprising promoters, enhancers, as well as the network of their regulatory interactions. More importantly, serendipitous observations of the FANTOM dataset led us to realize that the mammalian genome is pervasively transcribed, even from retrotransposon elements, which were previously considered junk DNA. The majority of products from the mammalian genome are long non-coding RNAs (lncRNAs), including sense-antisense, intergenic, and enhancer RNAs. While the biological function has been elucidated for some lncRNAs, more than 98 % of them remain without a known function. We argue that large-scale studies are urgently needed to address the functional role of lncRNAs.

ARMADA: Using motif activity dynamics to infer gene regulatory networks from gene expression data

Analysis of gene expression data remains one of the most promising avenues toward reconstructing genome-wide gene regulatory networks. However, the large dimensionality of the problem prohibits the fitting of explicit dynamical models of gene regulatory networks, whereas machine learning methods for dimensionality reduction such as clustering or principal component analysis typically fail to provide mechanistic interpretations of the reduced descriptions. To address this, we recently developed a general methodology called motif activity response analysis (MARA) that, by modeling gene expression patterns in terms of the activities of concrete regulators, accomplishes dramatic dimensionality reduction while retaining mechanistic biological interpretations of its predictions (Balwierz, 2014). Here we extend MARA by presenting ARMADA, which models the activity dynamics of regulators across a time course, and infers the causal interactions between the regulators that drive the dynamics of their activities across time. We have implemented ARMADA as part of our ISMARA webserver, ismara.unibas.ch, allowing any researcher to automatically apply it to any gene expression time course. To illustrate the method, we apply ARMADA to a time course of human umbilical vein endothelial cells treated with TNF. Remarkably, ARMADA is able to reproduce the complex observed motif activity dynamics using a relatively small set of interactions between the key regulators in this system. In addition, we show that ARMADA successfully infers many of the key regulatory interactions known to drive this inflammatory response and discuss several novel interactions that ARMADA predicts. In combination with ISMARA, ARMADA provides a powerful approach to generating plausible hypotheses for the key interactions between regulators that control gene expression in any system for which time course measurements are available.

Repositioning monocyte TFRN into fibroblasts

Repositioning of cell-specific Transcription Factor Regulatory Networks (TFRNs) into another cell type is a pragmatic approach to elicit desired functions in impaired cells and also to investigate transcription regulations underlying cellular functions. However, a systematic method to identify a set of transcription factors (TFs) which can efficiently promote TFRN in a foreign cell type is lacking. Here, we demonstrate a perturbation-matrix method for the identification and repositioning of monocyte TFRNs in dermal fibroblasts.

The emerging role of human PYHIN proteins in innate immunity: implications for health and disease

The innate immune response depends on the ability of immune cells to detect pathogens through germline-encoded pattern recognition receptors (PRRs). Recently discovered PRRs include some members of the Pyrin and HIN domain (PYHIN) family, which are encoded on an interferon-inducible gene cluster located on chromosome 1q23. There are five human PYHIN proteins; Absent in melanoma 2 (AIM2), IFN-γ inducible protein 16 (IFI16), Myeloid cell nuclear differentiation antigen (MNDA), Pyrin and HIN domain family member 1 (PYHIN1) and the recently identified Pyrin domain only protein 3 (POP3). Early studies reported roles for these proteins in cell cycle control, tumour suppression and transcriptional regulation. AIM2 and IFI16 have now been shown to be immune sensors of non-self DNA, such as that produced by viruses in infected cells. AIM2 binds DNA to activate the inflammasome, while IFI16 detection of DNA can lead to the up-regulation of type I interferons or inflammasome activation. Recent studies have shown how IFI16 senses DNA viruses, and also how viruses evade detection by IFI16, while structural studies have greatly advanced our understanding of how AIM2 and IFI16 bind DNA to activate these immune responses. Furthermore, following the identification of POP3, interplay between members of this gene cluster has been established, with POP3 acting as a negative regulator of the AIM2 and IFI16 inflammasomes. In this review we discuss the current understanding of how PYHIN proteins function in innate immunity, their role in disease and the therapeutic possibilities that arise as a result.

A transient disruption of fibroblastic transcriptional regulatory network facilitates trans-differentiation

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.

ISMARA: automated modeling of genomic signals as a democracy of regulatory motifs

Accurate reconstruction of the regulatory networks that control gene expression is one of the key current challenges in molecular biology. Although gene expression and chromatin state dynamics are ultimately encoded by constellations of binding sites recognized by regulators such as transcriptions factors (TFs) and microRNAs (miRNAs), our understanding of this regulatory code and its context-dependent read-out remains very limited. Given that there are thousands of potential regulators in mammals, it is not practical to use direct experimentation to identify which of these play a key role for a particular system of interest. We developed a methodology that models gene expression or chromatin modifications in terms of genome-wide predictions of regulatory sites and completely automated it into a web-based tool called ISMARA (Integrated System for Motif Activity Response Analysis). Given only gene expression or chromatin state data across a set of samples as input, ISMARA identifies the key TFs and miRNAs driving expression/chromatin changes and makes detailed predictions regarding their regulatory roles. These include predicted activities of the regulators across the samples, their genome-wide targets, enriched gene categories among the targets, and direct interactions between the regulators. Applying ISMARA to data sets from well-studied systems, we show that it consistently identifies known key regulators ab initio. We also present a number of novel predictions including regulatory interactions in innate immunity, a master regulator of mucociliary differentiation, TFs consistently disregulated in cancer, and TFs that mediate specific chromatin modifications.

Computational modeling identifies key gene regulatory interactions underlying phenobarbital-mediated tumor promotion

Gene regulatory interactions underlying the early stages of non-genotoxic carcinogenesis are poorly understood. Here, we have identified key candidate regulators of phenobarbital (PB)-mediated mouse liver tumorigenesis, a well-characterized model of non-genotoxic carcinogenesis, by applying a new computational modeling approach to a comprehensive collection of in vivo gene expression studies. We have combined our previously developed motif activity response analysis (MARA), which models gene expression patterns in terms of computationally predicted transcription factor binding sites with singular value decomposition (SVD) of the inferred motif activities, to disentangle the roles that different transcriptional regulators play in specific biological pathways of tumor promotion. Furthermore, transgenic mouse models enabled us to identify which of these regulatory activities was downstream of constitutive androstane receptor and β-catenin signaling, both crucial components of PB-mediated liver tumorigenesis. We propose novel roles for E2F and ZFP161 in PB-mediated hepatocyte proliferation and suggest that PB-mediated suppression of ESR1 activity contributes to the development of a tumor-prone environment. Our study shows that combining MARA with SVD allows for automated identification of independent transcription regulatory programs within a complex in vivo tissue environment and provides novel mechanistic insights into PB-mediated hepatocarcinogenesis.

Genome-scale technology driven advances to research into normal and malignant haematopoiesis

Haematopoiesis or blood development has long served as a model system for adult stem cell biology. Moreover, when combined, the various cancers of the blood represent one of the commonest human malignancies. Large numbers of researchers have therefore dedicated their scientific careers to studying haematopoiesis for more than a century. Throughout this period, many new technologies have first been applied towards the study of blood cells, and the research fields of normal and malignant haematopoiesis have also been some of the earliest adopters of genome-scale technologies. This has resulted in significant new insights with implications ranging from basic biological mechanisms to patient diagnosis and prognosis and also produced lessons likely to be relevant for many other areas of biomedical research. This paper discusses the current state of play for a range of genome-scale applications within haemopoiesis research, including gene expression profiling, ChIP-sequencing, genomewide association analysis, and cancer genome sequencing. A concluding outlook section explores likely future areas of progress as well as potential technological and educational bottlenecks.

Myeloid nuclear differentiation antigen, neutrophil apoptosis and sepsis

Sepsis and septic shock are characterized by prolonged inflammation and delayed resolution, which are associated with suppression of neutrophil apoptosis. The role of the intrinsic apoptotic pathway and intracellular factors in regulation of neutrophil apoptosis remain incompletely understood. We previously reported that the nuclear factor MNDA (myeloid nuclear differentiation antigen) is fundamental to execution of the constitutive neutrophil death program. During neutrophil apoptosis MNDA is cleaved by caspases and relocated to the cytoplasm. However, when challenged with known mediators of sepsis, human neutrophils of healthy donors or neutrophils from patients with sepsis exhibited impaired MNDA relocation/cleavage parallel with myeloid cell leukemia-1 (MCL-1) accumulation and suppression of apoptosis. MNDA knockdown in a model cell line indicated that upon induction of apoptosis, MNDA promotes proteasomal degradation of MCL-1, thereby aggravating mitochondrial dysfunction. Thus, MNDA is central to a novel nucleus-mitochondrion circuit that promotes progression of apoptosis. Disruption of this circuit contributes to neutrophil longevity, thereby identifying MNDA as a potential therapeutic target in sepsis and other inflammatory pathologies.

Identification of ZNF395 as a novel modulator of adipogenesis

Adipogenesis is the process of cell differentiation by which mesenchymal stem cells (MSC) become adipocytes. Investigating the transcriptional regulatory process during adipogenesis may provide strategies to prevent obesity and other metabolic disorders. In recent years, numerous zinc finger proteins (ZFPs) have been implicated in regulating differentiation and cell fate determination. To investigate the regulatory role of ZFPs involved in adipogenesis, we performed genome-wide microarray expression profiling of an adipogenesis time series. Particularly focusing on the transiently responsive ZFPs, we identified and characterized the functional role of ZNF395 in adipogenesis. A systematic ablation of the ZNF395 transcript during adipogenesis revealed 40% reduction of adipocytes when compared to control. Furthermore, the number of adipocytes as well as the expression of key adipocyte markers were greatly induced when MSC were co-transduced with ZNF395 and PPARG2. To further elucidate the functional role of ZNF395 during adipogenesis, we attempted to trans-differentiate human dermal fibroblasts with PPARG2. The test remarkably revealed that ZNF395 in conjunction with PPARG2 greatly induced adipogenesis from dermal fibroblasts when compared to PPARG2 alone. These loss and gain of function experiments firmly establish that ZNF395 coordinate the transcriptional regulatory pathway with PPARG2, which may be necessary for the genesis of adipocytes.

Establishment of single-cell screening system for the rapid identification of transcriptional modulators involved in direct cell reprogramming

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.

[Involvement and paralysis of the ulnar nerve in elbow arthrosis]

The FANTOM5 project investigates transcription initiation activities in more than 1,000 human and mouse primary cells, cell lines and tissues using CAGE. Based on manual curation of sample information and development of an ontology for sample classification, we assemble the resulting data into a centralized data resource (http://fantom.gsc.riken.jp/5/). This resource contains web-based tools and data-access points for the research community to search and extract data related to samples, genes, promoter activities, transcription factors and enhancers across the FANTOM5 atlas.