Past Roadblocks and New Opportunities in Transcription Factor Network Mapping

Computational Reconstruction of the Transcription Factor Regulatory Network Induced by Auxin in Arabidopsis thaliana L

In plant hormone signaling, transcription factor regulatory networks (TFRNs), which link the master transcription factors to the biological processes under their control, remain insufficiently characterized despite their crucial function. Here, we identify a TFRN involved in the response to the key plant hormone auxin and define its impact on auxin-driven biological processes. To reconstruct the TFRN, we developed a three-step procedure, which is based on the integrated analysis of differentially expressed gene lists and a representative collection of transcription factor binding profiles. Its implementation is available as a part of the CisCross web server. With the new method, we distinguished two transcription factor subnetworks. The first operates before auxin treatment and is switched off upon hormone application, the second is switched on by the hormone. Moreover, we characterized the functioning of the auxin-regulated TFRN in control of chlorophyll and lignin biosynthesis, abscisic acid signaling, and ribosome biogenesis.

NetProphet 3: a machine learning framework for transcription factor network mapping and multi-omics integration

MOTIVATION

Many methods have been proposed for mapping the targets of transcription factors (TFs) from gene expression data. It is known that combining outputs from multiple methods can improve performance. To date, outputs have been combined by using either simplistic formulae, such as geometric mean, or carefully hand-tuned formulae that may not generalize well to new inputs. Finally, the evaluation of accuracy has been challenging due to the lack of genome-scale, ground-truth networks.

RESULTS

We developed NetProphet3, which combines scores from multiple analyses automatically, using a tree boosting algorithm trained on TF binding location data. We also developed three independent, genome-scale evaluation metrics. By these metrics, NetProphet3 is more accurate than other commonly used packages, including NetProphet 2.0, when gene expression data from direct TF perturbations are available. Furthermore, its integration mode can forge a consensus network from gene expression data and TF binding location data.

AVAILABILITY AND IMPLEMENTATION

All data and code are available at https://zenodo.org/record/7504131#.Y7Wu3i-B2x8.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The MITF regulatory network in melanoma

Bidirectional interactions between plastic tumor cells and the microenvironment critically impact tumor evolution and metastatic dissemination by enabling cancer cells to adapt to microenvironmental stresses by switching phenotype. In melanoma, a key determinant of phenotypic identity is the microphthalmia‐associated transcription factor MITF that promotes proliferation, suppresses senescence, and anticorrelates with immune infiltration and therapy resistance. What determines whether MITF can activate or repress genes associated with specific phenotypes, or how signaling regulating MITF might impact immune infiltration is poorly understood. Here, we find that MITF binding to genes associated with high MITF is via classical E/M‐box motifs, but genes downregulated when MITF is high contain FOS/JUN/AP1/ATF3 sites. Significantly, the repertoire of MITF‐interacting factors identified here includes JUN and ATF3 as well as many previously unidentified interactors. As high AP1 activity is a hallmark of MITF^Low, invasive, slow‐cycling, therapy resistant cells, the ability of MITF to repress AP1‐regulated genes provides an insight into how MITF establishes and maintains a pro‐proliferative phenotype. Moreover, although β‐catenin has been linked to immune exclusion, many Hallmark β‐catenin signaling genes are associated with immune infiltration. Instead, low MITF together with Notch signaling is linked to immune infiltration in both mouse and human melanoma tumors.

A Hierarchical Transcriptional Regulatory Network Required for Long-Term Thermal Stress Tolerance in an Industrial Saccharomyces cerevisiae Strain

Yeast cells suffer from continuous and long-term thermal stress during high-temperature ethanol fermentation. Understanding the mechanism of yeast thermotolerance is important not only for studying microbial stress biology in basic research but also for developing thermotolerant strains for industrial application. Here, we compared the effects of 23 transcription factor (TF) deletions on high-temperature ethanol fermentation and cell survival after heat shock treatment and identified three core TFs, Sin3p, Srb2p and Mig1p, that are involved in regulating the response to long-term thermotolerance. Further analyses of comparative transcriptome profiling of the core TF deletions and transcription regulatory associations revealed a hierarchical transcriptional regulatory network centered on these three TFs. This global transcriptional regulatory network provided a better understanding of the regulatory mechanism behind long-term thermal stress tolerance as well as potential targets for transcriptome engineering to improve the performance of high-temperature ethanol fermentation by an industrial Saccharomyces cerevisiae strain.

TIMEOR: a web-based tool to uncover temporal regulatory mechanisms from multi-omics data

Uncovering how transcription factors regulate their targets at DNA, RNA and protein levels over time is critical to define gene regulatory networks (GRNs) and assign mechanisms in normal and diseased states. RNA-seq is a standard method measuring gene regulation using an established set of analysis stages. However, none of the currently available pipeline methods for interpreting ordered genomic data (in time or space) use time-series models to assign cause and effect relationships within GRNs, are adaptive to diverse experimental designs, or enable user interpretation through a web-based platform. Furthermore, methods integrating ordered RNA-seq data with protein–DNA binding data to distinguish direct from indirect interactions are urgently needed. We present TIMEOR (Trajectory Inference and Mechanism Exploration with Omics data in R), the first web-based and adaptive time-series multi-omics pipeline method which infers the relationship between gene regulatory events across time. TIMEOR addresses the critical need for methods to determine causal regulatory mechanism networks by leveraging time-series RNA-seq, motif analysis, protein–DNA binding data, and protein-protein interaction networks. TIMEOR’s user-catered approach helps non-coders generate new hypotheses and validate known mechanisms. We used TIMEOR to identify a novel link between insulin stimulation and the circadian rhythm cycle. TIMEOR is available at https://github.com/ashleymaeconard/TIMEOR.git and http://timeor.brown.edu.

Inferring TF activities and activity regulators from gene expression data with constraints from TF perturbation data

Motivation

The activity of a transcription factor (TF) in a sample of cells is the extent to which it is exerting its regulatory potential. Many methods of inferring TF activity from gene expression data have been described, but due to the lack of appropriate large-scale datasets, systematic and objective validation has not been possible until now.

Results

We systematically evaluate and optimize the approach to TF activity inference in which a gene expression matrix is factored into a condition-independent matrix of control strengths and a condition-dependent matrix of TF activity levels. We find that expression data in which the activities of individual TFs have been perturbed are both necessary and sufficient for obtaining good performance. To a considerable extent, control strengths inferred using expression data from one growth condition carry over to other conditions, so the control strength matrices derived here can be used by others. Finally, we apply these methods to gain insight into the upstream factors that regulate the activities of yeast TFs Gcr2, Gln3, Gcn4 and Msn2.

Availability and implementation

Evaluation code and data are available at https://doi.org/10.5281/zenodo.4050573.

Supplementary information

Supplementary data are available at Bioinformatics online.

Network motif-based analysis of regulatory patterns in paralogous gene pairs

Current high-throughput experimental techniques make it feasible to infer gene regulatory interactions at the whole-genome level with reasonably good accuracy. Such experimentally inferred regulatory networks have become available for a number of simpler model organisms such as S. cerevisiae, and others. The availability of such networks provides an opportunity to compare gene regulatory processes at the whole genome level, and in particular, to assess similarity of regulatory interactions for homologous gene pairs either from the same or from different species. We present here a new technique for analyzing the regulatory interaction neighborhoods of paralogous gene pairs. Our central focus is the analysis of S. cerevisiae gene interaction graphs, which are of particular interest due to the ancestral whole-genome duplication (WGD) that allows to distinguish between paralogous transcription factors that are traceable to this duplication event and other paralogues. Similar analysis is also applied to E. coli and C. elegans networks. We compare paralogous gene pairs according to the presence and size of bi-fan arrays, classically associated in the literature with gene duplication, within other network motifs. We further extend this framework beyond transcription factor comparison to obtain topology-based similarity metrics based on the overlap of interaction neighborhoods applicable to most genes in a given organism. We observe that our network divergence metrics show considerably larger similarity between paralogues, especially those traceable to WGD. This is the case for both yeast and C. elegans, but not for E. coli regulatory network. While there is no obvious cross-species link between metrics, different classes of paralogues show notable differences in interaction overlap, with traceable duplications tending toward higher overlap compared to genes with shared protein families. Our findings indicate that divergence in paralogous interaction networks reflects a shared genetic origin, and that our approach may be useful for investigating structural similarity in the interaction networks of paralogous genes.

Runx1 and Runx3 drive progenitor to T-lineage transcriptome conversion in mouse T cell commitment via dynamic genomic site switching

Runt domain-related (Runx) transcription factors are essential for early T cell development in mice from uncommitted to committed stages. Single and double Runx knockouts via Cas9 show that target genes responding to Runx activity are not solely controlled by the dominant factor, Runx1. Instead, Runx1 and Runx3 are coexpressed in single cells; bind to highly overlapping genomic sites; and have redundant, collaborative functions regulating genes pivotal for T cell development. Despite stable combined expression levels across pro-T cell development, Runx1 and Runx3 preferentially activate and repress genes that change expression dynamically during lineage commitment, mostly activating T-lineage genes and repressing multipotent progenitor genes. Furthermore, most Runx target genes are sensitive to Runx perturbation only at one stage and often respond to Runx more for expression transitions than for maintenance. Contributing to this highly stage-dependent gene regulation function, Runx1 and Runx3 extensively shift their binding sites during commitment. Functionally distinct Runx occupancy sites associated with stage-specific activation or repression are also distinguished by different patterns of partner factor cobinding. Finally, Runx occupancies change coordinately at numerous clustered sites around positively or negatively regulated targets during commitment. This multisite binding behavior may contribute to a developmental "ratchet" mechanism making commitment irreversible.

Dynamic proteomic profiling of human periodontal ligament stem cells during osteogenic differentiation

Background

Human periodontal ligament stem cells (hPDLSCs) are ideal seed cells for periodontal regeneration. A greater understanding of the dynamic protein profiles during osteogenic differentiation contributed to the improvement of periodontal regeneration tissue engineering.

Methods

Tandem Mass Tag quantitative proteomics was utilized to reveal the temporal protein expression pattern during osteogenic differentiation of hPDLSCs on days 0, 3, 7 and 14. Differentially expressed proteins (DEPs) were clustered and functional annotated by Gene Ontology (GO) terms. Pathway enrichment analysis was performed based on the Kyoto Encyclopedia of Genes and Genomes database, followed by the predicted activation using Ingenuity Pathway Analysis software. Interaction networks of redox-sensitive signalling pathways and oxidative phosphorylation (OXPHOS) were conducted and the hub protein SOD2 was validated with western blotting.

Results

A total of 1024 DEPs were identified and clustered in 5 distinctive clusters representing dynamic tendencies. The GO enrichment results indicated that proteins with different tendencies show different functions. Pathway enrichment analysis found that OXPHOS was significantly involved, which further predicted continuous activation. Redox-sensitive signalling pathways with dynamic activation status showed associations with OXPHOS to various degrees, especially the sirtuin signalling pathway. SOD2, an important component of the sirtuin pathway, displays a persistent increase during osteogenesis. Data are available via ProteomeXchange with identifier PXD020908.

Conclusion

This is the first in-depth dynamic proteomic analysis of osteogenic differentiation of hPDLSCs. It demonstrated a dynamic regulatory mechanism of hPDLSC osteogenesis and might provide a new perspective for research on periodontal regeneration.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-020-02123-6.

An Early Stage Researcher's Primer on Systems Medicine Terminology

Background: Systems Medicine is a novel approach to medicine, that is, an interdisciplinary field that considers the human body as a system, composed of multiple parts and of complex relationships at multiple levels, and further integrated into an environment. Exploring Systems Medicine implies understanding and combining concepts coming from diametral different fields, including medicine, biology, statistics, modeling and simulation, and data science. Such heterogeneity leads to semantic issues, which may slow down implementation and fruitful interaction between these highly diverse fields. Methods: In this review, we collect and explain more than100 terms related to Systems Medicine. These include both modeling and data science terms and basic systems medicine terms, along with some synthetic definitions, examples of applications, and lists of relevant references. Results: This glossary aims at being a first aid kit for the Systems Medicine researcher facing an unfamiliar term, where he/she can get a first understanding of them, and, more importantly, examples and references for digging into the topic.

A time-resolved proteomic analysis of transcription factors regulating adipogenesis of human adipose derived stem cells

Adipogenesis is one of the key processes during obesity development. Better understanding of this process could advance our knowledge on obesity and its treatment. Transcription factors (TFs) are master regulators during adipogenesis, however, a system-wide analysis of TFs dynamic proteome during adipogenesis is lacking. Here, we profiled 472 TFs and systematically elucidated their roles during the first 7 days of adipogenesis of human adipose-derived stem cells (hADSCs) on proteome scale. We identified two main and four sub-phases during adipogenesis. The commitment phase (0 h-8 h) mainly mediated stem cell proliferation, differentiation and chromatin remodeling, in which proteins of SWI/SNF family are the key centroid nodes. The determination phase (1D-7D) predominately regulated fat cell differentiation and response to lipid and oxygen, which could be associated with terminal differentiation of adipocyte and responsible for maturation. PPARγ, CREB1 and MYC are the centroid nodes of this phase. Remarkably, we identified and verified three TFs (BATF3, MAFF and MXD4) as novel regulators of adipogenesis, whose over-expression could inhibit adipogenesis of hADSCs in vitro. Overall, our study provided a valuable TFs resource to understand the complex process of adipogenesis.

The Genome-Scale Integrated Networks in Microorganisms

The genome-scale cellular network has become a necessary tool in the systematic analysis of microbes. In a cell, there are several layers (i.e., types) of the molecular networks, for example, genome-scale metabolic network (GMN), transcriptional regulatory network (TRN), and signal transduction network (STN). It has been realized that the limitation and inaccuracy of the prediction exist just using only a single-layer network. Therefore, the integrated network constructed based on the networks of the three types attracts more interests. The function of a biological process in living cells is usually performed by the interaction of biological components. Therefore, it is necessary to integrate and analyze all the related components at the systems level for the comprehensively and correctly realizing the physiological function in living organisms. In this review, we discussed three representative genome-scale cellular networks: GMN, TRN, and STN, representing different levels (i.e., metabolism, gene regulation, and cellular signaling) of a cell’s activities. Furthermore, we discussed the integration of the networks of the three types. With more understanding on the complexity of microbial cells, the development of integrated network has become an inevitable trend in analyzing genome-scale cellular networks of microorganisms.

Genomic and molecular control of cell type and cell type conversions

Organisms are made of a limited number of cell types that combine to form higher order tissues and organs. Cell types have traditionally been defined by their morphologies or biological activity, yet the underlying molecular controls of cell type remain unclear. The onset of single cell technologies, and more recently genomics (particularly single cell genomics), has substantially increased the understanding of the concept of cell type, but has also increased the complexity of this understanding. These new technologies have added a new genome wide molecular dimension to the description of cell type, with genome-wide expression and epigenetic data acting as a cell type ‘fingerprint’ to describe the cell state. Using these genomic fingerprints cell types are being increasingly defined based on specific genomic and molecular criteria, without necessarily a distinct biological function. In this review, we will discuss the molecular definitions of cell types and cell type control, and particularly how endogenous and exogenous transcription factors can control cell types and cell type conversions.

Mammalian Transcription Factor Networks: Recent Advances in Interrogating Biological Complexity

Transcription factor (TF) networks are a key determinant of cell fate decisions in mammalian development and adult tissue homeostasis and are frequently corrupted in disease. However, our inability to experimentally resolve and interrogate the complexity of mammalian TF networks has hampered the progress in this field. Recent technological advances, in particular large-scale genome-wide approaches, single-cell methodologies, live-cell imaging, and genome editing, are emerging as important technologies in TF network biology. Several recent studies even suggest a need to re-evaluate established models of mammalian TF networks. Here, we provide a brief overview of current and emerging methods to define mammalian TF networks. We also discuss how these emerging technologies facilitate new ways to interrogate complex TF networks, consider the current open questions in the field, and comment on potential future directions and biomedical applications.

NetProphet 2.0: mapping transcription factor networks by exploiting scalable data resources

Motivation

Cells process information, in part, through transcription factor (TF) networks, which control the rates at which individual genes produce their products. A TF network map is a graph that indicates which TFs bind and directly regulate each gene. Previous work has described network mapping algorithms that rely exclusively on gene expression data and ‘integrative’ algorithms that exploit a wide range of data sources including chromatin immunoprecipitation sequencing (ChIP-seq) of many TFs, genome-wide chromatin marks, and binding specificities for many TFs determined in vitro. However, such resources are available only for a few major model systems and cannot be easily replicated for new organisms or cell types.

Results

We present NetProphet 2.0, a ‘data light’ algorithm for TF network mapping, and show that it is more accurate at identifying direct targets of TFs than other, similarly data light algorithms. In particular, it improves on the accuracy of NetProphet 1.0, which used only gene expression data, by exploiting three principles. First, combining multiple approaches to network mapping from expression data can improve accuracy relative to the constituent approaches. Second, TFs with similar DNA binding domains bind similar sets of target genes. Third, even a noisy, preliminary network map can be used to infer DNA binding specificities from promoter sequences and these inferred specificities can be used to further improve the accuracy of the network map.

Availability and implementation

Source code and comprehensive documentation are freely available at https://github.com/yiming-kang/NetProphet_2.0.

Supplementary information

Supplementary data are available at Bioinformatics online.

Lineage specification of early embryos and embryonic stem cells at the dawn of enabling technologies

Establishment of progenitor cell populations and lineage diversity during embryogenesis and the differentiation of pluripotent stem cells is a fascinating and intricate biological process. Conceptually, an understanding of this developmental process provides a framework to integrate stem-cell pluripotency, cell competence and differentiating potential with the activity of extrinsic and intrinsic molecular determinants. The recent advent of enabling technologies of high-resolution transcriptome analysis at the cellular, population and spatial levels proffers the capability of gaining deeper insights into the attributes of the gene regulatory network and molecular signaling in lineage specification and differentiation. In this review, we provide a snapshot of the emerging enabling genomic technologies that contribute to the study of development and stem-cell biology.

Transcriptional events defining plant immune responses

Rapid and massive transcriptional reprogramming upon pathogen recognition is the decisive step in plant-phytopathogen interactions. Plant transcription factors (TFs) are key players in this process but they require a suite of other context-specific co-regulators to establish sensory transcription regulatory networks to bring about host immunity. Molecular, genetic and biochemical studies, particularly in the model plants Arabidopsis and rice, are continuously uncovering new components of the transcriptional machinery that can selectively impact host resistance toward a diverse range of pathogens. Moreover, detailed studies on key immune regulators, such as WRKY TFs and NPR1, are beginning to reveal the underlying mechanisms by which defense hormones influence the function of these factors. Here we provide a short update on such recent developments.

A cross-species approach to identify transcriptional regulators exemplified for Dnajc22 and Hnf4a

There is an enormous need to make better use of the ever increasing wealth of publicly available genomic information and to utilize the tremendous progress in computational approaches in the life sciences. Transcriptional regulation of protein-coding genes is a major mechanism of controlling cellular functions. However, the myriad of transcription factors potentially controlling transcription of any given gene makes it often difficult to quickly identify the biological relevant transcription factors. Here, we report on the identification of Hnf4a as a major transcription factor of the so far unstudied DnaJ heat shock protein family (Hsp40) member C22 (Dnajc22). We propose an approach utilizing recent advances in computational biology and the wealth of publicly available genomic information guiding the identification of potential transcription factor candidates together with wet-lab experiments validating computational models. More specifically, the combined use of co-expression analyses based on self-organizing maps with sequence-based transcription factor binding prediction led to the identification of Hnf4a as the potential transcriptional regulator for Dnajc22 which was further corroborated using publicly available datasets on Hnf4a. Following this procedure, we determined its functional binding site in the murine Dnajc22 locus using ChIP-qPCR and luciferase assays and verified this regulatory loop in fruitfly, zebrafish, and humans.

Developmental Control of NRAMP1 (SLC11A1) Expression in Professional Phagocytes

NRAMP1 (SLC11A1) is a professional phagocyte membrane importer of divalent metals that contributes to iron recycling at homeostasis and to nutritional immunity against infection. Analyses of data generated by several consortia and additional studies were integrated to hypothesize mechanisms restricting NRAMP1 expression to mature phagocytes. Results from various epigenetic and transcriptomic approaches were collected for mesodermal and hematopoietic cell types and compiled for combined analysis with results of genetic studies associating single nucleotide polymorphisms (SNPs) with variations in NRAMP1 expression (eQTLs). Analyses establish that NRAMP1 is part of an autonomous topologically associated domain delimited by ubiquitous CCCTC-binding factor (CTCF) sites. NRAMP1 locus contains five regulatory regions: a predicted super-enhancer (S-E) key to phagocyte-specific expression; the proximal promoter; two intronic areas, including 3' inhibitory elements that restrict expression during development; and a block of upstream sites possibly extending the S-E domain. Also the downstream region adjacent to the 3' CTCF locus boundary may regulate expression during hematopoiesis. Mobilization of the locus 14 predicted transcriptional regulatory elements occurs in three steps, beginning with hematopoiesis; at the onset of myelopoiesis and through myelo-monocytic differentiation. Basal expression level in mature phagocytes is further influenced by genetic variation, tissue environment, and in response to infections that induce various epigenetic memories depending on microorganism nature. Constitutively associated transcription factors (TFs) include CCAAT enhancer binding protein beta (C/EBPb), purine rich DNA binding protein (PU.1), early growth response 2 (EGR2) and signal transducer and activator of transcription 1 (STAT1) while hypoxia-inducible factors (HIFs) and interferon regulatory factor 1 (IRF1) may stimulate iron acquisition in pro-inflammatory conditions. Mouse orthologous locus is generally conserved; chromatin patterns typify a de novo myelo-monocytic gene whose expression is tightly controlled by TFs Pu.1, C/ebps and Irf8; Irf3 and nuclear factor NF-kappa-B p 65 subunit (RelA) regulate expression in inflammatory conditions. Functional differences in the determinants identified at these orthologous loci imply that species-specific mechanisms control gene expression.

Model-based transcriptome engineering promotes a fermentative transcriptional state in yeast

The ability to rationally manipulate the transcriptional states of cells would be of great use in medicine and bioengineering. We have developed an algorithm, NetSurgeon, which uses genome-wide gene-regulatory networks to identify interventions that force a cell toward a desired expression state. We first validated NetSurgeon extensively on existing datasets. Next, we used NetSurgeon to select transcription factor deletions aimed at improving ethanol production in Saccharomyces cerevisiae cultures that are catabolizing xylose. We reasoned that interventions that move the transcriptional state of cells using xylose toward that of cells producing large amounts of ethanol from glucose might improve xylose fermentation. Some of the interventions selected by NetSurgeon successfully promoted a fermentative transcriptional state in the absence of glucose, resulting in strains with a 2.7-fold increase in xylose import rates, a 4-fold improvement in xylose integration into central carbon metabolism, or a 1.3-fold increase in ethanol production rate. We conclude by presenting an integrated model of transcriptional regulation and metabolic flux that will enable future efforts aimed at improving xylose fermentation to prioritize functional regulators of central carbon metabolism.