Estimating the activity of transcription factors by the effect on their target genes

Genomic analysis of progenitors in viral infection implicates glucocorticoids as suppressors of plasmacytoid dendritic cell generation

Plasmacytoid Dendritic cells (pDCs) are the most potent producers of interferons, which are critical antiviral cytokines. pDC development is, however, compromised following a viral infection, and this phenomenon, as well as its relationship to conventional (c)DC development is still incompletely understood. By using lymphocytic choriomeningitis virus (LCMV) infection in mice as a model system, we observed that DC progenitors skewed away from pDC and toward cDC development during in vivo viral infection. Subsequent characterization of the transcriptional and epigenetic landscape of fms-like tyrosine kinase 3+ (Flt3+) DC progenitors and follow-up studies revealed increased apoptosis and reduced proliferation in different individual DC-progenitors as well as a profound type I interferon (IFN-I)-dependent ablation of pre-pDCs, but not pre-DC precursors, after both acute and chronic LCMV infections. In addition, integrated genomic analysis identified altered activity of 34 transcription factors in Flt3+ DC progenitors from infected mice, including two regulators of Glucocorticoid (GC) responses. Subsequent studies demonstrated that addition of GCs to DC progenitors led to downregulated pDC-primed-genes while upregulating cDC-primed-genes, and that endogenous GCs selectively decreased pDC, but not cDC, numbers upon in vivo LCMV infection. These findings demonstrate a significant ablation of pre-pDCs in infected mice and identify GCs as suppressors of pDC generation from early progenitors. This provides a potential explanation for the impaired pDC development following viral infection and links pDC numbers to the hypothalamic-pituitary-adrenal axis.

Multi-Omics Atlas-Assisted Discovery of Transcription Factors for Selective T Cell State Programming

Transcription factors (TFs) regulate the differentiation of T cells into diverse states with distinct functionalities. To precisely program desired T cell states in viral infections and cancers, we generated a comprehensive transcriptional and epigenetic atlas of nine CD8 + T cell differentiation states for TF activity prediction. Our analysis catalogued TF activity fingerprints of each state, uncovering new regulatory mechanisms that govern selective cell state differentiation. Leveraging this platform, we focused on two critical T cell states in tumor and virus control: terminally exhausted T cells (TEX term ), which are dysfunctional, and tissue-resident memory T cells (T RM ), which are protective. Despite their functional differences, these states share significant transcriptional and anatomical similarities, making it both challenging and essential to engineer T cells that avoid TEX term differentiation while preserving beneficial T RM characteristics. Through in vivo CRISPR screening combined with single-cell RNA sequencing (Perturb-seq), we validated the specific TFs driving the TEX term state and confirmed the accuracy of TF specificity predictions. Importantly, we discovered novel TEX term -specific TFs such as ZSCAN20, JDP2, and ZFP324. The deletion of these TEX term -specific TFs in T cells enhanced tumor control and synergized with immune checkpoint blockade. Additionally, this study identified multi-state TFs like HIC1 and GFI1, which are vital for both TEX term and T RM states. Furthermore, our global TF community analysis and Perturb-seq experiments revealed how TFs differentially regulate key processes in T RM and TEX term cells, uncovering new biological pathways like protein catabolism that are specifically linked to TEX term differentiation. In summary, our platform systematically identifies TF programs across diverse T cell states, facilitating the engineering of specific T cell states to improve tumor control and providing insights into the cellular mechanisms underlying their functional disparities.

Genomic Analysis of Progenitors in Viral Infection Implicates Glucocorticoids as Suppressors of Plasmacytoid Dendritic Cell Generation

Plasmacytoid Dendritic cells (pDCs) are the most potent producers of interferons, which are critical antiviral cytokines. pDC development is, however, compromised following a viral infection, and this phenomenon, as well as its relationship to conventional (c)DC development is still incompletely understood. By using lymphocytic choriomeningitis virus (LCMV) infection in mice as a model system, we observed that DC progenitors skewed away from pDC and towards cDC development during in vivo viral infection. Subsequent characterization of the transcriptional and epigenetic landscape of fms-like tyrosine kinase 3 + (Flt3 + ) DC progenitors and follow-up studies revealed increased apoptosis and reduced proliferation in different individual DC-progenitors as well as a profound IFN-I-dependent ablation of pre-pDCs, but not pre-DC precursor, after both acute and chronic LCMV infections. In addition, integrated genomic analysis identified altered activity of 34 transcription factors in Flt3 + DC progenitors from infected mice, including two regulators of Glucocorticoid (GC) responses. Subsequent studies demonstrated that addition of GCs to DC progenitors led to downregulated pDC-primed-genes while upregulating cDC-primed-genes, and that endogenous GCs selectively decreased pDC, but not cDC, numbers upon in-vivo LCMV infection. These findings demonstrate a significant ablation of pre-pDCs in infected mice and identify GCs as suppressors of pDC generation from early progenitors. This provides an explanation for the impaired pDC development following viral infection and links pDC generation to the hypothalamic-pituitary-adrenal axis.

Significance Statement

Plasmacytoid dendritic cells (pDCs) play critical roles in antiviral responses. However, adaptations of DC progenitors lead to compromised pDC generation after viral infection. Here, we characterized the transcriptional and epigenetic landscapes of DC progenitors after infection. We observed widespread changes in gene expression and chromatin accessibility, reflecting shifts in proliferation, apoptosis, and differentiation potential into various DC subsets. Notably, we identified alterations in the predicted activity of 34 transcription factors, including two regulators of glucocorticoid responses. Our data demonstrate that glucocorticoids inhibit pDC generation by reprogramming DC progenitors. These findings establish a molecular framework for understanding how DC progenitors adapt to infection and highlight the role of glucocorticoid signaling in this process.

Gene regulatory networks in disease and ageing

The precise control of gene expression is required for the maintenance of cellular homeostasis and proper cellular function, and the declining control of gene expression with age is considered a major contributor to age-associated changes in cellular physiology and disease. The coordination of gene expression can be represented through models of the molecular interactions that govern gene expression levels, so-called gene regulatory networks. Gene regulatory networks can represent interactions that occur through signal transduction, those that involve regulatory transcription factors, or statistical models of gene-gene relationships based on the premise that certain sets of genes tend to be coexpressed across a range of conditions and cell types. Advances in experimental and computational technologies have enabled the inference of these networks on an unprecedented scale and at unprecedented precision. Here, we delineate different types of gene regulatory networks and their cell-biological interpretation. We describe methods for inferring such networks from large-scale, multi-omics datasets and present applications that have aided our understanding of cellular ageing and disease mechanisms.

The Blueprint of Logical Decisions in a NF-κB Signaling System

Nearly identical cells can exhibit substantially different responses to the same stimulus that causes phenotype diversity. Such interplay between phenotype diversity and the architecture of regulatory circuits is crucial since it determines the state of a biological cell. Here, we theoretically analyze how the circuit blueprints of NF-κB in cellular environments are formed and their role in determining the cells' metabolic state. The NF-κB is a collective name for a developmental conserved family of five different transcription factors that can form homodimers or heterodimers and often promote DNA looping to reprogram the inflammatory gene response. The NF-κB controls many biological functions, including cellular differentiation, proliferation, migration, and survival. Our model shows that nuclear localization of NF-κB differentially promotes logic operations such as AND, NAND, NOR, and OR in its regulatory network. Through the quantitative thermodynamic model of transcriptional regulation and systematic variation of promoter-enhancer interaction modes, we can account for the origin of various logic gates as formed in the NF-κB system. We further show that the interconversion or switching of logic gates yielded under systematic variations of the stimuli activity and DNA looping parameters. Such computation occurs in regulatory and signaling pathways in individual cells at a molecular scale, which one can exploit to design a biomolecular computer.

Network-based drug repurposing for schizophrenia

Despite recent progress, the challenges in drug discovery for schizophrenia persist. However, computational drug repurposing has gained popularity as it leverages the wealth of expanding biomedical databases. Network analyses provide a comprehensive understanding of transcription factor (TF) regulatory effects through gene regulatory networks, which capture the interactions between TFs and target genes by integrating various lines of evidence. Using the PANDA algorithm, we examined the topological variances in TF-gene regulatory networks between individuals with schizophrenia and healthy controls. This algorithm incorporates binding motifs, protein interactions, and gene co-expression data. To identify these differences, we subtracted the edge weights of the healthy control network from those of the schizophrenia network. The resulting differential network was then analysed using the CLUEreg tool in the GRAND database. This tool employs differential network signatures to identify drugs that potentially target the gene signature associated with the disease. Our analysis utilised a large RNA-seq dataset comprising 532 post-mortem brain samples from the CommonMind project. We constructed co-expression gene regulatory networks for both schizophrenia cases and healthy control subjects, incorporating 15,831 genes and 413 overlapping TFs. Through drug repurposing, we identified 18 promising candidates for repurposing as potential treatments for schizophrenia. The analysis of TF-gene regulatory networks revealed that the TFs in schizophrenia predominantly regulate pathways associated with energy metabolism, immune response, cell adhesion, and thyroid hormone signalling. These pathways represent significant targets for therapeutic intervention. The identified drug repurposing candidates likely act through TF-targeted pathways. These promising candidates, particularly those with preclinical evidence such as rimonabant and kaempferol, warrant further investigation into their potential mechanisms of action and efficacy in alleviating the symptoms of schizophrenia.

HIF-1α Mediates Immunosuppression and Chemoresistance in Colorectal Cancer by Inhibiting CXCL9, -10 and -11

Uncertainty exists regarding the mechanisms by which hypoxia-inducible factors (HIFs) control CD8+T-cell migration into tumor microenvironments. Here, we found that HIF-1α knockdown or overexpression resulted in increased or decreased CXCL9, -10, and -11 expression in vitro, respectively. Gene Set Variation Analysis revealed that elevated HIF-1α levels correlated with a poor prognosis, severe pathological stage, and an absence of CD8+ T cells in the tumor microenvironment in colorectal cancer (CRC) patients. HIF-1α was inversely associated with pathways beneficial to anti-tumor immunotherapy and cytokine/chemokine function. In vivo, inhibiting HIF-1α or its upstream regulator BIRC2 significantly suppressed tumor growth and promoted CD8+ T-cell infiltration. CXCR3 neutralizing antibodies reversed these effects, implicating the involvement of CXCL9, -10, and -11/CXCR3 axis. The presence of HIF-1α weakened the upregulation of CXCL9, -10, and -11 by bleomycin and doxorubicin. Combining HIF-1α inhibition with bleomycin promoted CD8+ T-cell infiltration and tumor suppression in vivo. Moreover, doxorubicin could upregulate CXCL9, -10 and -11 by suppressing HIF-1α. Our findings highlight the potential of HIF-1α inhibition to improve CRC microenvironments and increase chemotherapy sensitivity.

Structure-primed embedding on the transcription factor manifold enables transparent model architectures for gene regulatory network and latent activity inference

BACKGROUND

Modeling of gene regulatory networks (GRNs) is limited due to a lack of direct measurements of genome-wide transcription factor activity (TFA) making it difficult to separate covariance and regulatory interactions. Inference of regulatory interactions and TFA requires aggregation of complementary evidence. Estimating TFA explicitly is problematic as it disconnects GRN inference and TFA estimation and is unable to account for, for example, contextual transcription factor-transcription factor interactions, and other higher order features. Deep-learning offers a potential solution, as it can model complex interactions and higher-order latent features, although does not provide interpretable models and latent features.

RESULTS

We propose a novel autoencoder-based framework, StrUcture Primed Inference of Regulation using latent Factor ACTivity (SupirFactor) for modeling, and a metric, explained relative variance (ERV), for interpretation of GRNs. We evaluate SupirFactor with ERV in a wide set of contexts. Compared to current state-of-the-art GRN inference methods, SupirFactor performs favorably. We evaluate latent feature activity as an estimate of TFA and biological function in S. cerevisiae as well as in peripheral blood mononuclear cells (PBMC).

CONCLUSION

Here we present a framework for structure-primed inference and interpretation of GRNs, SupirFactor, demonstrating interpretability using ERV in multiple biological and experimental settings. SupirFactor enables TFA estimation and pathway analysis using latent factor activity, demonstrated here on two large-scale single-cell datasets, modeling S. cerevisiae and PBMC. We find that the SupirFactor model facilitates biological analysis acquiring novel functional and regulatory insight.

A Bayesian noisy logic model for inference of transcription factor activity from single cell and bulk transcriptomic data

The advent of high-throughput sequencing has made it possible to measure the expression of genes at relatively low cost. However, direct measurement of regulatory mechanisms, such as transcription factor (TF) activity is still not readily feasible in a high-throughput manner. Consequently, there is a need for computational approaches that can reliably estimate regulator activity from observable gene expression data. In this work, we present a noisy Boolean logic Bayesian model for TF activity inference from differential gene expression data and causal graphs. Our approach provides a flexible framework to incorporate biologically motivated TF-gene regulation logic models. Using simulations and controlled over-expression experiments in cell cultures, we demonstrate that our method can accurately identify TF activity. Moreover, we apply our method to bulk and single cell transcriptomics measurements to investigate transcriptional regulation of fibroblast phenotypic plasticity. Finally, to facilitate usage, we provide user-friendly software packages and a web-interface to query TF activity from user input differential gene expression data: https://umbibio.math.umb.edu/nlbayes/.

Identifying genetic regulatory variants that affect transcription factor activity

Genetic variants affecting gene expression levels in humans have been mapped in the Genotype-Tissue Expression (GTEx) project. Trans-acting variants impacting many genes simultaneously through a shared transcription factor (TF) are of particular interest. Here, we developed a generalized linear model (GLM) to estimate protein-level TF activity levels in an individual-specific manner from GTEx RNA sequencing (RNA-seq) profiles. It uses observed differential gene expression after TF perturbation as a predictor and, by analyzing differential expression within pairs of neighboring genes, controls for the confounding effect of variation in chromatin state along the genome. We inferred genotype-specific activities for 55 TFs across 49 tissues. Subsequently performing genome-wide association analysis on this virtual trait revealed TF activity quantitative trait loci (aQTLs) that, as a set, are enriched for functional features. Altogether, the set of tools we introduce here highlights the potential of genetic association studies for cellular endophenotypes based on a network-based multi-omics approach. The transparent peer review record is available.

Characterization of cell-fate decision landscapes by estimating transcription factor dynamics

Time-specific modulation of gene expression during differentiation by transcription factors promotes cell diversity. However, estimating their dynamic regulatory activity at the single-cell level and in a high-throughput manner remains challenging. We present FateCompass, an integrative approach that utilizes single-cell transcriptomics data to identify lineage-specific transcription factors throughout differentiation. By combining a probabilistic framework with RNA velocities or differentiation potential, we estimate transition probabilities, while a linear model of gene regulation is employed to compute transcription factor activities. Considering dynamic changes and correlations of expression and activities, FateCompass identifies lineage-specific regulators. Our validation using in silico data and application to pancreatic endocrine cell differentiation datasets highlight both known and potentially novel lineage-specific regulators. Notably, we uncovered undescribed transcription factors of an enterochromaffin-like population during in vitro differentiation toward ß-like cells. FateCompass provides a valuable framework for hypothesis generation, advancing our understanding of the gene regulatory networks driving cell-fate decisions.

A Bayesian Noisy Logic Model for Inference of Transcription Factor Activity from Single Cell and Bulk Transcriptomic Data

The advent of high-throughput sequencing has made it possible to measure the expression of genes at relatively low cost. However, direct measurement of regulatory mechanisms, such as Transcription Factor (TF) activity is still not readily feasible in a high-throughput manner. Consequently, there is a need for computational approaches that can reliably estimate regulator activity from observable gene expression data. In this work, we present a noisy Boolean logic Bayesian model for TF activity inference from differential gene expression data and causal graphs. Our approach provides a flexible framework to incorporate biologically motivated TF-gene regulation logic models. Using simulations and controlled over-expression experiments in cell cultures, we demonstrate that our method can accurately identify TF activity. Moreover, we apply our method to bulk and single cell transcriptomics measurements to investigate transcriptional regulation of fibroblast phenotypic plasticity. Finally, to facilitate usage, we provide user-friendly software packages and a web-interface to query TF activity from user input differential gene expression data: https://umbibio.math.umb.edu/nlbayes/.

Integrated analysis of ovarian cancer patients from prospective transcription factor activity reveals subtypes of prognostic significance

Transcription factors are protein molecules that act as regulators of gene expression. Aberrant protein activity of transcription factors can have a significant impact on tumor progression and metastasis in tumor patients. In this study, 868 immune-related transcription factors were identified from the transcription factor activity profile of 1823 ovarian cancer patients. The prognosis-related transcription factors were identified through univariate Cox analysis and random survival tree analysis, and two distinct clustering subtypes were subsequently derived based on these transcription factors. We assessed the clinical significance and genomics landscape of the two clustering subtypes and found statistically significant differences in prognosis, response to immunotherapy, and chemotherapy among ovarian cancer patients with different subtypes. Multi-scale Embedded Gene Co-expression Network Analysis was used to identify differential gene modules between the two clustering subtypes, which allowed us to conduct further analysis of biological pathways that exhibited significant differences between them. Finally, a ceRNA network was constructed to analyze lncRNA-miRNA-mRNA regulatory pairs with differential expression levels between two clustering subtypes. We expected that our study may provide some useful references for stratifying and treating patients with ovarian cancer.

Structure primed embedding on the transcription factor manifold enables transparent model architectures for gene regulatory network and latent activity inference

The modeling of gene regulatory networks (GRNs) is limited due to a lack of direct measurements of regulatory features in genome-wide screens. Most GRN inference methods are therefore forced to model relationships between regulatory genes and their targets with expression as a proxy for the upstream independent features, complicating validation and predictions produced by modeling frameworks. Separating covariance and regulatory influence requires aggregation of independent and complementary sets of evidence, such as transcription factor (TF) binding and target gene expression. However, the complete regulatory state of the system, e.g. TF activity (TFA) is unknown due to a lack of experimental feasibility, making regulatory relations difficult to infer. Some methods attempt to account for this by modeling TFA as a latent feature, but these models often use linear frameworks that are unable to account for non-linearities such as saturation, TF-TF interactions, and other higher order features. Deep learning frameworks may offer a solution, as they are capable of modeling complex interactions and capturing higher-order latent features. However, these methods often discard central concepts in biological systems modeling, such as sparsity and latent feature interpretability, in favor of increased model complexity. We propose a novel deep learning autoencoder-based framework, StrUcture Primed Inference of Regulation using latent Factor ACTivity (SupirFactor), that scales to single cell genomic data and maintains interpretability to perform GRN inference and estimate TFA as a latent feature. We demonstrate that SupirFactor outperforms current leading GRN inference methods, predicts biologically relevant TFA and elucidates functional regulatory pathways through aggregation of TFs.

ZmCCT10-relayed photoperiod sensitivity regulates natural variation in the arithmetical formation of male germinal cells in maize

Extensive mutational screening studies have documented genes regulating anther and pollen development. Knowledge concerning how formation of male germinal cell is arithmetically controlled in natural populations, under different environmental conditions, is lacking. We counted pollen number within a single anther and a maize-teosinte BC₂ S₃ recombinant inbred line population to identify ZmCCT10 as a major determinant of pollen number variation. ZmCCT10 was originally identified as a photoperiod-sensitive negative regulator of flowering. ZmCCT10 inactivation, after transposon insertion within its promoter, is proposed to have accelerated maize spread toward higher latitudes, thus allowing temperate maize to flower under long-day conditions. We showed that the active ZmCCT10 allele decreased pollen formation. As different active and inactive ZmCCT10 alleles have been found in natural maize populations, this represents the first report of a gene controlling pollen number in a crop natural population. These findings suggest that higher pollen number, which provides a competitive advantage in open-pollinated populations, may have been one of the major driving forces for the selection of an inactive ZmCCT10 allele during tropical maize domestication. We provide evidence that ZmCCT10 has opposite effects on cell proliferation of archesporial and tapetum cells and it modulates expression of key regulators during early anther development.

Alternative splicing: transcriptional regulatory network in agroforestry

Alternative splicing (AS) in plants plays a key role in regulating the expression of numerous transcripts from a single gene in a regulatory pathway. Variable concentrations of growth regulatory hormones and external stimuli trigger alternative splicing to switch among different growth stages and adapt to environmental stresses. In the AS phenomenon, a spliceosome causes differential transcriptional modifications in messenger RNA (mRNAs), resulting in partial or complete retention of one or more introns as compared to fully spliced mRNA. Differentially expressed proteins translated from intron-retaining messenger RNA (mRNA^ir) perform vital functions in the feedback mechanism. At the post-transcriptional level, AS causes the remodeling of transcription factors (TFs) by the addition or deletion of binding domains to activate and/or repress transcription. In this study, we have summarized the specific role of AS in the regulation of gene expression through repression and activation of the transcriptional regulatory network under external stimuli and switch among developmental stages.

Coordinated control of neuronal differentiation and wiring by sustained transcription factors

The large diversity of cell types in nervous systems presents a challenge in identifying the genetic mechanisms that encode it. Here, we report that nearly 200 distinct neurons in the Drosophila visual system can each be defined by unique combinations of on average 10 continuously expressed transcription factors. We show that targeted modifications of this terminal selector code induce predictable conversions of neuronal fates that appear morphologically and transcriptionally complete. Cis-regulatory analysis of open chromatin links one of these genes to an upstream patterning factor that specifies neuronal fates in stem cells. Experimentally validated network models describe the synergistic regulation of downstream effectors by terminal selectors and ecdysone signaling during brain wiring. Our results provide a generalizable framework of how specific fates are implemented in postmitotic neurons.

Systems-level identification of key transcription factors in immune cell specification

Transcription factors (TFs) are crucial for regulating cell differentiation during the development of the immune system. However, the key TFs for orchestrating the specification of distinct immune cells are not fully understood. Here, we integrated the transcriptomic and epigenomic measurements in 73 mouse and 61 human primary cell types, respectively, that span the immune cell differentiation pathways. We constructed the cell-type-specific transcriptional regulatory network and assessed the global importance of TFs based on the Taiji framework, which is a method we have previously developed that can infer the global impact of TFs using integrated transcriptomic and epigenetic data. Integrative analysis across cell types revealed putative driver TFs in cell lineage-specific differentiation in both mouse and human systems. We have also identified TF combinations that play important roles in specific developmental stages. Furthermore, we validated the functions of predicted novel TFs in murine CD8+ T cell differentiation and showed the importance of Elf1 and Prdm9 in the effector versus memory T cell fate specification and Kdm2b and Tet3 in promoting differentiation of CD8+ tissue resident memory (Trm) cells, validating the approach. Thus, we have developed a bioinformatic approach that provides a global picture of the regulatory mechanisms that govern cellular differentiation in the immune system and aids the discovery of novel mechanisms in cell fate decisions.

PITX1 Is a Regulator of TERT Expression in Prostate Cancer with Prognostic Power

Simple Summary

Most prostate cancer is of an indolent form and is curable. However, some prostate cancer belongs to rather aggressive subtypes leading to metastasis and death, and immediate therapy is mandatory. However, for these, the therapeutic options are highly invasive, such as radical prostatectomy, radiation or brachytherapy. Hence, a precise diagnosis of these tumor subtypes is needed, and the thus far applied diagnostic means are insufficient for this. Besides this, for their endless cell divisions, prostate cancer cells need the enzyme telomerase to elongate their telomeres (chromatin endings). In this study, we developed a gene regulatory model based on large data from transcription profiles from prostate cancer and chromatin-immuno-precipitation studies. We identified the developmental regulator PITX1 regulating telomerase. Besides observing experimental evidence of PITX1′s functional role in telomerase regulation, we also found PITX1 serving as a prognostic marker, as concluded from an analysis of more than 15,000 prostate cancer samples.

Abstract

The current risk stratification in prostate cancer (PCa) is frequently insufficient to adequately predict disease development and outcome. One hallmark of cancer is telomere maintenance. For telomere maintenance, PCa cells exclusively employ telomerase, making it essential for this cancer entity. However, TERT, the catalytic protein component of the reverse transcriptase telomerase, itself does not suit as a prognostic marker for prostate cancer as it is rather low expressed. We investigated if, instead of TERT, transcription factors regulating TERT may suit as prognostic markers. To identify transcription factors regulating TERT, we developed and applied a new gene regulatory modeling strategy to a comprehensive transcriptome dataset of 445 primary PCa. Six transcription factors were predicted as TERT regulators, and most prominently, the developmental morphogenic factor PITX1. PITX1 expression positively correlated with telomere staining intensity in PCa tumor samples. Functional assays and chromatin immune-precipitation showed that PITX1 activates TERT expression in PCa cells. Clinically, we observed that PITX1 is an excellent prognostic marker, as concluded from an analysis of more than 15,000 PCa samples. PITX1 expression in tumor samples associated with (i) increased Ki67 expression indicating increased tumor growth, (ii) a worse prognosis, and (iii) correlated with telomere length.

Gene regulatory circuitry of plant-environment interactions: scaling from cells to the field

Plant growth and development is the product of layers of sensing and regulation that are modulated by multifactorial environmental cues. Innovations in genomics currently allow gene regulatory control to be quantified at multiple scales and high resolution in defined cell populations and even in individual cells or nuclei in plants. The application of these 'omic technologies in highly controlled, as well as field environments is revolutionizing the recognition of factors critical to spatial and temporal responses to single or multiple environmental cues. Within and pan-species comparisons illuminate deeply conserved circuitry and targets of selection. This knowledge can benefit the breeding and engineering of crops with greater resilience to climate variability and the ability to augment nutrition through plant-microbial interactions.

Profiling transcription factor activity dynamics using intronic reads in time-series transcriptome data

Activities of transcription factors (TFs) are temporally modulated to regulate dynamic cellular processes, including development, homeostasis, and disease. Recent developments of bioinformatic tools have enabled the analysis of TF activities using transcriptome data. However, because these methods typically use exon-based target expression levels, the estimated TF activities have limited temporal accuracy. To address this, we proposed a TF activity measure based on intron-level information in time-series RNA-seq data, and implemented it to decode the temporal control of TF activities during dynamic processes. We showed that TF activities inferred from intronic reads can better recapitulate instantaneous TF activities compared to the exon-based measure. By analyzing public and our own time-series transcriptome data, we found that intron-based TF activities improve the characterization of temporal phasing of cycling TFs during circadian rhythm, and facilitate the discovery of two temporally opposing TF modules during T cell activation. Collectively, we anticipate that the proposed approach would be broadly applicable for decoding global transcriptional architecture during dynamic processes.

Taiji-reprogram: a framework to uncover cell-type specific regulators and predict cellular reprogramming cocktails

Cellular reprogramming is a promising technology to develop disease models and cell-based therapies. Identification of the key regulators defining the cell type specificity is pivotal to devising reprogramming cocktails for successful cell conversion but remains a great challenge. Here, we present a systems biology approach called Taiji-reprogram to efficiently uncover transcription factor (TF) combinations for conversion between 154 diverse cell types or tissues. This method integrates the transcriptomic and epigenomic data to construct cell-type specific genetic networks and assess the global importance of TFs in the network. Comparative analysis across cell types revealed TFs that are specifically important in a particular cell type and often tightly associated with cell-type specific functions. A systematic search of TFs with differential importance in the source and target cell types uncovered TF combinations for desired cell conversion. We have shown that Taiji-reprogram outperformed the existing methods to better recover the TFs in the experimentally validated reprogramming cocktails. This work not only provides a comprehensive catalog of TFs defining cell specialization but also suggests TF combinations for direct cell conversion.

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.

Modeling transcriptional regulation using gene regulatory networks based on multi-omics data sources

Background

Transcriptional regulation is complex, requiring multiple cis (local) and trans acting mechanisms working in concert to drive gene expression, with disruption of these processes linked to multiple diseases. Previous computational attempts to understand the influence of regulatory mechanisms on gene expression have used prediction models containing input features derived from cis regulatory factors. However, local chromatin looping and trans-acting mechanisms are known to also influence transcriptional regulation, and their inclusion may improve model accuracy and interpretation. In this study, we create a general model of transcription factor influence on gene expression by incorporating both cis and trans gene regulatory features.

Results

We describe a computational framework to model gene expression for GM12878 and K562 cell lines. This framework weights the impact of transcription factor-based regulatory data using multi-omics gene regulatory networks to account for both cis and trans acting mechanisms, and measures of the local chromatin context. These prediction models perform significantly better compared to models containing cis-regulatory features alone. Models that additionally integrate long distance chromatin interactions (or chromatin looping) between distal transcription factor binding regions and gene promoters also show improved accuracy. As a demonstration of their utility, effect estimates from these models were used to weight cis-regulatory rare variants for sequence kernel association test analyses of gene expression.

Conclusions

Our models generate refined effect estimates for the influence of individual transcription factors on gene expression, allowing characterization of their roles across the genome. This work also provides a framework for integrating multiple data types into a single model of transcriptional regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04126-3.

Derivation and Application of Molecular Signatures to Prostate Cancer: Opportunities and Challenges

Prostate cancer is a high-incidence cancer that requires improved patient stratification to ensure accurate predictions of risk and treatment response. Due to the significant contributions of transcription factors and epigenetic regulators to prostate cancer progression, there has been considerable progress made in developing gene signatures that may achieve this. Some of these are aligned to activities of key drivers such as the androgen receptor, whilst others are more agnostic. In this review, we present an overview of these signatures, the strategies for their derivation, and future perspectives on their continued development and evolution.

Mathematical Programming for Modeling Expression of a Gene Using Gurobi Optimizer to Identify Its Transcriptional Regulators

The cell expresses various genes in specific contexts with respect to internal and external perturbations to invoke appropriate responses. Transcription factors (TFs) orchestrate and define the expression level of genes by binding to their regulatory regions. Dysregulated expression of TFs often leads to aberrant expression changes of their target genes and is responsible for several diseases including cancers. In the last two decades, several studies experimentally identified target genes of several TFs. However, these studies are limited to a small fraction of the total TFs encoded by an organism, and only for those amenable to experimental settings. Experimental limitations lead to many computational techniques having been proposed to predict target genes of TFs. Linear modeling of gene expression is one of the most promising computational approaches, readily applicable to the thousands of expression datasets available in the public domain across diverse phenotypes. Linear models assume that the expression of a gene is the sum of expression of TFs regulating it. In this chapter, I introduce mathematical programming for the linear modeling of gene expression, which has certain advantages over the conventional statistical modeling approaches. It is fast, scalable to genome level and most importantly, allows mixed integer programming to tune the model outcome with prior knowledge on gene regulation.

Decoding myofibroblast origins in human kidney fibrosis

Kidney fibrosis is the hallmark of chronic kidney disease progression; however, at present no antifibrotic therapies exist^1-3. The origin, functional heterogeneity and regulation of scar-forming cells that occur during human kidney fibrosis remain poorly understood^1,2,4. Here, using single-cell RNA sequencing, we profiled the transcriptomes of cells from the proximal and non-proximal tubules of healthy and fibrotic human kidneys to map the entire human kidney. This analysis enabled us to map all matrix-producing cells at high resolution, and to identify distinct subpopulations of pericytes and fibroblasts as the main cellular sources of scar-forming myofibroblasts during human kidney fibrosis. We used genetic fate-tracing, time-course single-cell RNA sequencing and ATAC-seq (assay for transposase-accessible chromatin using sequencing) experiments in mice, and spatial transcriptomics in human kidney fibrosis, to shed light on the cellular origins and differentiation of human kidney myofibroblasts and their precursors at high resolution. Finally, we used this strategy to detect potential therapeutic targets, and identified NKD2 as a myofibroblast-specific target in human kidney fibrosis.

Determinants of correlated expression of transcription factors and their target genes

While transcription factors (TFs) are known to regulate the expression of their target genes (TGs), only a weak correlation of expression between TFs and their TGs has generally been observed. As lack of correlation could be caused by additional layers of regulation, the overall correlation distribution may hide the presence of a subset of regulatory TF-TG pairs with tight expression coupling. Using reported regulatory pairs in the plant Arabidopsis thaliana along with comprehensive gene expression information and testing a wide array of molecular features, we aimed to discern the molecular determinants of high expression correlation of TFs and their TGs. TF-family assignment, stress-response process involvement, short genomic distances of the TF-binding sites to the transcription start site of their TGs, few required protein-protein-interaction connections to establish physical interactions between the TF and polymerase-II, unambiguous TF-binding motifs, increased numbers of miRNA target-sites in TF-mRNAs, and a young evolutionary age of TGs were found particularly indicative of high TF-TG correlation. The modulating roles of post-transcriptional, post-translational processes, and epigenetic factors have been characterized as well. Our study reveals that regulatory pairs with high expression coupling are associated with specific molecular determinants.

DrivAER: Identification of driving transcriptional programs in single-cell RNA sequencing data

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) unfolds complex transcriptomic datasets into detailed cellular maps. Despite recent success, there is a pressing need for specialized methods tailored towards the functional interpretation of these cellular maps.

FINDINGS

Here, we present DrivAER, a machine learning approach for the identification of driving transcriptional programs using autoencoder-based relevance scores. DrivAER scores annotated gene sets on the basis of their relevance to user-specified outcomes such as pseudotemporal ordering or disease status. DrivAER iteratively evaluates the information content of each gene set with respect to the outcome variable using autoencoders. We benchmark our method using extensive simulation analysis as well as comparison to existing methods for functional interpretation of scRNA-seq data. Furthermore, we demonstrate that DrivAER extracts key pathways and transcription factors that regulate complex biological processes from scRNA-seq data.

CONCLUSIONS

By quantifying the relevance of annotated gene sets with respect to specified outcome variables, DrivAER greatly enhances our ability to understand the underlying molecular mechanisms.

Prediction of single-cell gene expression for transcription factor analysis

BACKGROUND

Single-cell RNA sequencing is a powerful technology to discover new cell types and study biological processes in complex biological samples. A current challenge is to predict transcription factor (TF) regulation from single-cell RNA data.

RESULTS

Here, we propose a novel approach for predicting gene expression at the single-cell level using cis-regulatory motifs, as well as epigenetic features. We designed a tree-guided multi-task learning framework that considers each cell as a task. Through this framework we were able to explain the single-cell gene expression values using either TF binding affinities or TF ChIP-seq data measured at specific genomic regions. TFs identified using these models could be validated by the literature.

CONCLUSION

Our proposed method allows us to identify distinct TFs that show cell type-specific regulation. This approach is not limited to TFs but can use any type of data that can potentially be used in explaining gene expression at the single-cell level to study factors that drive differentiation or show abnormal regulation in disease. The implementation of our workflow can be accessed under an MIT license via https://github.com/SchulzLab/Triangulate.

Identification of putative master regulators in rheumatoid arthritis synovial fibroblasts using gene expression data and network inference

Rheumatoid arthritis (RA) is a systemic autoimmune disease that affects the synovial joints of the body. Rheumatoid arthritis fibroblast-like synoviocytes (RA FLS) are central players in the disease pathogenesis, as they are involved in the secretion of cytokines and proteolytic enzymes, exhibit invasive traits, high rate of self-proliferation and an apoptosis-resistant phenotype. We aim at characterizing transcription factors (TFs) that are master regulators in RA FLS and could potentially explain phenotypic traits. We make use of differentially expressed genes in synovial tissue from patients suffering from RA and osteoarthritis (OA) to infer a TF co-regulatory network, using dedicated software. The co-regulatory network serves as a reference to analyze microarray and single-cell RNA-seq data from isolated RA FLS. We identified five master regulators specific to RA FLS, namely BATF, POU2AF1, STAT1, LEF1 and IRF4. TF activity of the identified master regulators was also estimated with the use of two additional, independent software. The identified TFs contribute to the regulation of inflammation, proliferation and apoptosis, as indicated by the comparison of their differentially expressed target genes with hallmark molecular signatures derived from the Molecular Signatures Database (MSigDB). Our results show that TFs influence could be used to identify putative master regulators of phenotypic traits and suggest novel, druggable targets for experimental validation.

Significant association between host transcriptome-derived HPV oncogene E6* influence score and carcinogenic pathways, tumor size, and survival in head and neck cancer

BACKGROUND

Human papillomavirus (HPV) oncogenes E6, E7, and shorter isoforms of E6 (E6*) are known carcinogenic factors in head and neck squamous cell carcinoma (HNSCC). Little is known regarding E6* functions.

METHODS

We analyzed RNA-seq data from 68 HNSCC HPV type 16-positive tumors to determine host genes and pathways associated with E6+E7 expression (E6E7) or the percent of full-length E6 (E6%FL). Influence scores of E6E7 and E6%FL were used to test for associations with clinical variables.

RESULTS

For E6E7, we recapitulated all major known affected pathways and revealed additional pathways. E6%FL was found to affect mitochondrial processes, and E6%FL influence score was significantly associated with overall survival and tumor size.

CONCLUSIONS

HPV E6E7 and E6* result in extensive, dose-dependent compensatory effects and dysregulation of key cancer pathways. The switch from E6 to E6* promotes oxidative phosphorylation, larger tumor size, and worse prognosis, potentially serving as a prognostic factor for HPV-positive HNSCC.

Chimeric Activators and Repressors Define HY5 Activity and Reveal a Light-Regulated Feedback Mechanism

The first exposure to light marks a crucial transition in plant development. This transition relies on the transcription factor HY5 controlling a complex downstream growth program. Despite its importance, its function in transcription remains unclear. Previous studies have generated lists of thousands of potential target genes and competing models of HY5 transcription regulation. In this work, we carry out detailed phenotypic and molecular analysis of constitutive activator and repressor HY5 fusion proteins. Using this strategy, we were able to filter out large numbers of genes that are unlikely to be direct targets, allowing us to eliminate several proposed models of HY5's mechanism of action. We demonstrate that the primary activity of HY5 is promoting transcription and that this function relies on other, likely light-regulated, factors. In addition, this approach reveals a molecular feedback loop via the COP1/SPA E3 ubiquitin ligase complex, suggesting a mechanism that maintains low HY5 in the dark, primed for rapid accumulation to reprogram growth upon light exposure. Our strategy is broadly adaptable to the study of transcription factor activity. Lastly, we show that modulating this feedback loop can generate significant phenotypic diversity in both Arabidopsis (Arabidopsis thaliana) and tomato (Solanum lycopersicum).

Inferring TF activation order in time series scRNA-Seq studies

Methods for the analysis of time series single cell expression data (scRNA-Seq) either do not utilize information about transcription factors (TFs) and their targets or only study these as a post-processing step. Using such information can both, improve the accuracy of the reconstructed model and cell assignments, while at the same time provide information on how and when the process is regulated. We developed the Continuous-State Hidden Markov Models TF (CSHMM-TF) method which integrates probabilistic modeling of scRNA-Seq data with the ability to assign TFs to specific activation points in the model. TFs are assumed to influence the emission probabilities for cells assigned to later time points allowing us to identify not just the TFs controlling each path but also their order of activation. We tested CSHMM-TF on several mouse and human datasets. As we show, the method was able to identify known and novel TFs for all processes, assigned time of activation agrees with both expression information and prior knowledge and combinatorial predictions are supported by known interactions. We also show that CSHMM-TF improves upon prior methods that do not utilize TF-gene interaction.

An important attribute of time series single cell RNA-Seq (scRNA-Seq) data, is the ability to infer continuous trajectories of genes based on orderings of the cells. While several methods have been developed for ordering cells and inferring such trajectories, to date it was not possible to use these to infer the temporal activity of several key TFs. These TFs are are only post-transcriptionally regulated and so their expression does not provide complete information on their activity. To address this we developed the Continuous-State Hidden Markov Models TF (CSHMM-TF) methods that assigns continuous activation time to TFs based on both, their expression and the expression of their targets. Applying our method to several time series scRNA-Seq datasets we show that it correctly identifies the key regulators for the processes being studied. We analyze the temporal assignments for these TFs and show that they provide new insights about combinatorial regulation and the ordering of TF activation. We used several complementary sources to validate some of these predictions and discuss a number of other novel suggestions based on the method. As we show, the method is able to scale to large and noisy datasets and so is appropriate for several studies utilizing time series scRNA-Seq data.

Reprogramming of macrophages employing gene regulatory and metabolic network models

Upon exposure to different stimuli, resting macrophages undergo classical or alternative polarization into distinct phenotypes that can cause fatal dysfunction in a large range of diseases, such as systemic infection leading to sepsis or the generation of an immunosuppressive tumor microenvironment. Investigating gene regulatory and metabolic networks, we observed two metabolic switches during polarization. Most prominently, anaerobic glycolysis was utilized by M1-polarized macrophages, while the biosynthesis of inosine monophosphate was upregulated in M2-polarized macrophages. Moreover, we observed a switch in the urea cycle. Gene regulatory network models revealed E2F1, MYC, PPARγ and STAT6 to be the major players in the distinct signatures of these polarization events. Employing functional assays targeting these regulators, we observed the repolarization of M2-like cells into M1-like cells, as evidenced by their specific gene expression signatures and cytokine secretion profiles. The predicted regulators are essential to maintaining the M2-like phenotype and function and thus represent potential targets for the therapeutic reprogramming of immunosuppressive M2-like macrophages.

Modelling TERT regulation across 19 different cancer types based on the MIPRIP 2.0 gene regulatory network approach

Background

Reactivation of the telomerase reverse transcriptase gene TERT is a central feature for unlimited proliferation of the majority of cancers. However, the underlying regulatory processes are only partly understood.

Results

We assembled regulator binding information from serveral sources to construct a generic human and mouse gene regulatory network. Advancing our “Mixed Integer linear Programming based Regulatory Interaction Predictor” (MIPRIP) approach, we identified the most common and cancer-type specific regulators of TERT across 19 different human cancers. The results were validated by using the well-known TERT regulation by the ETS1 transcription factor in a subset of melanomas with mutations in the TERT promoter.

Our improved MIPRIP2 R-package and the associated generic regulatory networks are freely available at https://github.com/KoenigLabNM/MIPRIP.

Conclusion

MIPRIP 2.0 identified common as well as tumor type specific regulators of TERT. The software can be easily applied to transcriptome datasets to predict gene regulation for any gene and disease/condition under investigation.

Cancer classification from time series microarray data through regulatory Dynamic Bayesian Networks

Genomic profiling of cancer studies has generated comprehensive gene expression patterns for diverse phenotypes. Computational methods which employ transcriptomics datasets have been proposed to model gene expression data. Dynamic Bayesian Networks (DBNs) have been used for modeling time series datasets and for the inference of regulatory networks. Furthermore, cancer classification through DBN-based approaches could reveal the importance of exploiting knowledge from statistically significant genes and key regulatory molecules. Although microarray datasets have been employed extensively by several classification methods for decision making, the use of new knowledge from the pathway level has not been addressed adequately in the literature in terms of DBNs for cancer classification. In the present study, we identify the genes that act as regulators and mediate the activity of transcription factors that have been found in all promoters of our differentially expressed gene sets. These features serve as potential priors for distinguishing tumor from normal samples using a DBN-based classification approach. We employed three microarray datasets from the Gene Expression Omnibus (GEO) public functional repository and performed differential expression analysis. Promoter and pathway analysis of the identified genes revealed the key regulators which influence the transcription mechanisms of these genes. We applied the DBN algorithm on selected genes and identified the features that can accurately classify the samples into tumors and controls. Both accuracy and Area Under the Curve (AUC) were high for the gene sets comprising of the differentially expressed genes along with their master regulators (accuracy: 70.8%-98.5%; AUC: 0.562-0.985).

Melanoma plasticity and phenotypic diversity: therapeutic barriers and opportunities

An incomplete view of the mechanisms that drive metastasis, the primary cause of cancer-related death, has been a major barrier to development of effective therapeutics and prognostic diagnostics. Increasing evidence indicates that the interplay between microenvironment, genetic lesions, and cellular plasticity drives the metastatic cascade and resistance to therapies. Here, using melanoma as a model, we outline the diversity and trajectories of cell states during metastatic dissemination and therapy exposure, and highlight how understanding the magnitude and dynamics of nongenetic reprogramming in space and time at single-cell resolution can be exploited to develop therapeutic strategies that capitalize on nongenetic tumor evolution.

Benchmark and integration of resources for the estimation of human transcription factor activities

The prediction of transcription factor (TF) activities from the gene expression of their targets (i.e., TF regulon) is becoming a widely used approach to characterize the functional status of transcriptional regulatory circuits. Several strategies and data sets have been proposed to link the target genes likely regulated by a TF, each one providing a different level of evidence. The most established ones are (1) manually curated repositories, (2) interactions derived from ChIP-seq binding data, (3) in silico prediction of TF binding on gene promoters, and (4) reverse-engineered regulons from large gene expression data sets. However, it is not known how these different sources of regulons affect the TF activity estimations and, thereby, downstream analysis and interpretation. Here we compared the accuracy and biases of these strategies to define human TF regulons by means of their ability to predict changes in TF activities in three reference benchmark data sets. We assembled a collection of TF-target interactions for 1541 human TFs and evaluated how different molecular and regulatory properties of the TFs, such as the DNA-binding domain, specificities, or mode of interaction with the chromatin, affect the predictions of TF activity. We assessed their coverage and found little overlap on the regulons derived from each strategy and better performance by literature-curated information followed by ChIP-seq data. We provide an integrated resource of all TF-target interactions derived through these strategies, with confidence scores, as a resource for enhanced prediction of TF activities.

Pan-Cancer analysis of the expression and regulation of matrisome genes across 32 tumor types

The microenvironment plays a central role in cancer, and neoplastic cells actively shape it to their needs by complex arrays of extracellular matrix (ECM) proteins, enzymes, cytokines and growth factors collectively referred to as the matrisome. Studies on the cancer matrisome have been performed for single or few neoplasms, but a more systematic analysis is still missing. Here we present a Pan-Cancer study of matrisome gene expression in 10,487 patients across 32 tumor types, supplemented with transcription factors (TFs) and driver genes/pathways regulating each tumor's matrisome. We report on 919 TF-target pairs, either used specifically or shared across tumor types, and their prognostic significance, 40 master regulators, 31 overarching regulatory pathways and the potential for druggability with FDA-approved cancer drugs. These results provide a comprehensive transcriptional architecture of the cancer matrisome and suggest the need for development of specific matrisome-targeting approaches for future therapies.

•

In-depth characterization of matrisome gene expression and regulation in 10,487 patients across 32 human tumor types.

•

Identification of transcription factor (TF) and “master regulators” governing each cancer’s matrisome.

•

Analysis unveils therapeutic possibilities and suggests new treatments by repurposing of FDA-approved cancer drugs.

Taiji: System-level identification of key transcription factors reveals transcriptional waves in mouse embryonic development

Transcriptional regulation is pivotal to the specification of distinct cell types during embryonic development. However, it still lacks a systematic way to identify key transcription factors (TFs) orchestrating the temporal and tissue specificity of gene expression. Here, we integrated epigenomic and transcriptomic data to reveal key regulators from two cells to postnatal day 0 in mouse embryogenesis. We predicted three-dimensional chromatin interactions in 12 tissues across eight developmental stages, which facilitates linking TFs to their target genes for constructing transcriptional regulatory networks. To identify driver TFs, we developed a new algorithm, dubbed Taiji, to assess the global influence of each TF and systematically uncovered TFs critical for lineage-specific and stage-dependent tissue specification. We have also identified TF combinations that function in spatiotemporal order to form transcriptional waves regulating developmental progress. Furthermore, lacking stage-specific TF combinations suggests a distributed timing strategy to orchestrate the coordination between tissues during embryonic development.

Global mapping of transcription factor motifs in human aging

Biological aging is a complex process dependent on the interplay of cell autonomous and tissue contextual changes which occur in response to cumulative molecular stress and manifest through adaptive transcriptional reprogramming. Here we describe a transcription factor (TF) meta-analysis of gene expression datasets accrued from 18 tissue sites collected at different biological ages and from 7 different in-vitro aging models. In-vitro aging platforms included replicative senescence and an energy restriction model in quiescence (ERiQ), in which ATP was transiently reduced. TF motifs in promoter regions of trimmed sets of target genes were scanned using JASPAR and TRANSFAC. TF signatures established a global mapping of agglomerating motifs with distinct clusters when ranked hierarchically. Remarkably, the ERiQ profile was shared with the majority of in-vivo aged tissues. Fitting motifs in a minimalistic protein-protein network allowed to probe for connectivity to distinct stress sensors. The DNA damage sensors ATM and ATR linked to the subnetwork associated with senescence. By contrast, the energy sensors PTEN and AMPK connected to the nodes in the ERiQ subnetwork. These data suggest that metabolic dysfunction may be linked to transcriptional patterns characteristic of many aged tissues and distinct from cumulative DNA damage associated with senescence.

Systems-epigenomics inference of transcription factor activity implicates aryl-hydrocarbon-receptor inactivation as a key event in lung cancer development

Background

Diverse molecular alterations associated with smoking in normal and precursor lung cancer cells have been reported, yet their role in lung cancer etiology remains unclear. A prominent example is hypomethylation of the aryl hydrocarbon-receptor repressor (AHRR) locus, which is observed in blood and squamous epithelial cells of smokers, but not in lung cancer.

Results

Using a novel systems-epigenomics algorithm, called SEPIRA, which leverages the power of a large RNA-sequencing expression compendium to infer regulatory activity from messenger RNA expression or DNA methylation (DNAm) profiles, we infer the landscape of binding activity of lung-specific transcription factors (TFs) in lung carcinogenesis. We show that lung-specific TFs become preferentially inactivated in lung cancer and precursor lung cancer lesions and further demonstrate that these results can be derived using only DNAm data. We identify subsets of TFs which become inactivated in precursor cells. Among these regulatory factors, we identify AHR, the aryl hydrocarbon-receptor which controls a healthy immune response in the lung epithelium and whose repressor, AHRR, has recently been implicated in smoking-mediated lung cancer. In addition, we identify FOXJ1, a TF which promotes growth of airway cilia and effective clearance of the lung airway epithelium from carcinogens.

Conclusions

We identify TFs, such as AHR, which become inactivated in the earliest stages of lung cancer and which, unlike AHRR hypomethylation, are also inactivated in lung cancer itself. The novel systems-epigenomics algorithm SEPIRA will be useful to the wider epigenome-wide association study community as a means of inferring regulatory activity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1366-0) contains supplementary material, which is available to authorized users.

Roles of alternative splicing in modulating transcriptional regulation

Background

The ability of a transcription factor to regulate its targets is modulated by a variety of genetic and epigenetic mechanisms. Alternative splicing can modulate gene function by adding or removing certain protein domains, and therefore affect the activity of protein. Reverse engineering of gene regulatory networks using gene expression profiles has proven valuable in dissecting the logical relationships among multiple proteins during the transcriptional regulation. However, it is unclear whether alternative splicing of certain proteins affects the activity of other transcription factors.

Results

In order to investigate the roles of alternative splicing during transcriptional regulation, we constructed a statistical model to infer whether the alternative splicing events of modulator proteins can affect the ability of key transcription factors in regulating the expression levels of their transcriptional targets. We tested our strategy in KIRC (Kidney Renal Clear Cell Carcinoma) using the RNA-seq data downloaded from TCGA (the Cancer Genomic Atlas). We identified 828of modulation relationships between the splicing levels of modulator proteins and activity levels of transcription factors. For instance, we found that the activity levels of GR (glucocorticoid receptor) protein, a key transcription factor in kidney, can be influenced by the splicing status of multiple proteins, including TP53, MDM2 (mouse double minute 2 homolog), RBM14 (RNA-binding protein 14) and SLK (STE20 like kinase). The influenced GR-targets are enriched by key cancer-related pathways, including p53 signaling pathway, TR/RXR activation, CAR/RXR activation, G1/S checkpoint regulation pathway, and G2/M DNA damage checkpoint regulation pathway.

Conclusions

Our analysis suggests, for the first time, that exon inclusion levels of certain regulatory proteins can affect the activities of many transcription factors. Such analysis can potentially unravel a novel mechanism of how splicing variation influences the cellular function and provide important insights for how dysregulation of splicing outcome can lead to various diseases.

Electronic supplementary material

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

Phosphoproteomics-Based Profiling of Kinase Activities in Cancer Cells

Cellular signaling, predominantly mediated by phosphorylation through protein kinases, is found to be deregulated in most cancers. Accordingly, protein kinases have been subject to intense investigations in cancer research, to understand their role in oncogenesis and to discover new therapeutic targets. Despite great advances, an understanding of kinase dysfunction in cancer is far from complete.A powerful tool to investigate phosphorylation is mass-spectrometry (MS)-based phosphoproteomics, which enables the identification of thousands of phosphorylated peptides in a single experiment. Since every phosphorylation event results from the activity of a protein kinase, high-coverage phosphoproteomics data should indirectly contain comprehensive information about the activity of protein kinases.In this chapter, we discuss the use of computational methods to predict kinase activity scores from MS-based phosphoproteomics data. We start with a short explanation of the fundamental features of the phosphoproteomics data acquisition process from the perspective of the computational analysis. Next, we briefly review the existing databases with experimentally verified kinase-substrate relationships and present a set of bioinformatic tools to discover novel kinase targets. We then introduce different methods to infer kinase activities from phosphoproteomics data and these kinase-substrate relationships. We illustrate their application with a detailed protocol of one of the methods, KSEA (Kinase Substrate Enrichment Analysis). This method is implemented in Python within the framework of the open-source Kinase Activity Toolbox (kinact), which is freely available at http://github.com/saezlab/kinact/ .

Fungal biomarker discovery by integration of classifiers

Background

The human immune system is responsible for protecting the host from infection. However, in immunocompromised individuals the risk of infection increases substantially with possible drastic consequences. In extreme, systemic infection can lead to sepsis which is responsible for innumerous deaths worldwide. Amongst its causes are infections by bacteria and fungi. To increase survival, it is mandatory to identify the type of infection rapidly. Discriminating between fungal and bacterial pathogens is key to determine if antifungals or antibiotics should be administered, respectively. For this, in situ experiments have been performed to determine regulation mechanisms of the human immune system to identify biomarkers. However, these studies led to heterogeneous results either due different laboratory settings, pathogen strains, cell types and tissues, as well as the time of sample extraction, to name a few.

Methods

To generate a gene signature capable of discriminating between fungal and bacterial infected samples, we employed Mixed Integer Linear Programming (MILP) based classifiers on several datasets comprised of the above mentioned pathogens.

Results

When combining the classifiers by a joint optimization we could increase the consistency of the biomarker gene list independently of the experimental setup. An increase in pairwise overlap (the number of genes that overlap in each cross-validation) of 43% was obtained by this approach when compared to that of single classifiers. The refined gene list was composed of 19 genes and ranked according to consistency in expression (up- or down-regulated) and most of them were linked either directly or indirectly to the ERK-MAPK signalling pathway, which has been shown to play a key role in the immune response to infection. Testing of the identified 12 genes on an unseen dataset yielded an average accuracy of 83%.

Conclusions

In conclusion, our method allowed the combination of independent classifiers and increased consistency and reliability of the generated gene signatures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4006-x) contains supplementary material, which is available to authorized users.

Estimating genome-wide regulatory activity from multi-omics data sets using mathematical optimization

Background

Gene regulation is one of the most important cellular processes, indispensable for the adaptability of organisms and closely interlinked with several classes of pathogenesis and their progression. Elucidation of regulatory mechanisms can be approached by a multitude of experimental methods, yet integration of the resulting heterogeneous, large, and noisy data sets into comprehensive and tissue or disease-specific cellular models requires rigorous computational methods. Recently, several algorithms have been proposed which model genome-wide gene regulation as sets of (linear) equations over the activity and relationships of transcription factors, genes and other factors. Subsequent optimization finds those parameters that minimize the divergence of predicted and measured expression intensities. In various settings, these methods produced promising results in terms of estimating transcription factor activity and identifying key biomarkers for specific phenotypes. However, despite their common root in mathematical optimization, they vastly differ in the types of experimental data being integrated, the background knowledge necessary for their application, the granularity of their regulatory model, the concrete paradigm used for solving the optimization problem and the data sets used for evaluation.

Results

Here, we review five recent methods of this class in detail and compare them with respect to several key properties. Furthermore, we quantitatively compare the results of four of the presented methods based on publicly available data sets.

Conclusions

The results show that all methods seem to find biologically relevant information. However, we also observe that the mutual result overlaps are very low, which contradicts biological intuition. Our aim is to raise further awareness of the power of these methods, yet also to identify common shortcomings and necessary extensions enabling focused research on the critical points.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0419-z) contains supplementary material, which is available to authorized users.

A novel statistical approach for identification of the master regulator transcription factor

BACKGROUND

Transcription factors are known to play key roles in carcinogenesis and therefore, are gaining popularity as potential therapeutic targets in drug development. A 'master regulator' transcription factor often appears to control most of the regulatory activities of the other transcription factors and the associated genes. This 'master regulator' transcription factor is at the top of the hierarchy of the transcriptomic regulation. Therefore, it is important to identify and target the master regulator transcription factor for proper understanding of the associated disease process and identifying the best therapeutic option.

METHODS

We present a novel two-step computational approach for identification of master regulator transcription factor in a genome. At the first step of our method we test whether there exists any master regulator transcription factor in the system. We evaluate the concordance of two ranked lists of transcription factors using a statistical measure. In case the concordance measure is statistically significant, we conclude that there is a master regulator. At the second step, our method identifies the master regulator transcription factor, if there exists one.

RESULTS

In the simulation scenario, our method performs reasonably well in validating the existence of a master regulator when the number of subjects in each treatment group is reasonably large. In application to two real datasets, our method ensures the existence of master regulators and identifies biologically meaningful master regulators. An R code for implementing our method in a sample test data can be found in http://www.somnathdatta.org/software .

CONCLUSION

We have developed a screening method of identifying the 'master regulator' transcription factor just using only the gene expression data. Understanding the regulatory structure and finding the master regulator help narrowing the search space for identifying biomarkers for complex diseases such as cancer. In addition to identifying the master regulator our method provides an overview of the regulatory structure of the transcription factors which control the global gene expression profiles and consequently the cell functioning.

Systematic Analysis of Transcriptional and Post-transcriptional Regulation of Metabolism in Yeast

Cells react to extracellular perturbations with complex and intertwined responses. Systematic identification of the regulatory mechanisms that control these responses is still a challenge and requires tailored analyses integrating different types of molecular data. Here we acquired time-resolved metabolomics measurements in yeast under salt and pheromone stimulation and developed a machine learning approach to explore regulatory associations between metabolism and signal transduction. Existing phosphoproteomics measurements under the same conditions and kinase-substrate regulatory interactions were used to in silico estimate the enzymatic activity of signalling kinases. Our approach identified informative associations between kinases and metabolic enzymes capable of predicting metabolic changes. We extended our analysis to two studies containing transcriptomics, phosphoproteomics and metabolomics measurements across a comprehensive panel of kinases/phosphatases knockouts and time-resolved perturbations to the nitrogen metabolism. Changes in activity of transcription factors, kinases and phosphatases were estimated in silico and these were capable of building predictive models to infer the metabolic adaptations of previously unseen conditions across different dynamic experiments. Time-resolved experiments were significantly more informative than genetic perturbations to infer metabolic adaptation. This difference may be due to the indirect nature of the associations and of general cellular states that can hinder the identification of causal relationships. This work provides a novel genome-scale integrative analysis to propose putative transcriptional and post-translational regulatory mechanisms of metabolic processes.

Phosphorylation is a broad regulatory mechanism with implications in nearly all processes of the cell. However, a global understanding of possible regulatory mechanisms remains elusive. In this study, we examined the potential regulatory role of kinases, phosphatases and transcription-factors in yeast metabolism across a variety of steady-state and dynamic conditions. The main novelty of our analysis was to infer putative regulatory interactions from in silico estimated activity of transcription-factors and kinases/phosphatases. This provided functional information about the proteins important for the experimental conditions at hand that had not been uncovered before. We showed that activity profiles are predictive features to estimate metabolite changes in dynamic experiments, while the same was not visible in steady-state conditions. We also showed that dynamic experiments could be used to recapitulate and provide novel TFs-metabolite and K/Ps-metabolite regulatory associations. We believe these findings illustrates the usefulness of this approach for future integrative studies interested in studying metabolic regulation.